Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 974808, 7 pages

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

Research Article

A New Numerical Method of Particular Solutions for Inhomogeneous Burgers’ Equation

Huantian Xie,1,2 Dingfang Li,² and Feng Li¹

¹School of Science, Linyi University, Linyi 276005, China

²School of Mathematics and Statistics, Wuhan University, Wuhan 430072, China

Received 9 May 2013; Accepted 28 June 2013

Academic Editor: Fazal M. Mahomed

Copyright © 2013 Huantian Xie 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.

Based on the finite difference scheme in time, the method of particular solutions using radial basis functions is proposed to solve one-dimensional time-dependent inhomogeneous Burgers’ equations. Two numerical examples with good accuracy are given to validate the proposed method.

1. Introduction

In this paper we consider the one-dimensional nonlinear evolutionary partial differential equation

The corresponding homogeneous equation

was first introduced by Bateman [1] who considered its steady state solutions. Later, Burgers [2, 3] treated it as a mathematical model for free turbulence, and subsequently this equation is widely referred to as Burgers’ equation. Burgers’ equation can model several physical phenomena such as traffic, shock waves, and continuous stochastic processes. It can also be used to test various numerical algorithms. Due to its wide range of applicability, several researchers have been interested in the properties of its solution.

Burgers’ equation has been solved analytically for a restricted set of arbitrary initial and boundary conditions [4, 5]. Benton and Platzman [6] surveyed about 35 distinct exact solutions of the one-dimensional Burgers-like equations and their classifications. It is well known that the exact solution of Burgers’ equation can only be computed for restricted values of . Therefore, various numerical methods were employed to obtain the solution of Burgers’ equation. It is not our purpose to exhaust all these numerical schemes. Nevertheless, the solution methodologies commonly fall into the following classes: finite difference method (FDM), finite element method (FEM), and spectral methods. A survey of these techniques is given in [7, 8]. The previous numerical methods all depend on the mesh of the studied domain [9–12].

To alleviate the difficulty of mesh generation, various meshless techniques have been introduced during the past two decades. In a meshless (meshfree) method, a set of scattered nodes are selected in the computational domain. Meshless schemes include the method of fundamental solutions (MFS) [13], the method of particular solutions (MPS) [14–17], the element-free Galerkin method [18], local point interpolation [12], and boundary knot method [19]. It is known that the MFS is a boundary-type meshless method which is highly accurate for solving homogeneous equations if the fundamental solution of the given differential operator is known [20]. However, the fundamental solution of a given differential equation is not always available and often very difficult to derive. The ill-conditioning of the matrix resulting from the formulation of using the MFS and the location of source points are still outstanding research problems. To extend the MFS to inhomogeneous equations or time-dependent problems, the MPS has been introduced to evaluate the particular solution of the given differential equation. Since the particular solution is not unique, there is a rich variety of numerical techniques developed for this purpose.

Radial basis functions (RBFs), polynomial functions, trigonometric functions, and so forth [20–24], have been employed as the basis functions to approximate the particular solutions for the given differential equation. Once a particular solution has been evaluated, the given inhomogeneous equations can be reduced to the homogeneous equation. The original differential equation can be recovered by adding the homogeneous solution and the particular solution. This is a two-stage numerical scheme and is a well-known procedure for solving linear partial differential equations.

In general, the fundamental solution can be viewed as a special type of particular solution. When the inhomogeneous term is replaced by the delta function, a particular solution becomes a fundamental solution. The main idea of the MFS is that a fundamental solution satisfies the homogeneous equation inside the domain, and one only needs to enforce the fundamental solution on the boundary conditions to obtain the solution of the given homogeneous problem. Motivated by a similar idea, the particular solution can be used to solve inhomogeneous problems; that is, since the particular solution satisfies the given inhomogeneous equation through the domain without satisfying the boundary conditions, one only needs to impose the boundary conditions to obtain the solution of the given inhomogeneous problem.

It is the purpose of this paper to extend the MPS to a one-stage numerical scheme for solving one-dimensional time-dependent Burgers’ equations through the use of RBFs and then solve the time-dependent problems without the need of a two-stage numerical scheme [25] for obtaining homogeneous solution.

2. The Method of Particular Solutions

Consider the following boundary value problem:

