Dear Meep developers, dear users,
I am struggling with unexpected field distributions in a geometry that
should actually support plane waves and thus a homogenous energy density
(tot-pwr).
I use (0 0 0) k-point, complex fields and double-wavelength thickness for
the PMLs in Y-direction, a flat cw source some distance below the top pml
of the 2D (sx * sy) cell and an empty computational cell:
;Source:
(set! sources
(list
(make source
(src (make continuous-src (frequency (/ 1 wvl))))
(component Ez)
(center 0 (- (* 2.89 Adpml) (/ sy 2)))
(size sx 0 infinity)
)
)
)
;Geometry
(set! geometry
(list
(make block
(center 0 0)
(size infinity sy)
(material (make dielectric (epsilon1)))
)
)
)
;Boundaries/PML
(set! k-point (vector3 0 0 0))
(set! force-complex-fields? true)
(set! pml-layers (list (make pml (thickness Adpml) (direction Y)))
Normally, I would expect even wave fronts and at most some fluctuations of
the energy density in the Y direction. However, the energy density drops
significantly in the spots where the source collides with the cell edges,
and also shows a lattice-like fluctuation that seems to be an interference
pattern caused by some x-component in the wave vectors.
Can somebody give me a hint on what might cause this? I would be very
grateful for any help.
Kind regards,
Sebastian Wolf
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
