This paper proposes a numerical model that can predict pulse-echo and pitch-catch signals in angle beam ultrasonic testing (UT). This modeling is based on the theory of a time-shift invariant system and each factor of ultrasonic measurement is individually characterized. We here treat an entire modeling of the angle beam UT including ultrasonic wave generation from an angle beam transducer, transmission from the transducer to a specimen and vice versa, and scattering by flaws in the specimen. A multi-Gaussian beam model is introduced into mathematical expressions of wave propagations in both transmission and reception processes. Scattering analyses are performed by means of the boundary element method. In this paper, the performance of this modeling is shown by simulations of scattered waveforms from an isolated internal crack and a surface-breaking crack.