where and are the Laplace operator and boundary differential operators, respectively, is the solution domain, is its boundary, and and are given functions.

Approximate by a finite series of RBFs through interpolation, and interpolants to can be constructed as

in which is Euclidean distance, is a set of interpolation points, and the real coefficients are to be determined by solving

if the real coefficient matrix is positive definite.

Therefore, from (6) an approximate particular solution to (3) is given by

where is obtained analytically by solving

If we impose in (7) to satisfy the governing equation in (3) and boundary conditions in (4), then becomes an approximate solution of the original partial differential equations (3)-(4). To be more specific, we have

For the numerical implementation, we let be the interior points, the boundary points, and . From (9) we have

The above system of equations can be easily solved using a standard matrix solver. Once the are determined, the approximate particular solution becomes the approximate solution of (3)-(4); that is,

Note that an accurate approximation of the particular solution depends on the appropriate choice of radial basis function . In the RBF literature [20, 26], some of the globally defined RBFs are only conditionally positive definite [26]. The unique solvability of the interpolation problem can be obtained by adding a polynomial term to the interpolation (5), giving

along with the constraints

where is a basis of , the space of -variate polynomials of order not exceeding , and

is the dimension of .

There are many types of globally defined RBFs [20], and the most popular RBFs are

 Inverse multiquadric (IMQ),

 Multiquadric (MQ),

 Gaussian (G),

 Polyharmonic (PH),

 Polyharmonic (PH),

3. Methodology

In order to solve such whole-space problems by numerical methods, we limit our consideration to a finite subdomain . In other words, function satisfies the following general nonlinear one-dimensional time-dependent Burgers’ equation:

with the initial condition

and the Dirichlet boundary condition

where is interpreted as the Reynolds number, and is the kinematic viscosity, and ,  ,  , and are known functions.

In the following section, a generalized trapezoidal method (-method) is used to approximate the time derivative in (20). Let be the time step and . For any and ,   can be approximated as follows:

Then,

where . For simplicity, we denote ,  , and . Substituting (23)–(24) into (20)–(22), we obtain the following equation:

Rewrite the previous equation as follows:

Then,

Assume that is a sought solution to the elliptic PDE. We can represent the right hand side of (27) as a function . This means that (27) is a standard Poisson-type differential equation

Therefore, if the fictitious function is known, (27) is equivalent to the Poisson-type equation (28) under the same boundary conditions.

Uniformly choose collocation points in the interior of domain and two boundary points and . For implementation, let ,  . Approximating the function by RBFs , we have

where . Then we can approximate at time step as follows:

where is obtained by solving

We use two RBFs, namely, IMQ and PH, where

It is easy to obtain the following :

Note that (27) is a recursion formula, and we can solve each elliptic PDE step by step starting with initial condition (21). As is well known, it is difficult to obtain an accurate numerical derivative from scattered data. Therefore, we choose in our method in order to avoid evaluating in (27). In this case, we can reformulate (27) as follows:

The nonlinear term is linearized as follows:

Substituting (35) into (34) and rearranging, we obtain

Write (30) together with boundary condition (22) in matrix form

where ,  , and ,  ,  . There are internal points and boundary points. The matrix can be split into , where

Applying this to the domain points and boundary points, (37) and (22) can be reformulated in the following matrix form:

where is the gradient differential operator, is a diagonal matrix with as its main diagonal, , and the accent “” means component by component multiplication of two vectors.

4. Numerical Results

Two different problems are used to test the accuracy of our method. In order to evaluate the numerical errors, we adopt three kinds of errors defined by

where is the exact analytical solution, and is the numerical solution of .

Example 1. Consider the following nonlinear one-dimensional time-dependent Burgers’ equation with a large Reynolds number , in the square domain ,

with the initial condition

and the boundary condition

The analytical solution is given as

We choose two RBFs, namely, IMQ and PH, as defined in (32). The ,  , and root-mean-square (RMS) errors for our numerical solutions are shown in Table 1 for ,  ,  , and 2. It can be seen that the accuracy of the method using IMQ is much higher than that using PH.

The graph of the analytical and estimated solutions for is shown in Figure 1. The absolute error graph is shown in Figure 2. We also show the space-time graph of the estimated solution in Figure 3.

Table 1: , , and RMS errors, with , , and [0, 2].

RMS

