On Aug 1, 2008, at 6:05 PM, Zhu, Lei wrote:
First, if I want to get scattered fields, I did the simulation
twice. One has object, one without object.
Then I use the fields (Ez, Hx, Hy) that I got with object subtract
the fields without the object. These fields are the scattered
fields, right?
Yes.
Second, for my ntf code, I want to make the code match the data
gotten by meep.
I set the parameters like this:
Eps0=1
Mu0=1
Z (intrinsic impedance of free space) =1
K=1 (the wave vector of incident wave)
Omiga (angular frequency)=1
Wavelength=a=1 (also a=wavelength)
Wl=40 (resolution)
Are these parameters right?
Especially for the omiga, for my ntf code I used to set eps0, mu0…
as the real value (8.85e-12…) and omiga=pi/wl, so now how should I
set the omiga then it can match the data from meep?
Your parameters are inconsistent - if you have a wavelength = 1, then
your wavevector and omega should be 1/2pi in vacuum for c=1 units.
In Meep, when you specify frequency or wavevector it already divided
by 2pi, so a wavelength of 1 in vacuum corresponds to a wavevector
with magnitude 1 and a frequency 1. So the above is fine and
consistent for Meep, but not for your NTF code unless your NTF code
also divides by 2pi implicitly like Meep does.
Steven
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
