This work is motivated by the need for realistic ultrasonic probability of detection (POD) models in nondestructive evaluation (NDE). Past POD models have utilized flaw farfield scattering amplitudes along with other system parameters to predict the expected signal in postulated measurement geometries [1}. However, numerical evaluations of scattering amplitudes have generally been restricted to idealized fl aw shapes and , t o our knowledge , no scheme to calculate scattering amplitude s of arbitrary shape has ever been implemented in 3D. Volumetric shapes with an axis of symmetry have been examined with T-matrix and MOOT [2,3] but the axisymmetric limitation precludes a large portion of all expected flaw shapes. Furthermore, a quasi-plane wave assumption is often made. This assumption can become inappropriate for critical flaw sizes on the order of the beam size . A truly general POD model needs to have these assumptions removed.