IMQ
PH
IMQ
PH
IMQ
PH

0.5

1.0

1.5

2.0

Figure 1: Analytical and numerical solution at , with IMQ as RBF, ,  ,  and , for Example 1.

Figure 2: Absolute error at , with IMQ as RBF, ,  ,  and , for Example 1.

Figure 3: Space-time graph of the solution up to , with IMQ as RBF, ,  ,  and , for Example 1.

Example 2. We consider the second-order nonlinear Burgers’ equation

where depends upon the exact solution of (45) as follows:

We take the required initial and boundary functions from the exact solution in the domain . Similar to Example 1, we also choose two RBFs, namely, IMQ and PH, as shown in (32). In Table 2, we compute the ,   errors and RMS errors at ,  ,  , and 2, with the Reynolds number . It indicates that the case of IMQ also has higher accuracy than the case of PH. In Figure 4, we draw the graph of the analytical and estimated solutions, and in Figure 5 the absolute error graph is shown. Figures 4 and 5 are shown at . The space-time graph of the numerical solution is shown in Figure 6.

Table 2: , , and RMS errors, with , , and 0, 1.

t

RMS

IMQ
PH
IMQ
PH
IMQ
PH

0.5

1.0

1.5

2.0

Figure 4: Analytical and numerical solution at , with PH as RBF, ,  ,  and , for Example 2.

Figure 5: Absolute error at , with PH as RBF, ,  ,  and , for Example 2.

Figure 6: Space-time graph of the solution up to , with PH as RBF, ,  ,  and , for Example 2.

Remark 3. Note that, for the two examples shown previously, the RBFs are unconditionally positive definite to guarantee the solvability of the resulting systems. However, some RBFs are conditionally positive definite. These types of interpolation problems can be obtained by adding a polynomial term to the interpolation (12), and it is also easy to verify the efficiency of the proposed schemes for these cases.

5. Conclusions

In this paper we proposed and implemented the method of particular solutions to solve the one-dimensional time-dependent Burgers’ equation. The effectiveness of the computational scheme is well demonstrated. It must be emphasized that the choice of radial basis functions is a flexible feature of these methods. The radial basis functions can be globally supported, infinitely differentiable and contain free parameters, namely, shape parameters, which affect both accuracy of the solutions and conditioning of the collocation matrix. The optimal shape parameters in Examples 1 and 2 using IMQ and PH, respectively, for all the calculations were found experimentally. The optimal choice of the shape parameters in RBFs is still an outstanding research problem [27–29]. A similar approach can be extended to solving 2D or 3D time-dependent partial differential equations. These research topics will be the focus of future investigation.

Acknowledgments

The work has been supported by the National Natural Science Foundation of China (Grants nos. 51190094, 61271337, 11201211, and 11201212) and AMEP of Linyi University. Special thanks are given to the reviewers for constructive advice and thoughtful comments.

References

H. Bateman, “Some recent researches on the motion of fluids,” Monthly Weather Review, vol. 43, no. 4, pp. 163–170, 1915.

J. M. Burgers, “Mathematical examples illustrating relations occurring in the theory of turbulent fluid motion,” vol. 17, no. 2, pp. 1–53, 1939.

J. M. Burgers, “A mathematical model illustrating the theory of turbulence,” Advances in Applied Mathematics, vol. 1, pp. 171–199, 1948.

J. D. Cole, “On a quasi-linear parabolic equation occurring in aerodynamics,” Quarterly of Applied Mathematics, vol. 9, pp. 225–236, 1951.

E. Hopf, “The partial differential equation ,” Communications on Pure and Applied Mathematics, vol. 3, pp. 201–230, 1950.

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

I. A. Hassanien, A. A. Salama, and H. A. Hosham, “Fourth-order finite difference method for solving Burgers' equation,” Applied Mathematics and Computation, vol. 170, no. 2, pp. 781–800, 2005.

S. Kutluay, A. Esen, and I. Dag, “Numerical solutions of the Burgers' equation by the least-squares quadratic B-spline finite element method,” Journal of Computational and Applied Mathematics, vol. 167, no. 1, pp. 21–33, 2004.

M. Dehghan and M. Tatari, “Determination of a control parameter in a one-dimensional parabolic equation using the method of radial basis functions,” Mathematical and Computer Modelling, vol. 44, no. 11-12, pp. 1160–1168, 2006.

