The goal of using numerical methods in this study is two-fold: to replicate a set of measured, individualized HRTFs by a computer simulation, and also to visualise the resultant sound field around the head. Two methods can be wed: the Boundary Element Method (BEM) and the Infinite-Finite Element Method (IFEM). This paper presents the results of a preliminary study carried out on a KEMAR dummy-head, the geometry of which was captured with a high accuracy 3-D laser scanner and digitiser. The scanned computer model was converted to a few valid BEM and IFEM meshes with different polygon resolutions, enabling us to optimise the simulation for different frequency ranges. The results show a good agreement between simulations and measurements of the sound pressure at the blocked ear-canal of the dummy-head. The principle of reciprocity provides an effect method to simulate HRTF database. The BEM was also used to investigate the total sound field around the head, providing a tool to visualise the sound field for different arrangements of virtual acoustic imaging systems.