Numerical modeling of a nondestructive testing system plays an important role in many aspects of quantitative nondestructive evaluation (QNDE). The ultimate goal of a model is to predict test results for a specific flaw in a material. Thus, in ultrasonic testing, a system model should include the transducer, its radiation pattern, the beam reflection and propagation, and scattering from defects. In this paper attention is focused on the scattering model and the scattered fields by defects are observed by an elastodynamic boundary element method. Flaw types addressed are void-like and crack-like flaws. When transverse ultrasonic waves are obliquely incident on the flaw, the angular distribution of far-field scattered displacements are calculated and presented in the form of A-scan mode. The component signals obtained from each scattering problem are identified and their differences are addressed. The numerical results are also compared with those obtained by high frequency approximate solutions.