Received Date: February 10, 2015; Accepted Date: April 29, 2015; Published Date: May 08, 2015
Citation: Su LD, Jiang ZW, Jiang TS (2015) Numerical Method for One-Dimensional Convection-diffusion Equation Using Radical Basis Functions. J Phys Math 6:136. doi:10.4172/2090-0902.1000136
Copyright: © 2015 Su LD, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Visit for more related articles at Journal of Physical Mathematics
In this paper, the meshless method is employed for the numerical solution of the one-dimensional (1D)
convection-diffusion equation based on radical basis functions (RBFs). Coupled with the time discretization and the collocation method, the proposed method is a truly meshless method which requires neither domain nor boundary discretization. The algorithm is very simple so it is very easy to implement. The results of numerical experiments are presented, and are compared with analytical solutions to confirm the good accuracy of the presented scheme
Meshless method; Radical basis function (RBF); Numerical solution; Convection-diffusion equation
Whenever we consider mass transport of a dissolved species (solute species) or a component in a gas mixture, concentration gradients will cause diffusion. If there is bulk fluid motion, convection will also contribute to the flux of chemical species. Therefore, we are often interested in solving for the combined effect of both convection and diffusion.
The convection-diffusion equation is a combination of the diffusion and convection (advection) equations, and describes physical phenomena where particles, energy, or other physical quantities are transferred inside a physical system due to two processes: diffusion and convection .
In the above equation, four terms represents transient, convection, diffusion and source term respectively. Where u(x, t) is the variable of interest (species concentration for mass transfer, temperature for heat transfer), D is the diffusivity (also called diffusion coefficient), ν is the average velocity that the quantity is moving, R(x, t) is source term represents capacity of internal sources, ∇ represents gradient and ∇. represents divergence.
This paper is devoted to the numerical computation of the onedimension
(1D) convection-diffusion equation:
With the initial conditions:
With the initial conditions:
And Dirichlet boundary conditions:
Where α, β and ε are known constant coefficients, h(x) and gi (t) (i =0,1) are known continuous functions.
Recently, much attention has been given to the development, analysis, and implementation of stable methods for the numerical solution of the convection-diffusion equations (see  and the reference therein). Jim Douglas, et al.  combine definite element and finite difference methods based on the method of characteristic for solving the convection-diffusion problems. Chen and Hon  consider the 2D and 3D Helmholtz and convection-diffusion equation using boundary knot method. The meshless local Petro-Galerkin method for convection-diffusion equation was considered in . A new finite difference method described by Ram P. Manohar and John W. Stephenson .
In this article, we present a numerical scheme to solve the convection-diffusion equation using the collocation method with Radial Basis Function (RBF). The results of numerical experiments are presented, and are compared with analytical solutions to confirm the good accuracy of the presented scheme.
In last 25 years, the radial basis functions (RBFs) method is known as a powerful tool for scattered data interpolation problem. The use of RBFs as a meshless procedure for numerical solution of partial differential equations is based on the collocation scheme. Because of the collection technique, this method does not need to evaluate any integral. The main advantage of numerical procedures which use RBFs over traditional techniques is mesh-less property of these methods. RBFs are used actively for solving partial differential equations. The examples see [9-11]. In the last decade, the development of the RBFs as a truly meshless method for approximating the solutions of PDEs has drawn the attention of many researchers in science and engineering [12-14]. Meshless method has become an important numerical computation method, and there are many academic monographs are published [15-17].
The layout of the article is as follows: In section 2, we introduce the collocation method and apply this method on the convection-diffusion equation. The results of numerical experiments are presented in section 3. Section 4 is dedicated to a brief conclusion. Finally, some references are introduced at the end.
Radial basis function approximation
The approximation of a distribution u(x), using RBF may be written as a linear combination of N radial functions, usually it takes the following form:
, for (2.1)
Where N is the number of data points, x=(x1, x2…xd), d is the dimension of the problem, the λ’s are coefficients to be determined and φ is the radial basis function. Eq. (2.1.1) can be written without the polynomial ψ. In that case, φ must be unconditional positive definite to guarantee the solvability of the resulting system (e. g. Gaussian or Inverse Multi quadrics). However, ψ is usually required when φ is conditionally positive definite, i. e, when φ has a polynomial growth towards infinity. We will use the Multi quadrics (MQ), which defined as:
Where is the Euclidean norm. Since φ given by (2.2) is continuous, we can use it directly.
If denotes the space of d-variate polynomial of order not exceeding than q, and letting the polynomials (P1,P2, …, Pm) be the basis of in Rd, then the polynomial ψ(x) in Eq. (2.1) is usually written in the following form:
Where m = (q-1+d)!/(d!(q-1)!). To get the coefficients x=(λ1, λ2… λN), and the collocation method is used. However, in addition to the N equations resulting from collecting Eq. (2.1.1) at N points and extra m equations are required. This is ensured by the m conditions for Eq. (2.1),
In a similar representation as Eq. (2.1), for any linear partial differential operator u can be approximated by:
The convection-diffusion equation
Let us consider the 1D convection-diffusion equation Eq. (1.1), with the initial conditions Eq. (1.2) and the Dirichlet boundary conditions Eq. (1.3).
First, let us discretize Eq. (1.1) according to the following θ-weighted scheme:
where 0 ≤θ ≤1 , and τ is the time step size, and using the notation where , we get:
Assuming that there are N − 2 interpolation points, can be approximated by:
To guarantee the positive definition, here we use the following approximation:
Where is the Euclidean norm. The additional conditions due to Eq. (2.4) are written as:
Writing Eq. (2.9) together with Eq. (2.10) in a matrix form we have:
Where , and is given by:
Assuming that there are p < N − 2 internal points and N −2 − p boundary points, then the N ×N matrix A can be split into: , where
Using the notation to designate the matrix of the same dimension as A and containing the elements where , 1≤ i, j ≤ N , then Eq. (2.2.1) together with the boundary conditions Eq. (1.3) can be written in matrix form as:
Eq. (2.14) is obtained by combining Eq. (2.6), which applies to the domain points, while Eq. (1.3) applies to the boundary points. Together with the initial condition Eq. (1.2) and Eq. (2.14), we can get all λ’s, thus we can get the numerical solutions.
Since the coefficient matrix is unchanged in time steps, we use the LU factorization to the coefficient matrix only once and use this factorization in our algorithm.
Remark: Although Eq. (2.14) is valid for any value of , we will use (The famous Crank-Nicolson scheme).
In this section, we present several numerical results to confirm the efficiency of our algorithm for solving the 1D convection-diffusion equation.
In this example, we consider the convection-diffusion Eq. (1.1) in with , , , with the boundary conditions:
And the initial condition
Then the analytical solution of the equation is . The right side functions .
We use MQ radical basis function for the computation, the and RMS errors and Root-Mean-Square (RMS) of errors are obtained in Table 1 for t = 0.1, 0.3, 0.5, 0.7 and 1.0 with time steps τ = 0.001 and dx = 0.001 .
Table 1: For T=0.1, 0.3, 0.5, 0.7 and 1.0 with time steps T=0.001 and dx=0.001.
The space-time graph of analytical and numerical solutions for t=1 are given in Figure 1. Note that we cannot distinguish the exact solution from the estimated solution in Figure 1.
In this example, we consider Eq. (1.1) with α =1 , β = 0 , ε =1 and the boundary conditions:
And the analytical solution of the equation is given as:
We get the initial conditions from the exact solution. The right side functions
The and RMS errors and Root-Mean-Square (RMS) of errors are obtained in Table 2 for T=0.1, 0.25, 0.5, 0.75 and 1.0 with time steps τ = 0.001 and dx 0.01 .
Table 2: Numerical Errors at different times with T=0.001 and dx=0.001 for Example 2.
Similar to the previous example, the space-time graph of analytical and estimated solutions for t=1 are presented in Figure 2.
We consider the convection-diffusion equation Eq. (1.1) with α = −1 , β =10 and ε =1 in the interval , the exact solution is given as . The boundary conditions are:
The right side functions of , and we extract the initial conditions from the exact solution.
These results are obtained with dx = 0.001 , τ = 0.001 . Similar to the previous examples, the L∞ and L2 error and RMS errors for t =0.5, 0.75, 1.0, 1.25 and 1.5 are presented in Table 3.
Table 3: Numerical Errors at different times with T=0.001 and dx=0.001 for Example 3.
The space-time graph of analytical and numerical solution for t=1 is presented in Figure 3.
In this example, we consider the convection-diffusion equation in [0,1] has the following form:
The right side functions , with the boundary condition:
Then the analytical solution of the equation is we get the initial conditions from the exact solution.
In this case, we use the radial basis functions MQ for the discussed scheme. These results are obtained for dx = 0.001 , τ = 0.001. The graph of analytical and numerical solution for t=0.1, 0.2, 0.3, 0.4 and 0.5 is given in Figure 4. The results obtained show the very good accuracy and efficiency of the new approximate scheme. Note that we cannot distinguish the exact solution from the estimated solution in Figure 4.
We also give the difference between exact solutions and numerical solutions in Figure 5.
In this paper, the collocation method is employed for the numerical solution of convection-diffusion equation based on radical basis functions (RBFs). Coupled with the time discretization, the proposed method is a truly meshless method which requires neither domain nor boundary discretization. The results of numerical experiments are presented, and are compared with analytical solutions confirmed the good accuracy of the presented scheme.
This work was supported by the national Natural Science Foundation of China (Grant Nos.11301252, 11201212) and Applied Mathematics Enhancement Program of Linyi University.