TY - JOUR
T1 - Spectrally accurate solution of nonperiodic differential equations by the Fourier-Gegenbauer method
AU - Vozovoi, L.
AU - Weill, A.
AU - Israeli, M.
PY - 1997
Y1 - 1997
N2 - It is well known that the Fourier partial sum of an analytic nonperiodic function, supported on a finite interval, converges slowly inside the interval and exhibits O(1) spurious oscillations near the boundaries (the Gibbs phenomenon). An effective algorithm which allows one to completely overcome the Gibbs phenomenon was developed in [J. Comput. Appl. Math., 43 (1992), pp. 81-92], The basic concept of this approach consists of the re-expansion of the Fourier partial sums into the rapidly convergent Gegenbauer series. In this paper we extend the Fourier-Gegenbauer (F-G) method of [J. Comput. Appl. Math., 43 (1992), pp. 81-92] to the evaluation of the spatial derivatives of a piecewise analytic function. Also, we apply this method to the solution of nonperiodic boundary value problems. Although the derivatives of a discontinuous function are not in L2, the exponential convergence of the truncated Gegenbauer series can be proved, and the rate of convergence can be estimated. The solution of differential equations is accomplished in two steps. First, a particular solution with arbitrary boundary conditions is constructed using the F-G method. This particular solution is then corrected to satisfy the prescribed boundary conditions of the problem by adding a proper linear combination of homogeneous solutions. For boundary layer problems the intermediate (particular) solution has steep profiles near the boundaries. These steep regions introduce a large error into the final solution, which presumably has a smooth profile on the whole interval. A method which compensates for this loss of accuracy by using the appropriately constructed homogeneous solutions is proposed.
AB - It is well known that the Fourier partial sum of an analytic nonperiodic function, supported on a finite interval, converges slowly inside the interval and exhibits O(1) spurious oscillations near the boundaries (the Gibbs phenomenon). An effective algorithm which allows one to completely overcome the Gibbs phenomenon was developed in [J. Comput. Appl. Math., 43 (1992), pp. 81-92], The basic concept of this approach consists of the re-expansion of the Fourier partial sums into the rapidly convergent Gegenbauer series. In this paper we extend the Fourier-Gegenbauer (F-G) method of [J. Comput. Appl. Math., 43 (1992), pp. 81-92] to the evaluation of the spatial derivatives of a piecewise analytic function. Also, we apply this method to the solution of nonperiodic boundary value problems. Although the derivatives of a discontinuous function are not in L2, the exponential convergence of the truncated Gegenbauer series can be proved, and the rate of convergence can be estimated. The solution of differential equations is accomplished in two steps. First, a particular solution with arbitrary boundary conditions is constructed using the F-G method. This particular solution is then corrected to satisfy the prescribed boundary conditions of the problem by adding a proper linear combination of homogeneous solutions. For boundary layer problems the intermediate (particular) solution has steep profiles near the boundaries. These steep regions introduce a large error into the final solution, which presumably has a smooth profile on the whole interval. A method which compensates for this loss of accuracy by using the appropriately constructed homogeneous solutions is proposed.
KW - Exponential accuracy
KW - Fourier series
KW - Gegenbauer series
KW - Gibbs phenomenon
KW - Green's functions
KW - Helmholtz equation
UR - http://www.scopus.com/inward/record.url?scp=0000626401&partnerID=8YFLogxK
U2 - 10.1137/S0036142994278814
DO - 10.1137/S0036142994278814
M3 - ???researchoutput.researchoutputtypes.contributiontojournal.article???
AN - SCOPUS:0000626401
SN - 0036-1429
VL - 34
SP - 1451
EP - 1471
JO - SIAM Journal on Numerical Analysis
JF - SIAM Journal on Numerical Analysis
IS - 4
ER -