This paper introduces a generalised 3rd-order Spectral Representation Method for the simulation of multi-dimensional stochastic fields with asymmetric non-linearities. The simulated random fields satisfy a prescribed Power Spectrum and Bispectrum. The general d-dimensional simulation equations are presented and the method is applied to simulate 2D and 3D random fields. The differences between samples generated by the proposed methodology and the existing classical Spectral Representation Method are analysed. An important feature of this methodology is that the formula can be implemented efficiently with the Fast Fourier Transform, details of which are presented. Computational savings are shown to grow exponentially with dimensionality as a testament of the scalability of the simulation methodology.