Complex coordinate methods for hydrodynamic instabilities and Sturm-Liouville eigenproblems with an interior singularity
dc.contributor.author | Boyd, John P. | en_US |
dc.date.accessioned | 2006-04-07T19:09:43Z | |
dc.date.available | 2006-04-07T19:09:43Z | |
dc.date.issued | 1985-02 | en_US |
dc.identifier.citation | Boyd, John P. (1985/02)."Complex coordinate methods for hydrodynamic instabilities and Sturm-Liouville eigenproblems with an interior singularity." Journal of Computational Physics 57(3): 454-471. <http://hdl.handle.net/2027.42/25770> | en_US |
dc.identifier.uri | http://www.sciencedirect.com/science/article/B6WHY-4DDR4YF-W9/2/75fa39208412eef5f2c42a9a1611c858 | en_US |
dc.identifier.uri | https://hdl.handle.net/2027.42/25770 | |
dc.description.abstract | Calculations of inviscid, linearized waves in fluids are very difficult when a mean wind or current U(y) is included because the differential equation is singular wherever U(y)=c, the phase speed. These "critical latitude," "critical level," or "critical point" singularities are particularly severe for Chebyshev methods since these global expansion algorithms are very sensitive to the analytic properties of the solution. A simple remedy is described: by making a change of coordinates y=f(x) where y is the original variable and x is the new coordinate with f(x) a complex function, one can solve the problem on an arc in the complex plane that makes a wide detour around the singularity. Specific guidelines for choosing f(x) for different problems are given in the text. Results are impressive: for an eigenvalue problem with a pole in the middle of the original real interval (a "Sturm-Liouville problem of the fourth kind"), just six basis functions suffice to calculate the real and imaginary parts of the lowest eigenvalue to within 1.4%. For strong instability, i.e., modes whose phase speeds have large imaginary parts, the complex mapping is unnecessary because the critical latitudes are complex and distant from the real axis. Even so, the mapping is useful for instability problems because it can be used to make calculations for very slowly growing modes to follow the changes in c right up to the "neutral curve" where the imaginary part of c=0. Although especially valuable for spectral algorithms, the same trick can be applied with finite difference methods also. The main disadvantage of the algorithm is that the eigenfunction must be calculated in a second, separate step, but this is usually a minor flaw in comparison to the complex mapping's virtues for coping with singular eigenvalue problems. | en_US |
dc.format.extent | 1164222 bytes | |
dc.format.extent | 3118 bytes | |
dc.format.mimetype | application/pdf | |
dc.format.mimetype | text/plain | |
dc.language.iso | en_US | |
dc.publisher | Elsevier | en_US |
dc.title | Complex coordinate methods for hydrodynamic instabilities and Sturm-Liouville eigenproblems with an interior singularity | en_US |
dc.type | Article | en_US |
dc.rights.robots | IndexNoFollow | en_US |
dc.subject.hlbsecondlevel | Physics | en_US |
dc.subject.hlbsecondlevel | Mathematics | en_US |
dc.subject.hlbtoplevel | Science | en_US |
dc.description.peerreviewed | Peer Reviewed | en_US |
dc.contributor.affiliationum | Department of Atmospheric and Oceanic Science, University of Michigan, 2455 Hayward Avenue, Ann Arbor, Michigan 48109, USA | en_US |
dc.description.bitstreamurl | http://deepblue.lib.umich.edu/bitstream/2027.42/25770/1/0000331.pdf | en_US |
dc.identifier.doi | http://dx.doi.org/10.1016/0021-9991(85)90190-1 | en_US |
dc.identifier.source | Journal of Computational Physics | en_US |
dc.owningcollname | Interdisciplinary and Peer-Reviewed |
Files in this item
Remediation of Harmful Language
The University of Michigan Library aims to describe library materials in a way that respects the people and communities who create, use, and are represented in our collections. Report harmful or offensive language in catalog records, finding aids, or elsewhere in our collections anonymously through our metadata feedback form. More information at Remediation of Harmful Language.
Accessibility
If you are unable to use this file in its current format, please select the Contact Us link and we can modify it to make it more accessible to you.