Journal
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS
Volume 56, Issue 3, Pages 233-255Publisher
JOHN WILEY & SONS LTD
DOI: 10.1002/fld.1510
Keywords
high-order scheme; finite difference; piecewise polynomials; stability; Runge phenomenon; pseudospectra
Ask authors/readers for more resources
It is well known that high-order finite-difference methods may become unstable due to the presence of boundaries and the imposition of boundary conditions. For uniform grids, Gustafsson, Kreiss, and Sundstrom theory and the summation-by-parts method provide sufficient conditions for stability. For non-uniform grids, clustering of nodes close to the boundaries improves the stability of the resulting finite-difference operator. Several heuristic explanations exist for the goodness of the clustering, and attempts have been made to link it to the Runge phenomenon present in polynomial interpolations of high degree. By following the philosophy behind the Chebyshev polynomials, a non-uniform grid for piecewise polynomial interpolations of degree q <= N is introduced in this paper, where N + 1 is the total number of grid nodes. It is shown that when q = N, this polynomial interpolation coincides with the Chebyshev interpolation, and the resulting finite-difference schemes are equivalent to Chebyshev collocation methods. Finally, test cases are run showing how stability and correct transient behaviours are achieved for any degree q < N through the use of the proposed non-uniform grids. Discussions are complemented by spectra and pseudospectra of the finite-difference operators. Copyright (C) 2007 John Wiley & Sons, Ltd.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available