We simulate the measurements of an active bifocal terahertz imaging system to reproduce the ability of the system to detect the internal structure of foams having embedded defects. Angular spectrum theory and geometric optics tracing are used to calculate the incident and received electric fields of the system and the scattered light distribution of the measured object. The finite-element method is also used to calculate the scattering light distribution of the measured object for comparison with the geometric optics model. The simulations are consistent with the measurements at the central axis of the horizontal stripe defects.