A Shannon-Rugge-Kutta-Gill method for solving convection-diffusion equations is discussed. This approach transforms convection-diffusion equations into one-dimensional equations at collocations points, which we solve by Runge-Kutta-Gill method. A concrete example solved is used to examine the method’s feasibility.

1. Introduction

Most of the physics phenomenon are stated in terms of partial differential equations (PDEs). Convection-diffusion equation is a kind of PDE which can be used in many science and technology fields, especially in image and signal procession such as image segmentation and the quickly stability of image processing. The numerical solution of convection-diffusion equations as an important subject has always attracted the attentions of the researchers for a long time.

The standard Galerkin finite-element method can solve the solution of the equations, but it is numerically unstable for small values of the diffusion parameter. So on the basis of this method, in [1] King and Krueger investigated the effect of a stabilized finite-element approximation and drew a conclusion that the stabilized system can provide accurate controllers. A novel multilevel particle methods and two complementary approaches are researched in [2]. In this paper, a new class of particle based on mapping functions is introduced, and particle remeshing is used as a key element in overlapping domains in the particle-AMR method. For the two-dimensional convection-diffusion equation, Gupta et al. proposed a fourth-order nine-point compact finite-difference formulae, which is shown to be computationally efficient and stable and yield highly accurate numerical solutions in [3, 4]. The resulting linear system is solved by classical iterative methods for large values of the Reynolds number in [4]. With the wavelet method, Shi et al. solved the solution of covection-diffusion equations by Haar wavelet method in [5]. In short, convection-diffusion equations are studied by scholars via different methods.

Recently, Wavelet analysis as a new subject has attracted a lot of attention. As a mathematical tool, it has been widely used in numerical analysis, signal processing [6, 7], image processing, and so forth. Many years ago, wavelet methods were used for numerical analysis, particularly the numerical solution of PDEs. Up to now, researchers have utilized the simplest Haar wavelets to solve kinds of PDEs. Chen and Hsiao, in [8], established an operational matrix of integration based on Haar wavelets and used a procedure for applying the matrix to obtain wavelet solution of PDEs. In [9, 10], Cattani solved Poisson's problem and Fredholm type integral equations by Harmonic wavelet method. Other wavelets are also extensively used to solve the kinds of PDEs, in which Shannon wavelet is applied in the numerical solution of some equations, such as [11, 12]. Shannon scaling function and Sinc function combined with other methods (Galerkin, etc.) are used to solve some PDEs [13, 14]. In light of the above description, we are enlightened that Shannon wavelet is a useful tool to obtain the solution of convection-diffusion equations, which combined with Rugge-Kutta-Gill method.

In this paper, the content is assigned as follows. In Section 2, Shannon wavelet is introduced. We elaborate the concrete method solving convection-diffusion equation in Section 3. In Section 4 the viability of Shannon wavelet collocation is tested by a listed example.

2. Preliminaries

2.1. Shannon Wavelet

Wavelets are classified as families with names, such as Haar wavelet, Meyer wavelet, and Shannon wavelet. Shannon wavelets are the real part of harmonic wavelets. They have a slow decay in the time domain but a very sharp compact support in the frequency (Fourier) domain. This fact, together with the Parseval's identity, is used to compute the inner product and the expansion coefficients of the Shannon wavelets easily. A set of Shannon scaling functions in the subspace is defined as and the mother wavelets are In (1) and (2), the scaling function and mother wavelet for (see Figure 1) are as follows:

To some properties of Shannon scaling and wavelet functions, Cattani has detailedly researched in [1518]. So in this paper, we will not narrate the properties again.

To (3), its Fourier transform . It is very easy to see that

According to this equation, the sequence of function is orthonormal. A reproducing kernel is generated [19] as follows: Recomposing (6), we obtain a new reproducing kernel where is the spatial mesh size.

In one-dimensional function , we make the domain be discrete and set the grid size So we obtain collocation points where is a number of nodes, which used in the discretization and also is the maximum wavelet index number. Now, a basis function of Shannon wavelet will be constructed by (7) It has some properties as follows.(i)To the random , the function fulfills interpolation property: (ii)We have noticed that the constructed basis functions are orthogonal to each other as follows: (iii)If we make the integral with the basis functions and their derived functions, we obtain.

Both and its associated wavelet play an important part in signal processing. Unfortunately, when , the reduction of is very slow. So our paper only researches the case which belongs to finite interval.

2.2. Function Approximation

According to Shannon's sampling theorem, any function can be denoted as [19] where the coefficients is the value of the function at the point . is the Paley-Wiener reproducing kernel Hilbert space which is a subspace of the Hilbert space . In , can be approximated by . So we get

3. Method of Solution of Convection-Diffusion Equation

In this section, let us consider the one-dimensional convection-diffusion equation with constant coefficients: with initial condition and boundary conditions: Like (8) and (9), we will also divide the interval into equal parts of length and denote . We know that can be approximated by expanded in terms of the constructed basis function as formula (15) as follows: We multiply formula (16) with the constructed basis function , then we obtain Integrate that formula (19) with respect to from to as follows: The left expression of formula (20) is as follows: The right expression of formula (20) is as follows: Via (21) = (22), we obtain the following expression:

To any , we will get one equation. So equations will be obtained. In order to simplify the equations, we define the matrices and as follows: Combining with (24), the formula (23) is evolved into a matrix equation: Now we will use Runge-Kutta-Gill method to solve formula (25) as follows: where is the time interval. From (9) and (17), the initial value is obtained, and then we can evaluate the numerical solution at any collocation point within the different parameter .

4. Test of Example

A concrete convection-diffusion equation has known exact solution will be considered, and we observe how well the Shannon wavelet solution approximates the exact solution. We assume that , For which the exact solution is

In the course of the experiment, we got and setted . We got the approximate charts in the case of and and obtained the conclusion that the wavelet solution is approximate to the exact solution more precisely (see Figure 2).

5. Conclusion

In this paper, the theory of Shannon wavelet combined with Runge-Kutta-Gill method is used to solve the approximation of convection-diffusion equations. It has been shown that the key idea of shannon wavelet collocation method is to transform convection-diffusion equations into one-dimensional equations at collocations points and to solve the problem via Rugge-Kutta-Gill method.


The authors are thankful to the referees for their valuable comments and suggestions that improved the presentation of this paper. This work was supported by the National Natural Science Foundation of China (11271001).