You can generate a plane wave of oblique incidence by using a couple
different tricks in meep. All of these have been discussed before on
previous threads:
http://ab-initio.mit.edu/pipermail/meep-discuss/2007-February/000700.html
http://ab-initio.mit.edu/pipermail/meep-discuss/2006-July/000228.html
http://article.gmane.org/gmane.comp.science.electromagnetism.meep.general/706
http://article.gmane.org/gmane.comp.science.electromagnetism.meep.general/245
http://article.gmane.org/gmane.comp.science.electromagnetism.meep.general/229
I recall from an older thread (which I can't find), that there's a way
to make a quasi-planewave from two point sources, placed at the corner
of the computation domain.
You can't arrange a line source to be oblique to the Yee grid, but you
can make a line source which is parallel to the grid and creates an
oblique wave.
If you want to ignore the wave on the opposite side of your source line,
place the source line at a PML interface. The PML will absorb the fields
on the other side.
Good Luck,
Matt
On Mon, 19 Feb 2007, MUHAMMAD Amin wrote:
Hi,
I want to discuss two issues:
1) I want to tilt line source at an angle to produce plane wave, I have
tried axis but it seems that it is not supported with source.
2) Making a single line source and radiating it at an angle e.g, 45
degree. The line source also radiates at an angle -45 degree. How can I
avoid it? i.e. propagation in +ve direction and not in -ve direction
(set! eps-averaging? false)
(define-param sxy 20)
(set! geometry-lattice (make lattice (size sxy sxy no-size)))
(set! pml-layers (list (make pml (thickness 1.0))))
(define ((my-amp k x0) x)
(exp (* 0+1i (vector3-dot k (vector3+ x x0)))))
(define-param fcen .5)
(define-param df 0.02)
(define-param kdir (vector3 1 1))
(define k (vector3-scale (* 2 pi fcen)
(unit-vector3 kdir)))
(set! k-point (vector3 0 0 0))
(set! sources
(list
;left head
(make source
(src (make continuous-src (frequency fcen)(fwidth no-size)))
(component Ez)(center (* s -0.3) 0)
(size 8 0)
(amp-func (my-amp k (vector3 (* s -0.3) 0))))))
Thanks
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
