The aim of the present work was to develop a method for simulating breast lesions in digital mammographic images. Based on the visual appearance of real masses, three dimensional masses were created using a 3D random walk method where the choice of parameters (number of walks and number of steps) enables one to control the appearance of the simulated structure. This work is the first occasion that the random walk results have been combined with a model of digital mammographic imaging systems. This model takes into account appropriate physical image acquisition processes representing a particular digital X-ray mammography system. The X-ray spectrum, local glandularity above the insertion site and scatter were all taken account during the insertion procedure. A preliminary observer study was used to validate the realism of the masses. Seven expert readers each viewed 60 full field mammograms and rated the realism of the masses they contained. Half of the images contained real, histologically-confirmed masses, and half contained simulated lesions. The ROC analysis of the study (average AUC of 0.58±0.06) suggests that, on the average, there is evidence that the radiologists could distinguish, somewhat, between real and simulated masses.