Silicon photonics integrates both photonic and electronic components on the same silicon chip and promises ultra-dense, high-bandwidth interconnects via wavelength division multiplexing (WDM). However, when verifying such silicon photonic systems, the existing IC simulators face challenges due to the WDM signals containing multiple frequency tones at ~200-THz with ~50-GHz spacing. This paper presents a systematic approach to modeling the silicon photonic elements and devices as equivalent multi-port transmission lines using XMODEL primitives and simulating the WDM link models in an efficient, event-driven fashion in SystemVerilog. The 5Gb/s, 3-channel WDM link models with micro-ring, Mach-Zehnder, and electro-absorption modulators demonstrate the simulation speeds of 4.2, 8.3, and 8.3 symbols/second, respectively.