The k-Wave MATLAB toolbox is widely used to conduct medical ultrasound simulations. It uses a Fourier collocation method to numerically solve the governing model equations, and introduces sources by adding acoustic pressure at points on an orthogonal grid. This approach introduces two errors when sources don't exactly align with the grid. These are phase errors arising from shifting source points to nearby grid nodes, and amplitude errors arising from an angular dependence in the density of source points. These two errors are collectively referred to as `staircasing'. Staircasing errors can be overcome by considering the band-limited representation of sources that arises from the use of a Fourier collocation method. To do so, sources are discretised by convolving a band-limited point source with the desired source geometry. To validate this approach, a comparison is made with current k-Wave source algorithms and with the FOCUS ultrasound simulation code. The new sources are shown to eliminate staircasing errors
