Lie Group Method of the Diffusion Equations

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 MuntheKaas, Iserles, Nørsett, and their collaborators. In general, if a Lie groupG acts freely and transitively on amanifoldM, then a differential equation onM uniquely determines a differential equation on G. Thus, replacing the original differential on M by the equivalent differential equation on Lie algebra g 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-Kuttamethods. In Section 4, theAllen-Cahn equation is solved by the explicit exponential integrator method and the explicit Runge-Kutta methods. At last, some conclusions are obtained.


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][2][3][4][5][6][7][8][9][10][11][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.

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 V ∈  and  0 () ∈  the solution of at time  = 1 with the initial condition  0 () and some boundary conditions.There is freedom to choose the function , the corresponding differential equation in the Lie algebra 2

Advances in Mathematical Physics
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   = ∑  =1  , .The algorithm integrates (1) from  = 0 to  = ℎ: Set  0 = .
For  = 1, 2, . . ., End where ).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  +1 =   + ℎ, 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. exp −1  (  , ) is the  order truncation error of the  exp −1  (  ).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.

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 Re > 0 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][16][17].With the initial condition (, 0) = sin(), 0 ≤  ≤ 1, and  > 0, the exact solution of (18) is and   () is the Bessel function.Recently some numerical methods of Burgers' equation have been proposed [15][16][17].
In Table 1, we compare the stability property of ERK method and EEI methods with Re = 100.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  = 100 and ℎ = 0.002.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  ∈ [0, 3] are obtained by ERK method and the EEI method with  = 100, Re = 100, and ℎ = 0.005.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.
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  = −0.4 with ℎ = 0.01 obtained by ERK method and the numerical solution obtained by EEI method.The error is up to 10 −8 at different times; the error is (ℎ 4 ).We can conclude that they have the same accuracy.In Figure 2, the numerical solutions of the Allen-Cahn equation in  ∈ [0, 200] were obtained by ERK method and EEI method with  = 100, ℎ = 0.01, and  = −0.4. 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.

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.

Figure 2 :
Figure 2: The numerical solutions of the Allen-Cahn equation at  ∈ [0, 200] obtained by the two methods with  = 100 and  = 0.01.

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

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

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

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