Hi, Meep community,
I need to calculate the frequency-domain green function at the location of
a current source. Following one of the previous discussion threads, I
introduced a "unit" current density in the form of delayed Gaussian. After
an enough long simulation time, I did a Fourier transform on the recorded
electric field at the point dipole before normalizing to the incident wave.
However, the calculated Green function is too oscillatory to have any
physical meaning. Could someone take a look at my input file below, and let
me know what went wrong?
Thanks.
Hanning Chen
Northwestern University
Input file:
(set! geometry-lattice (make lattice (size 80.0 80.0 80.0)))
(set! geometry (list
(make sphere (center 0.0 0.0 0.0) (radius 10.0)
(material (make dielectric (epsilon 3.189)
(E-polarizations
(make polarizability (omega 1.0E-6) (gamma
1.4435E-5) (sigma 5.4837E7) )
(make polarizability (omega 3.7643E-3)
(gamma 3.3385E-4) (sigma 0.432304) )
(make polarizability (omega 3.4192E-3)
(gamma 2.9998E-4) (sigma 0.223696) )
)
)
)
)
)
)
(set! sources (list
(make source
(src (make gaussian-src (frequency 1.6675E-3) (start-time
3000) (width 148.4924) ))
(component Ex)
(center 0 0 -15.0)
(size 0 0 0)
)))
(set! pml-layers (list (make pml (thickness 3.0))))
(set! resolution 4)
(run-until 30000
(at-beginning output-epsilon)
(at-every 30 (in-volume (volume (center 0 0 -15.0) (size 3.0 3.0
3.0) ) output-efield-x) )
(at-every 30 (in-volume (volume (center 0 0 -15.0) (size 3.0 3.0
3.0) ) output-efield-y) )
(at-every 30 (in-volume (volume (center 0 0 -15.0) (size 3.0 3.0
3.0) ) output-efield-z) )
)
_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
