Advances in Mathematical Physics

Volume 2014 (2014), Article ID 641918, 6 pages

http://dx.doi.org/10.1155/2014/641918

Research Article

Lie Group Method of the Diffusion Equations

Jian-Qiang Sun,¹ Rong-Fang Huang,¹ Xiao-Yan Gu,² and Ling Yu¹

¹Department of Mathematics, College of Information Science and Technology, Hainan University, Haikou 570228, China

²Department of Physics, East China University of Science and Technology, Shanghai 200237, China

Received 1 January 2014; Accepted 7 April 2014; Published 29 April 2014

Academic Editor: Shi-Hai Dong

Copyright © 2014 Jian-Qiang Sun et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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 [1–12]. 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 [15–17]. 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 [15–17].

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.

Table 1: The stability comparison of ERK method and EEI method.

ERK
EEI

0.028

0.007

0.0016

0.0005

Table 2: The errors of the numerical solutions and the exact solutions at different times by ERK method and EEI method.

ERK
EEI

(a)

(b)

Figure 1: The numerical solutions of Burgers’ equation at obtained by ERK method and EEI method.

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.

Table 3: The stability comparison of ERK method and EEI method.

ERK
EEI

Table 4: The errors of the numerical solutions at different times obtained by ERK method and EEI method.

ERK
EEI
Err

(a)

(b)

Figure 2: The numerical solutions of the Allen-Cahn equation at obtained by the two methods with and .

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).

References

H. Munthe-Kaas, “Runge-Kutta methods on Lie groups,” BIT Numerical Mathematics, vol. 38, no. 1, pp. 92–111, 1998.

A. Iserles, H. Munthe-Kaas, S. P. Nørsett, and A. Zanna, “Lie group methods,” Acta Numerica, vol. 9, pp. 215–365, 2000.

G. J. Cooper, “Stability of Runge-Kutta methods for trajectory problems,” IMA Journal of Numerical Analysis, vol. 7, no. 1, pp. 1–13, 1987.

P. E. Crouch and R. Grossman, “Numerical integration of ordinary differential equations on manifolds,” Journal of Nonlinear Science, vol. 3, no. 1, pp. 1–33, 1993.

A. Iserles, “Solving linear ordinary differential equations by exponentials of iterated commutators,” Numerische Mathematik, vol. 45, no. 2, pp. 183–199, 1984.

M. P. Calvo, A. Iserles, and A. Zanna, “Conservative methods for the Toda lattice equations,” IMA Journal of Numerical Analysis, vol. 19, no. 4, pp. 509–523, 1999.

S. M. Cox and P. C. Matthews, “Exponential time differencing for stiff systems,” Journal of Computational Physics, vol. 176, no. 2, pp. 430–455, 2002.

L. Q. Chen, C. Wolverton, V. Vaithyanathan, and Z. Liu, “Modeling solid-state phase transformations and microstructure evolution,” MRS Bulletion, vol. 26, pp. 197–292, 2001.

G. Akrivis, M. Crouzeix, and C. Makridakis, “Implicit-explicit multistep methods for quasilinear parabolic equations,” Numerische Mathematik, vol. 82, no. 4, pp. 521–541, 1999.

E. Weinan, “Analysis of the heterogeneous multiscale method for ordinary differential equations,” Communications in Mathematical Sciences, vol. 3, pp. 423–436, 2003.

A. Zanna, “Collocation and relaxed collocation for the Fer and the Magnus expansions,” SIAM Journal on Numerical Analysis, vol. 36, no. 4, pp. 1145–1182, 1999.

L. R. Petzold, L. O. Jay, and J. Yen, “Numerical solution of highly oscillatory ordinary differential equations,” Acta Numerica, vol. 6, pp. 437–483, 1997.

E. Hairer and G. Wanner, Solving Ordinary Differential Equations II: Stiff and Differential Algebraic Problems, Springer, New York, NY, USA, 1999.

M. Z. Qin and Y. S. Wang, Structure-Preserving Algorithm for Partial Differential Equation, Zhejiang Science and Technology Publishing, 2012.

R. C. Mittal and P. Singhal, “Numerical solution of Burger's equation,” Communications in Numerical Methods in Engineering, vol. 9, pp. 163–175, 1993.

S. Kutluay, A. R. Bahadir, and A. Özdeş, “Numerical solution of one-dimensional Burgers' equation: explicit and exact-explicit finite difference methods,” Journal of Computational and Applied Mathematics, vol. 103, no. 2, pp. 251–261, 1999.

E. R. Benton and G. W. Platzman, “A table of solution of the one-deimensional Burgers' equation,” Quarterly of Applied Mathematics, vol. 30, pp. 195–212, 1972.