As the science and engineering associated with terahertz time-domain spectroscopy and imaging evolves past the use of conventional free-space optics, the continued development of waveguides for terahertz pulses is increasingly relevant. The ability to model and simulate terahertz wave propagation aids in the development, visualization, and understanding of novel terahertz devices and phenomena. We discuss the use of the finite-element method, a powerful computational tool for the modeling of guided wave phenomena and devices at terahertz frequencies.