Chebyshev solution of the nearly-singular one-dimensional Helmholtz equation and related singular perturbation equations: multiple scale series and the boundary layer rule-of-thumb

摘要

The one-dimensional Helmholtz equation, epsilon(2)u(xx) - u = f(x), arises in many applications, often as a component of three-dimensional fluids codes. Unfortunately, it is difficult to solve for epsilon 1 because the homogeneous solutions are exp(+/-x/E), which have boundary layers of thickness O(1/epsilon). By analyzing the asymptotic Chebyshev coefficients of exponentials, we rederive the Orszag-Israeli rule [16] that N approximate to 3/root3 Chebyshev polynomials are needed to obtain an accuracy of 1% or better for the homogeneous solutions. (Interestingly, this is identical with the boundary layer rule-of-thumb in [5], which was derived for singular functions like tanh([x - 1]/epsilon.) Two strategies for small epsilon are described. The first is the method of multiple scales, which is very general, and applies to variable coefficient differential equations, too. The second, when f(x) is a polynomial, is to compute an exact particular integral of the Helmholtz equation as a polynomial of the same degree in the form of a Chebyshev series by solving triangular pentadiagonal systems. This can be combined with the analytic homogeneous solutions to synthesize the general solution. However, the multiple scales method is more efficient than the Chebyshev algorithm when epsilon is very, very tiny.