  Anisotropic and spatially dispersive graphene patches of arbitrary shape
embedded in planarly layered uniaxial media are analyzed using spectral-domain
method of moments. Formulation and computational methods for the
spectral-domain method of moments using the Rao-Wilton-Glisson subdomain basis
function and incorporating the full-wavevector Bhatnagar-Gross-Krook
formulation of graphene surface conductivity tensor are proposed. The impedance
matrix is efficiently evaluated by a novel numerical method which firstly
approximates the spectral-domain Green function, basis function and
conductivity tensor with Chebyshev polynomials, and then sums up the Fourier
transformed coefficients. Blue-shift of the resonant frequency and variation of
the current distribution due to spatial dispersion are observed in various
structures demonstrated.
