Hi,
I have been recently working on a metamaterial simulation control file it seems
to work for transmission but not for reflection. The reflection seems to
produce a reflectivity that is over unity. I am not sure what could be causing
this issue. Posted below is my control file. Let me know if you have a solution.
(reset-meep)
(define Cu (make medium (epsilon -1439.8) (D-conductivity (/ (* 2 pi 0.15
1066.7) 1))))
(define lattice-size 0.35) ; size of x - y square of computational space
(define-param lattice-length 10) ; length of computational space
(define-param pml-thickness 0.8) ; thickness of pml layer
(define-param substrate-index 3.42) ; refractive index of substrate
(define-param fcen 0.833) ; center wavelength of incident pulse
(define-param df 1.666) ; parameter related to the spectral width of the pulse
(define-param t 0.002) ; metal depth (z dir)
(define-param l 0.30) ; side-length
(define-param w 0.05) ; microsrtip width
(define TMO (make medium (epsilon 14.167) (D-conductivity 0.419316)))
(define-param g 0.05) ; gap-length
(set! resolution 100)
(set! eps-averaging? false)
(set! geometry-lattice (make lattice (size lattice-size lattice-size
lattice-length)))
(set! geometry
(list
(make block
(center 0 0 (* lattice-length .75)) (size lattice-size lattice-size
(/ lattice-length 2))
(material (make medium (index substrate-index)))
)
(make block
(center 0 0 0)
(size lattice-size lattice-size t)
(material
(make dielectric (epsilon 1))
)
)
)
)
(set! k-point (vector3 0 0 0)); enable periodic boundary condition
(set! pml-layers (list (make pml (thickness pml-thickness) (direction Z))));
set pml
(set! sources (list
(make source
(src(make gaussian-src (frequency fcen) (fwidth df)))
(component Ex)
(center 0 0 (- (- (/ lattice-length 2) pml-thickness .5)))
(size lattice-size lattice-size 0)
)
)
)
(run-until 160
(at-beginning output-epsilon)
(to-appended "ref35"
(in-volume (volume (center 0 0 (- (/ lattice-length 2) pml-thickness
.5)) (size 0 0 0)) output-efield-x))
)
(reset-meep)
(define Cu (make medium (epsilon -13536) (D-conductivity 402199)))
(define lattice-size 0.35) ; size of x - y square of computational space
(define-param lattice-length 10) ; length of computational space
(define-param pml-thickness 0.8) ; thickness of pml layer
(define-param substrate-index 3.42) ; refractive index of substrate
(define-param fcen 0.833) ; center wavelength of incident pulse
(define-param df 1.666) ; parameter related to the spectral width of the pulse
(define TMO (make medium (epsilon 14.167) (D-conductivity 0.419316)))
(define-param t 0.002) ; metal depth (z dir)
(define-param w 0.05) ; microsrtip width
(define-param d 0.3) ; dielectric thickness
(set! resolution 100)
(set! eps-averaging? false)
(set! geometry-lattice (make lattice (size lattice-size lattice-size
lattice-length)))
(set! geometry
(list
(make block
(center 0 0 (* lattice-length .75)) (size lattice-size lattice-size
(/ lattice-length 2))
(material (make medium (index substrate-index)))
)
(make block
(center 0 0 0)
(size lattice-size lattice-size t)
(material
(make dielectric (epsilon 1))
)
)
(make block
(center 0 0 0)
(size w l t)
(material Cu)
)
(make block
(center 0 0 0)
(size l w t)
(material Cu)
)
(make block
(center (/ l 4) (/ (- l w) 2) 0)
(size (/ l 2) w t)
(material Cu)
)
(make block
(center (/ (- l w) 2) (/ l -4) 0)
(size w (/ l 2) t)
(material Cu)
)
(make block
(center (/ l -4) (/ (- w l) 2) 0)
(size (/ l 2) w t)
(material Cu)
)
(make block
(center (/ (- w l) 2) (/ l 4) 0)
(size w (/ l 2) t)
(material Cu)
)
)
)
(set! k-point (vector3 0 0 0)); enable periodic boundary condition
(set! pml-layers (list (make pml (thickness pml-thickness) (direction Z))));
set pml
(set! sources (list
(make source
(src(make gaussian-src (frequency fcen) (fwidth df)))
(component Ex)
(center 0 0 (- (- (/ lattice-length 2) pml-thickness .5)))
(size lattice-size lattice-size 0)
)
)
)
(run-until 160
(to-appended "trangc1layer3035"
(in-volume (volume (center 0 0 (- (/ lattice-length 2) pml-thickness
.5)) (size 0 0 0)) output-efield-x))
(to-appended "reflgc1layer3035"
(in-volume (volume (center 0 0 (* -1 (- (/ lattice-length 2)
pml-thickness .5))) (size 0 0 0)) output-efield-x))
(at-beginning output-epsilon))
_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
