This paper presents a fast algorithm for transient simulation of power grids in VLSI systems using waveform relaxation techniques. Novel partitioning methods and convergence accelerators are developed for fast convergence of waveform relaxation iterations when applied to power grid networks. Unlike the direct solvers, the new method is highly parallelizable and yields significant speedups. Numerical examples are presented to demonstrate the validity and efficiency of the proposed method.