M. Dehghan, “Finite difference procedures for solving a problem arising in modeling and design of certain optoelectronic devices,” Mathematics and Computers in Simulation, vol. 71, no. 1, pp. 16–30, 2006.

S. Chantasiriwan, “Cartesian grid methods using radial basis functions for solving Poisson, Helmholtz, and diffusion-convection equations,” Engineering Analysis with Boundary Elements, vol. 28, no. 12, pp. 1417–1425, 2004.

G. R. Liu and Y. T. Gu, “Boundary meshfree methods based on the boundary point interpolation methods,” Engineering Analysis with Boundary Elements, vol. 28, no. 5, pp. 475–487, 2004.

G. Fairweather and A. Karageorghis, “The method of fundamental solutions for elliptic boundary value problems,” Advances in Computational Mathematics, vol. 9, no. 1-2, pp. 69–95, 1998.

A. S. Muleshkov, M. A. Golberg, and C. S. Chen, “Particular solutions of Helmholtz-type operators using higher order polyharmonic splines,” Computational Mechanics, vol. 24, no. 5-6, pp. 411–419, 1999.

C. S. Chen, C. M. Fan, and P. H. Wen, “The method of approximate particular solutions for solving elliptic problems with variable coefficients,” International Journal of Computational Methods, vol. 8, no. 3, pp. 545–559, 2011.

T. Jiang, M. Li, and C. S. Chen, “The method of particular solutions for solving inverse problems of a nonhomogeneous convection-diffusion equation with variable coefficients,” Numerical Heat Transfer A, vol. 61, no. 5, pp. 338–352, 2012.

T. Jiang, X. Wang, and Z. Zhang, “Solutions for one-dimensional time-dependent Schrodinger equations,” in Information Computing and Applications, pp. 371–378, Springer, Berlin, Germany, 2012.

T. Belytschko, Y. Y. Lu, and L. Gu, “Element-free Galerkin methods,” International Journal for Numerical Methods in Engineering, vol. 37, no. 2, pp. 229–256, 1994.

W. Chen, “Symmetric boundary knot method,” Engineering Analysis with Boundary Elements, vol. 26, no. 6, pp. 489–494, 2002.

M. D. Buhmann, Radial Basis Functions, Cambridge University Press, Cambridge, UK, 2003.

K. E. Atkinson, “The numerical evaluation of particular solutions for Poisson's equation,” IMA Journal of Numerical Analysis, vol. 5, no. 3, pp. 319–338, 1985.

A. H.-D. Cheng, “Particular solutions of Laplacian, Helmholtz-type, and polyharmonic operators involving higher order radial basis functions,” Engineering Analysis with Boundary Elements, vol. 24, no. 7-8, pp. 531–538, 2000.

H. A. Cho, M. A. Golberg, A. S. Muleshkov, and X. Li, “Trefftz methods for time dependent partial differential equations,” Computers, Materials and Continua, vol. 1, no. 1, pp. 1–38, 2004.

S. Yu. Reutskiy, C. S. Chen, and H. Y. Tian, “A boundary meshless method using Chebyshev interpolation and trigonometric basis function for solving heat conduction problems,” International Journal for Numerical Methods in Engineering, vol. 74, no. 10, pp. 1621–1644, 2008.

X. Wang, T. Jiang, and Z. Jiang, “PDEs solution based on two-stage MFS-MPS and onestage MFS-MPS,” in Information Computing and Applications, pp. 356–363, Springer, Berlin, Germany, 2012.

C. A. Micchelli, “Interpolation of scattered data: distance matrices and conditionally positive definite functions,” Constructive Approximation, vol. 2, no. 1, pp. 11–12, 1986.

S. Rippa, “An algorithm for selecting a good value for the parameter c in radial basis function interpolation,” Advances in Computational Mathematics, vol. 11, no. 2-3, pp. 193–210, 1999.

J. Wertz, E. J. Kansa, and L. Ling, “The role of the multiquadric shape parameters in solving elliptic partial differential equations,” Computers & Mathematics with Applications, vol. 51, no. 8, pp. 1335–1348, 2006.

H. Xie and D. Li, “A meshless method for Burgers' equation using MQ-RBF and highorder temporal approximation,” Applied Mathematical Modelling, 2013.