Hi, Steven and all.
I understood the source of the problem (thanks to this forum previous
posts). Nevertheless the problem still exists.
First of all a few words about what exactly I do: I make a computational
volume with periodic borders in X and Y direction, in Z dircetion I put pml.
However to make a plane wave under some angle, I make boundary conditions.
All of it looks like as:
(define angle 45)
(define kx 0)
(set! k-point (vector3 kx 0 0))
//-------------------------smth
else-------------------------------------------------
(set! pml-layers (list (make pml (thickness pml_width) (direction Z))))
;; establish periodic boundaries
(set! kx (* fcen (sin (deg->rad angle))))
(change-k-point! (vector3 kx 0 0))
Well after it I put continuous source:
(set! sources
(list
(make source
(src (make continuous-src
(frequency fcen)))
(amp-func (lambda (p) (exp (* 0+2i pi kx (vector3-x p)))))
(component Hy)
(center 0 0 (- (/ cell_numZ 2) pml_width) )
(size cell_numX cell_numY 0))))
So, I created a source but it looks like splited, i.e. it consists of parts
of different sign which lengths are equal to the half of the wavelength used
in the simulation. The point is that this is not a continuous source anymore
with uniform wavefront. How can it be transformed to normal state? One of
the possible variants is to make cell in X direction smaller but this does
not suits me.
Another question arose due to Matt's discussion last month. What is the
general recipe to calculate transmition/reflection spectrum in periodic
system like mine when the wave propagates under some angle (any angle). Am I
correct that: due to the phase boundary condition (where
kx=2pi/lambda, ->exp(i kx x) )
I need to make a simulation in the following way. I am going to compute
spectrum for one angle, thus I make multiple numerical experiments where I
make my narrowband gaussian source run in different ranges in order to cover
spectral range I want, of course I change my k-point kx in this case. Are
there any different ways to do this?
Many thanks in advance,
Anton.
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
