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

Figure 15

The exterior Neumann problem for Helmholtz equation with RCIP applied to (63). The program demo12.m is used with Γ as in (1), and 𝜃=𝜋/2. The boundary condition 𝑔(𝑟) of (59) 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)