We model the transport of particles present in a fluid steadily flowing through a porous medium. The porous medium is described by a representative three-dimensional network. The particles are subjected to advection by the flow and to thermal diffusion. We propose to calculate their trajectories with the continuous time random walk framework. This enables us to efficiently sample disordered networks with realistic topology. The method proposed in this paper is general and can be adapted to model dispersion of tracers. It is applied here to simulate the measurement of the flow propagator [Formula: see text] which is defined as the ensemble density distribution of tracer displacements [Formula: see text], in a given time interval delta t. It can be extracted from pulsed magnetic field gradient spin echo NMR experiments carried out on porous media while fluid is flowing. Preliminary numerical results show good qualitative agreement with experiments.