The TLM method, being a time domain method, cannot deal directly with frequency dispersive media. A way for modelling such media, for instance, dielectric dispersive ones, is proposed, starting from the causality relationship between field vectors D and E, which is discretized over the same space–time mesh as the three-dimensional TLM mesh (in our case, with HSCN node). This leads to adding, at each node on the network, three voltage sources (one for each field component), their value depending on the previous instants. The method is validated by computing the reflection (magnitude and phase) coefficient at normal incidence at an air–dielectric interface, for different dispersive behaviours of the dielectrics. Copyright © 2001 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]