Multiple precision pseudospectral computations of the radiation coefficient for weakly nonlocal solitary waves: Fifth-order Korteweg--DeVries equation
- Department of Atmospheric, Oceanic & Space Science, University of Michigan, Ann Arbor, Michigan 48109 (United States)
By performing multiple precision pseudospectral calculations using two different basis sets, we compute the radiation coefficient {alpha}({epsilon}) for very small {epsilon} to resolve discrepancies between earlier numerical work of the author`s and a small-{epsilon} perturbation theory of Grimshaw and Joshi. Multiple precision is needed because {alpha} is asymptotically proportional to exp({minus}{pi}/2{epsilon}), and is therefore below the single precision roundoff threshold when {epsilon}{lt}1/25. Because {alpha} decreases exponentially with 1/{epsilon}, we use a numerical method whose error decreases exponentially fast with the number of grid points {ital N}: a pseudospectral method or finite differences of as high as twenty-fourth order. Richardson iteration, preconditioning by high order finite differences, parity symmetry, Aitken and Richardson extrapolation, and the fast Fourier transform are all crucial in reducing our longest runs to about 24 h on a Unix workstation. Although preconditioning by second order differences is the norm, we find that high order preconditioning---as large as 14th order---is more efficient for our one-dimensional problem. We find that the discrepancies are mostly due to differences in convention for (i) definition of the parameter {epsilon} and (ii) choice of the far field phase. When these are accounted for, we obtain very good agreement with Grimshaw and Joshi`s first order term. However, their second order term, which is predicted to be {pi}{sup 2}/2{approx}4.94, is about {ital O}(0.0 to {minus}0.2) in our computations. The reason for this difference is still a mystery.
- OSTI ID:
- 46444
- Journal Information:
- Computers in Physics, Vol. 9, Issue 3; Other Information: PBD: May 1995
- Country of Publication:
- United States
- Language:
- English
Similar Records
Black hole evolution with the BSSN system by pseudospectral methods
Traps and snares in eigenvalue calculations with application to pseudospectral computations of ocean tides in a basin bounded by meridians