The Trefftz method for the Helmholtz equation with degeneracy

  • Authors:

  • Z. C. Li

  • Affiliations:

  • Department of Applied Mathematics, Department of Computer Science and Engineering, National Sun Yat-sen University, Kaohsiung, Taiwan 80424 and National Center for Theoretical Sciences, Taiwan

  • Venue:
  • Applied Numerical Mathematics
  • Year:
  • 2008

Quantified Score

Hi-index 0.00

Visualization

Abstract

The Trefftz method (TM) [E. Trefftz, Ein Gegenstuck zum Ritz'schen Verfahren, in: Proc. 2nd Ind. Congr. Appl. Mech., Zurich, 1926, pp. 131-137] (i.e., Boundary approximation method) is developed to solve the Helmholtz equation, @Du+k^2u=0, where k^2 is not exactly equal (but may be very close) to an eigenvalue of the operator -@D. Piecewise particular solutions are chosen and then matched together in order to satisfy the exterior and interior boundary conditions. Error analysis is presented to estimate error bounds for the Helmholtz solutions in the entire solution domain. Let @d be the smallest relative distance between k^2 and the eigenvalues of -@D. We prove that the error asymptote of the solutions by the TM is O(1@d) as @d-0, which is called degeneracy in this paper. Such an asymptote O(1@d) has been verified by the numerical computations. We also explain why the exponential convergence rates of solutions can be obtained easily by splitting the solution domain into smaller subdomains. Numerical experiments of both smooth and singular solutions are provided, to support the TM algorithms and the error analysis made.