Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 163734, 5 pages

http://dx.doi.org/10.1155/2013/163734

Research Article

A Shannon-Runge-Kutta-Gill Method for Convection-Diffusion Equations

Xiaoming Duan,¹ Jinsong Leng,¹ Carlo Cattani,² and Caiyun Li¹

¹School of Mathematical Sciences, University of Electronic Science and Technology of China, Chengdu, Sichuan 611731, China

²Department of Mathematics, University of Salerno, Via Ponte Don Melillo, 84084 Fisciano, Italy

Received 20 January 2013; Accepted 5 February 2013

Academic Editor: Shengyong Chen

Copyright © 2013 Xiaoming Duan 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.

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:

Figure 1: Shannon scaling function and mother function .

To some properties of Shannon scaling and wavelet functions, Cattani has detailedly researched in [15–18]. 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).

Figure 2: Comparison of the analytical solution and Shannon wavelet solution.

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.

Acknowledgment

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

References

B. B. King and D. A. Krueger, “The 1-D convection diffusion equation: galerkin least squares approximations and feedback control,” in Proceedings of the 43rd IEEE Conference on Decision and Control, pp. 1502–1507, 2004.

M. Bergdorf, G.-H. Cottet, and P. Koumoutsakos, “Multilevel adaptive particle methods for convection-diffusion equations,” Multiscale Modeling and Simulation, vol. 4, no. 1, pp. 328–357, 2005.

M. M. Gupta, R. P. Manohar, and J. W. Stephenson, “A fourth order, cost effective and stable finite difference scheme for the convection-diffusion equation,” in Numerical Properties and Methodologies in Heat Transfer, pp. 201–209, Hemisphere, 1983.

M. M. Gupta, R. P. Manohar, and J. W. Stephenson, “A single cell high order scheme for the convection-diffusion equation with variable coefficients,” International Journal for Numerical Methods in Fluids, vol. 4, no. 7, pp. 641–651, 1984.

Z. Shi, L. Y. Deng, and Q. J. Chen, “Numerical solution of differential equations by using Haar wavelets,” in Proceedings of the International Conference on Wavelet Analysis and Pattern Recognition, pp. 2–4, Beijing, China, 2007.

C. Cattani, “Wavelet based approach to fractals and fractal signal denoising,” Transactions on Computational Sciences Journal, vol. 5730, pp. 143–162, 2009.

J. S. Leng, T. Z. Hung, Y. F. Jing, and W. Jiang, “A study on conjugate quadrature filters,” EURASIP Journal on Advances in Singnal Processing, vol. 2011, Article ID 231754, 7 pages, 2011.

C. F. Chen and C. H. Hsiao, “Haar wavelet method for solving lumped and distributedparameter systems,” IEE Proceedings—Control Theory and Applications, vol. 144, pp. 87–94, 1997.

C. Cattani, “Harmonic wavelet solution of Poisson's problem,” Balkan Journal of Geometry and Its Applications, vol. 13, no. 1, pp. 27–37, 2008.

C. Cattani and A. Kudreyko, “Harmonic wavelet method towards solution of the Fredholm type integral equations of the second kind,” Applied Mathematics and Computation, vol. 215, no. 12, pp. 4164–4171, 2010.

C. Cattani, “Shannon wavelets for the solution of integrodifferential equations,” Mathematical Problems in Engineering, vol. 2010, Article ID 408418, 22 pages, 2010.

C. Cattani, “Second order Shannon wavelet approximation of -functions,” Politehnica University of Bucharest Scientific Bulletin A, vol. 73, no. 3, pp. 73–84, 2011.

V. G. Koures, “Solving the Coulomb Schrödinger equation in via sinc collocation,” Journal of Computational Physics, vol. 128, no. 1, pp. 1–5, 1996.

F. Stenger, “A “sinc-Galerkin” method of solution of boundary value problems,” Mathematics of Computation, vol. 33, no. 145, pp. 85–109, 1979.

C. Cattani, “Connection coefficients of Shannon wavelets,” Mathematical Modelling and Analysis, vol. 11, no. 2, pp. 117–132, 2006.

C. Cattani, “Shannon wavelets theory,” Mathematical Problems in Engineering, vol. 2008, Article ID 164808, 24 pages, 2008.

C. Cattani, “Signorini cylindrical waves and Shannon wavelets,” Advances in Numerical Analysis, vol. 2012, Article ID 731591, 24 pages, 2012.

C. Cattani, “Fractional calculus and Shannon wavelets,” Mathematical Problems in Engineering, vol. 2012, Article ID 502812, 26 pages, 2012.

G. W. Wei, “Quasi wavelets and quasi interpolating wavelets,” Chemical Physics Letters, vol. 296, pp. 215–222, 1998.