Rational Spectral Collocation Method for a Coupled System of Singularly Perturbed Boundary Value Problems

Suqin Chen, Yingwei Wang, Xionghua Wu
2011 Journal of Computational Mathematics  
A novel collocation method for a coupled system of singularly perturbed linear equations is presented. This method is based on rational spectral collocation method in barycentric form with sinh transform. By sinh transform, the original Chebyshev points are mapped into the transformed ones clustered near the singular points of the solution. The results from asymptotic analysis about the singularity solution are employed to determine the parameters in this sinh transform. Numerical experiments
more » ... e carried out to demonstrate the high accuracy and efficiency of our method. Mathematics subject classification: 65L10, 65M70. 459 method and Linβ [9] established a central difference scheme on certain layer-adapted meshes for the cases of m ≥ 2. Generally speaking, the solution to the problem of this type has two boundary layers with width O( √ ε) respectively at x = 0 and x = 1 under some assumptions. If B ≡ 0, the system is said to be of convection-diffusion type. Riordan et al presented a finite difference method consisting of upwinding on piecewise-uniform Shishkin meshes for the cases of m = 2 [10] and m > 2 [11] . They also used a Jacobi-type iteration to compute the solution [12] . Generally speaking, the problem of this type has a solution with a single boundary layer of width O(ε) at x = 0 (or x = 1) under proper assumptions. Furthermore, if (1.1) are coupled through their convective (first-order) terms (i.e. for each i = 1, · · · , m there exits a j = i such that b ij = 0), we say it is strongly coupled; otherwise, if B ≡ 0 or B is just a non-zero diagonal matrix, it is said to be weakly coupled. The spectral collocation method based on rational interpolants in barycentric form was proposed by Berrut and his collaborators [13] [14] [15] [16] . An advantage of which is that after transform, the derivatives in underlying differential equation are not required to be transformed correspondingly as is usual in other methods. Besides, Tee and Trefethen [17] devised a sinh transform that maps original Chebyshev points clustered near the boundaries of [−1, 1] into a new set of collocation points, say the transformed Chebyshev points, which are clustered near the singular point of a function. In order to determine the parameters in above sinh transform, singularities of the solution, including the location and width of the boundary layer, should be known. Hence, we resorted to the singularity location technique in asymptotic analysis and solved a parameterized singular perturbation problem [18] . Here we present a kind of numerical method based on rational spectral collocation in barycentric form with sinh transform (RSC-sinh method) for solving a coupled system of singularly perturbed problems in various types, both weakly coupled and strongly coupled. Numerical experiments illustrate that the RSC-sinh method enjoys improved spectral accuracy. This paper is organized as follows. The asymptotic analysis of the problem is outlined in Section 2. In Section 3, we construct the RSC-sinh method for problem (1.1)-(1.2). The numerical results of several examples are given in Section 4. Finally, we present some concluding remarks in Section 5.
doi:10.4208/jcm.1103-m3434 fatcat:c4i3rqhim5gjxigw6672q3gpii