Hi Joe
> two questions about Bloch bounday condition
> 1. Lets assume that the entire region is filled with a dielectric
material with eps 2. If I want to apply Bloch boundary condition should
I use f.use_bloch(K)) or f.use_bloch(K.sqrt(2)).
It depends on what you mean by K.
The bloch boundary condition is about the spatial dependence of the
fields, not the time dependence. So use wavevector, not frequency.
> In first one it is assumed that meep will take care of materials in
boundaries
What do you mean?
> 2. As far as I know, in Bloch theorem phase changes by K.r; so why
Steven used 2*pi*K.r in his sample code when he wanted to introduce a
new source:
Because Steven is smart and he defined the fourier transform the way
engineers and mathematicians do, i.e.
f(x)=int dk F(k) exp(2 pi i k x)
F(k)=int dx f(x) exp(2 pi i k x)
This connects more closely to the discrete fourier transform.
Physicists tend to write shorthand
f(x)=int dk F(k) exp(i k x)
but this is a bad idea, the 2 pi bites back:
F(k)=(1/(2pi)) int dx f(x) exp(i k x)
Good luck,
m.
joe ertaba wrote:
Dear Steven and Meep users
I have two questions about Bloch bounday condition
1. Lets assume that the entire region is filled with a dielectric
material with eps 2. If I want to apply Bloch boundary condition
should I use f.use_bloch(K)) or f.use_bloch(K.sqrt(2)). In first one
it is assumed that meep will take care of materials in boundaries
2. As far as I know, in Bloch theorem phase changes by K.r; so why
Steven used 2*pi*K.r in his sample code when he wanted to introduce a
new source:
http://www.mail-archive.com/[email protected]/msg00076.html
f.add_point_source(Ez, src, vec(xsrc2 = xsrc1, ysrc2 = ysrc1),
polar(1.0, 2*pi*(kx * xsrc1 + ky * ysrc1)));
I will appreciate responses.
Joe
------------------------------------------------------------------------
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
