Why Newton's method is hard for travelling waves: Small denominators, KAM theory, Arnold's linear Fourier problem, non-uniqueness, constraints and erratic failure

Authors:
John P. Boyd
Affiliations:
Department of Atmospheric, Oceanic and Space Science, 2455 Hayward Avenue, Ann Arbor, MI 48109, United States
Venue:
Mathematics and Computers in Simulation
Year:
2007

Citing 7
Cited 6

Solitons from sine waves: analytical and numerical methods for non-integrable solitary and cnoidal waves

Physica D
Numerical recipes in C (2nd ed.): the art of scientific computing

Numerical recipes in C (2nd ed.): the art of scientific computing
Numerical methods for bifurcations of dynamical equilibria

Numerical methods for bifurcations of dynamical equilibria
Numerical Methods, Software, and Analysis

Numerical Methods, Software, and Analysis
On backtracking failure in newton-GMRES methods with a demonstration for the navier-stokes equations

Journal of Computational Physics
Fourier pseudospectral method with Kepler mapping for travelling waves with discontinuous slope: application to corner waves of the Ostrovsky-Hunter equation and equatorial Kelvin waves in the four-mode approximation

Applied Mathematics and Computation
Newton Methods for Nonlinear Problems: Affine Invariance and Adaptive Algorithms

Newton Methods for Nonlinear Problems: Affine Invariance and Adaptive Algorithms

A generalized Petviashvili iteration method for scalar and vector Hamiltonian equations with arbitrary form of nonlinearity

Journal of Computational Physics
Iteration methods for stability spectra of solitary waves

Journal of Computational Physics
Adaptive rational spectral methods for the linear stability analysis of nonlinear fourth-order problems

Journal of Computational Physics
Newton-conjugate-gradient methods for solitary wave computations

Journal of Computational Physics
Simulation of coherent structures in nonlinear Schrödinger-type equations

Journal of Computational Physics
Numerical and perturbative computations of solitary waves of the Benjamin-Ono equation with higher order nonlinearity using Christov rational basis functions

Journal of Computational Physics

Quantified Score

Hi-index	0.03

Visualization

Abstract

Nonlinear travelling waves and standing waves can computed by discretizing the appropriate partial differential equations and then solving the resulting system of nonlinear algebraic equations. Here, we show that the ''small denominator'' problem of Kolmogorov-Arnold-Moser (KAM) theory is equally awkward for numerical algorithms. Furthermore, Newton's iteration combined with continuation in a parameter often exhibits ''erratic failure'' even in the absence of bifurcation. Wave resonances can interlock a countable infinity of branches in an extremely complex topology, as will be illustrated through the fifth-degree Korteweg-deVries equation. Continuation can easily jump, unsuspected, from one branch to another. Constraints, sometimes finite and sometimes infinite in number, are usually needed to specify a unique solution. This confluence of numerical difficulties can be overcome only by combining the latest numerical algorithms with a strong understanding of travelling wave physics.