Premium
Stable high‐order finite‐difference methods based on non‐uniform grid point distributions
Author(s) -
Hermanns Miguel,
Hernández Juan Antonio
Publication year - 2007
Publication title -
international journal for numerical methods in fluids
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.938
H-Index - 112
eISSN - 1097-0363
pISSN - 0271-2091
DOI - 10.1002/fld.1510
Subject(s) - mathematics , piecewise , interpolation (computer graphics) , finite difference , chebyshev filter , boundary (topology) , polynomial , degree (music) , chebyshev polynomials , grid , degree of a polynomial , stability (learning theory) , collocation (remote sensing) , vandermonde matrix , finite difference method , mathematical analysis , geometry , computer science , animation , eigenvalues and eigenvectors , physics , computer graphics (images) , quantum mechanics , machine learning , acoustics
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 Sundström 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 Chebyshevinterpolation, 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 © 2007 John Wiley & Sons, Ltd.