Deleted residuals, the QR-factored Newton iteration, and other methods for formally overdetermined determinate discretizations of nonlinear eigenproblems for solitary, cnoidal, and shock waves
Solitary waves, cnoidal waves, and shock waves can be computed by solving a nonlinear eigenvalue problem, which in discretized form is a system of nonlinear algebraic equations. Unfortunately, many such systems a-re singular because the solution is not unique until one or more additional constraints are imposed. For example, if the waves are translationally invariant and u (X) is a solution, then so also is u (X + Phi) for arbitrary (1). To obtain a unique solution, one must impose an additional condition to reduce the one-parameter family of solutions by constraining (D. We describe five methods for coping with such singular systems: (i) reformulation of the problem, (ii) deleting residuals, (iii) Keller's bordered matrix scheme, (iv) QR-factored. over-determined Newton iteration, and (v) pseudoinverse-Newton iteration. We illustrate these ideas using the cnoidal waves of the Korteweg-deVries equation, the traveling shocks of the Korteweg-deVries-Burgers equation. and the weakly nonlocal solitary waves of the nonlinear equatorial beta-plane equations. Finite difference, Fourier and rational Chebyshev pseudospectral methods. and spectrally upgraded finite difference,, are applied. Reformulation and deleting residuals are the cheapest strategies, but the QR-factored Newton iteration is needed for the shock waves, which lack the symmetry of the other two wave species. (C) 2002 Elsevier Science (USA).
solitary wave nonlinear algebraic equations spectrally upgraded finite differences newton's iteration