janklinux
Newbie level 6
3d dft
Hi!
I have a working 3D FDTD with PML and I have checked it for reflections at the boundaries, it works fine. Anyway, I have a problem applying DFT (discrete fourier transform) to this. Here's what I do:
Hertzian Dipole in Gaussian Envelope (no difference in DFT with either hard or soft source)
Free Space for checking DFT
When I do DFT only the source comes out with a nice gaussian shape in DFT. All other spacial coordinates that I do a DFT at are not gaussian shaped but have more or less deep 'dints' in the shape. I can supply a pic if needed. The problem now is that in a lossless material there should be nearly the same signal whereever I do the DFT, right? Why doens't that work?
I have read in Taflov's book about numerical dispersion and the differences to them for waves travelling along the axes and along 2D-3D diagonals. But even if I do the DFT along one axis with like 10 grid points distance, they all look defferent. I'm starting to get mad at this I am interested in spectra of abs(ETrans)**2/abs(ESource)**2.
I would be really glad for a helping idea.
Thanx, jank
Hi!
I have a working 3D FDTD with PML and I have checked it for reflections at the boundaries, it works fine. Anyway, I have a problem applying DFT (discrete fourier transform) to this. Here's what I do:
Hertzian Dipole in Gaussian Envelope (no difference in DFT with either hard or soft source)
Free Space for checking DFT
When I do DFT only the source comes out with a nice gaussian shape in DFT. All other spacial coordinates that I do a DFT at are not gaussian shaped but have more or less deep 'dints' in the shape. I can supply a pic if needed. The problem now is that in a lossless material there should be nearly the same signal whereever I do the DFT, right? Why doens't that work?
I have read in Taflov's book about numerical dispersion and the differences to them for waves travelling along the axes and along 2D-3D diagonals. But even if I do the DFT along one axis with like 10 grid points distance, they all look defferent. I'm starting to get mad at this I am interested in spectra of abs(ETrans)**2/abs(ESource)**2.
I would be really glad for a helping idea.
Thanx, jank