Reconfigurable intelligent surfaces (RISs) are planar, almost 2-D, structures that can intelligently manipulate electromagnetic waves by low-cost near passive reflecting elements. RISs are considered a potential enabling technology for the sixthgeneration (6G) wireless communication systems due to their capability to tune wireless signals, thus smartly controlling propagation environments. This work has a twofold objective: firstly, we present a systematic methodology for characterizing RISs in the time domain, which rigorously includes propagation delays and mutual coupling among RIS elements; secondly, we analyze the RISs through a convolution-based solver that is enriched by a specialized convolution scheme that enables the analysis of systems with complex linear, nonlinear as well terminations. The proposed approach is validated by extensive numerical evaluations against well-established methods based on both frequency-and time-domain analyses, with a varying number of RIS elements, interdistance, and frequency of operation. Index Terms-Partial elements equivalent circuit method, reconfigurable intelligent surfaces (RISs), time-domain (TD) analysis.