Back
Texte intégral
Journal article

Optimized point shifts and poles in the linear rational pseudospectral method for boundary value problems

    30.11.2004
Published in:
  • Journal of Computational Physics. - 2005, vol. 204, no. 1, p. 292-301
two-point boundary value problems
linear rational collocation
pole optimization
mesh generation
point shift optimization
English Due to their rapid – often exponential – convergence as the number N of interpolation/collocation points is increased, polynomial pseudospectral methods are very efficient in solving smooth boundary value problems. However, when the solution displays boundary layers and/or interior fronts, this fast convergence will merely occur with very large N. To address this difficulty, we present a method which replaces the polynomial ansatz with a rational function r and considers the physical domain as the conformal map g of a computational domain. g shifts the interpolation points from their classical position in the computational domain to a problem-dependent position in the physical domain. Starting from a map by Bayliss and Turkel we have constructed a shift that can in principle accomodate an arbitrary number of fronts. Its parameters as well as the poles of r are optimized. Numerical results demonstrate how g best accomodates interior fronts while the poles also handle boundary layers.
Faculty
Faculté des sciences et de médecine
Department
Département de Mathématiques
Language
  • English
Classification
Mathematics
License
License undefined
Identifiers
Persistent URL
https://folia.unifr.ch/unifr/documents/299656
Statistics
Document views: 53 File downloads:
  • 1_berrut_ops.pdf: 96