High order Nystrom method for acoustic scattering

While high frequency approximation methods are widely used to solve flaw scattering in ultrasonic nondestructive evaluation, full wave approaches based on integral equations have great potentials due to their high accuracy. In this work, boundary integral equations for acoustic wave scattering are solved using high order Nystrom method. Compared with boundary elements method, it features the coincidence of the samples for interpolation basis and quadrature, which makes the far-field interaction free from numerical integration. The singular integral is dealt with using the Duffy transformation, while efficient singularity subtraction techniques are employed to evaluate the near singular integrals. This approach has the ease to go high order so highly accurate results can be obtained with fewer unknowns and faster convergence, and it is also amenable to incorporate fast algorithms like the multi-level fast multi-pole algorithm. The convergence of the approach for different orders of elements and interpolati...