Recent studies have indicated that terahertz time-domain spectroscopy (THz-TDS) can stably and efficiently acquire output spectra using an affordable and compact multimode laser diode (MMLD) with delayed optical feedback as the light source. This research focused on a numerical analysis of the optimal conditions for employing an MMLD with delayed optical feedback (chaotic oscillating laser diode) in THz-TDS, utilizing multimode rate equations. The findings revealed that the intermittent chaotic output generated by the MMLD, characterized by concurrent picosecond pulse oscillations lasting several tens of picoseconds, proves to be highly effective for THz-TDS. By appropriately setting the amounts of injection current and optical feedback, and the delay time of optical feedback, intermittent chaotic oscillation can be attained within a considerably broad parameter range. Moreover, both the MMLD output spectrum and the THz-TDS output spectrum exhibit a consistently stable shape at the microsecond scale, demonstrating the attractor properties inherent in an MMLD with delayed optical feedback.