2D magnetic dipole - discrepancy with Green function results (FDTD)
I'm intending to simulate the magnetic field distribution emanating from a 2D magnetic dipole (line source) over a lossy substrate at wavelengths between 500 and 1000 nms. As a preliminary work, I'm trying to compare the the results sampled from frequency-domain monitor with those from analytically known Green's function expression. While the general shape of the line (Hankel function) fits perfectly, there is a wavelength-dependent strong discrepancy in terms of complex magnitudes between the obtained results; this does not converge or change no matter how I decrease the mesh size. I'm also using cwnorm and normalizing the (current-density dependent) Green function results since Lumerical uses unit dipole moments.
Am I doing something wrong with the calculations, or is there some sort of amplitude normalization that I'm missing out Thanks for the responses.
I've also attached the formula that I'm using; since Lumerical uses unit dipole moments I'm normalizing this one with (-1i*w).