A technique for accurate collocation residual calculations

A multiple-grid collocation method is presented that allows exact evaluatio n of residuals generated by truncated trial function expansion solutions to boundary-value problems with polynomial nonlinearities. The method is used to formulate a true, discrete analog to the Galerkin projection applicable to the same class of problems. The numerical techniques developed are used to study the convergence behavior of a nonlinear, reaction-diffusion probl em as a function of Thiele modulus (phi) and trial function truncation numb er (N). The convergence problems encountered at high phi values are found t o result from a second, physically meaningless solution to the modeling equ ations. This 'spurious' solution and the true solution are involved in a sa ddle-node bifurcation that limits the range of phi where solutions are foun d for most finite N; the solutions appear to asymptotically approach each o ther as phi, N --> infinity regardless of the discretization method. The sa ddle-stable manifold of the spurious solution also defines the boundary of the set of initial conditions that diverge during dynamic simulations prior to the saddle-node bifurcation; all initial conditions are found to diverg e after this bifurcation point. (C) 1998 Elsevier Science S.A. All rights r eserved.