Solving Integral Equations on Piecewise Smooth Boundaries Using the RCIP Method: A Tutorial

Figure 14

The exterior Dirichlet problem for Helmholtz equation with RCIP applied to (55). The program demo11.m is used with Γ as in (1), and 𝜃=𝜋/2. The boundary condition 𝑔(𝑟) of (51) is generated by a point source at (0.3,0.1). (a) The absolute error in 𝑈(𝑟) at 𝑟=(−0.1,0.2). (b) The number of GMRES iterations needed to meet an estimated relative residual of 𝜖_mach.

(a)

(b)