4 research outputs found
Spectral Solution with a Subtraction Method to Improve Accuracy
This work addresses the solution to a Dirichlet boundary value problem for the Poisson equation in 1-D, d2u/dx2 = f using a numerical Fourier collocation approach. The order of accuracy of this approach can be increased by modifying f so the periodic extension of the right hand side is suffciently smooth. A proof for the order is given by Sköllermo. This work introduces a subtraction technique to modify the function\u27s right hand side to reduce the discontinuities or improve the smoothness of its periodic extension. This subtraction technique consists of cosine polynomials found by using boundary derivatives. We subtract cosine polynomials to match boundary values and derivatives of f. The derivatives need only be calculated numerically and approximately represent derivatives at the boundaries. Increasing the number of cosine polynomials in the subtraction technique increases the order of accuracy of the solution. The use of cosine polynomials matches well with the Fourier transform approach and is computationally efficient. Implementation of this technique results in a solution with variable accuracy depending on the number of collocation points and approximated boundary derivatives. Results show that the technique can be up to 14th order accurate
Recommended from our members
Modified Fourier expansions: theory, construction and applications
Modified Fourier expansions present an alternative to more standard algorithms for the approximation of nonperiodic functions in bounded domains. This thesis addresses the theory of such expansions, their effective construction and computation, and their application to the numerical solution of partial differential equations.
As the name indicates, modified Fourier expansions are closely related to classical Fourier series. The latter are naturally defined in the d-variate cube, and, in an analogous fashion, we primarily study modified Fourier expansions in this domain. However, whilst Fourier coefficients are commonly computed with the Fast Fourier Transform (FFT), we use modern numerical quadratures instead. In contrast to the FFT, such schemes are adaptive, leading to great potential savings in computational cost.
Standard algorithms for the approximation of nonperiodic functions in -variate cubes exhibit complexities that grow exponentially with dimension. The aforementioned quadratures permit the design of approximations based on modified Fourier expansions that do not possess this feature. Consequently, such schemes are increasingly effective in higher dimensions. When applied to the numerical solution of boundary value problems, such savings in computational cost impart benefits over more commonly used polynomial-based methods. Moreover, regardless of the dimensionality of the problem, modified Fourier methods lead to well-conditioned matrices and corresponding linear systems that can be solved cheaply with standard iterative techniques.
The theoretical component of this thesis furnishes modified Fourier expansions with a convergence analysis in arbitrary dimensions. In particular, we prove uniform convergence of modified Fourier expansions under rather general conditions. Furthermore, it is known that the notion of modified Fourier expansions can be effectively generalised, resulting in a family of approximation bases sharing many of the features of the modified Fourier case. The purpose of such a generalisation is to obtain both faster rates and higher degrees of convergence. Having detailed the approximation-theoretic properties of modified Fourier expansions, we extend this analysis to the general case and thereby verify this improvement.
A central drawback of these expansions is that their convergence rate is both fixed and typically slow. This makes the construction of effective convergence acceleration techniques imperative. In the final part of this thesis, we design and analyse a robust method, applicable in arbitrary numbers of dimensions, for accelerating convergence of modified Fourier expansions. When employed in the approximation of multivariate functions, this culminates in efficient, high-order approximants comprising relatively small numbers of terms