Future Far-IR space telescopes, such as the SAFARI instrument of the proposed JAXA/ESA SPICA mission, will use horn antennas to couple to cavity bolometers to achieve high levels of sensitivity for Mid-IR astronomy. In the case of the SAFARI instrument the bolometric detectors susceptibility to currents coupling into the detector system and dissipating power within the bolometers is a particular concern of the class of detector technology considered.1 The simulation of such structures can prove challenging. At THz frequencies ray tracing no longer proves completely accurate for these partially coherent large electrical structures, which also present significant computational difficulties for the more generic EM approaches applied at longer microwave wavelengths. The Finite Difference Time Domain method and other similar commercially viable approaches result in excessive computational requirements, especially when a large number of modes propagate. Work being carried out at NUI-Maynooth is utilising a mode matching approach to the simulation of such devices. This approach is based on the already proven waveguide mode scattering code "Scatter"2 developed at NUI-Maynooth, which is a piece of mode matching code that operates by cascading a Smatrice while conserving power at each waveguide junction. This paper outlines various approaches to simulating such Antenna Horns and Cavities at THz frequencies, focusing primarily on the waveguide modal Scatter approach. Recently the code has been adapted to incorporate a rectangular waveguide basis mode set instead of the already established circular basis set.