I did the simulation by using Meep 1.1.1 with libctl 3.1 and Guile 1.6.8, and h5topng 1.10.
$ meep metal-black.ctl
$ h5ls metal-black.h5
ez Dataset {400, 400, 83/Inf}
$ h5topng metal-black.h5 -Zc dkbluered -o metal-black.png -t 82 -A metal-black
-eps-000000.00.h5
attached: metal-black.png
Best
----- Original Message -----
>>I tried to modify the bend waveguide example from the website.
>The modifications, which I made were change the material from dielectric to m
etal and chopping the metal waveguide with air waveguide to make it to be a ho
llow metal waveguide.
>But I got problem, I got two black square and grow as the simulation time run
.... how that can happen and how to fix those black things??? (image attached)
>
>-------------------snip-------------------?
>(reset-meep)
>
>(define-param sx 40)
>(define-param sy 40)
>(set! geometry-lattice (make lattice (size? sx sy no-size)))
>
>(define-param pad 5) ; padding distance between waveguide and cell
> edge
>(define-param w 6) ; width of waveguide
>
>(define-param lam 5) ; wavelength
>
>(define-param thk 1) ; thickness of the
> hollow metal waveguide
>
>(define-param ctime (* 10 lam) ) ; computational time 10 periods
>
>(define wvg-ycen (* 0.5 (- sy w (* 2 pad)))) ; y center of horiz. wvg
>(define wvg-xcen (* 0.5 (- sx w (* 2 pad)))) ; x center of vert. wvg
>
>(set! geometry
>??? (
>??? ??? list
>??? ??? ??? (make block (center (* -0.5 pad) wvg-ycen) (size (- sx p
ad) w infinity)
>??? ??? ??? ??? (material metal)
>??? ??? ??? )
>??? ??? ??? (make block (center? wvg-xcen (* -0.5 pad)) (size w (-
sy pad) infinity)
>??? ??? ??? ??? (material metal)
>??? ??? ??? )
>??? ???
> ??? (make block (center (- (* -0.5 pad) (* 0.5 thk)) wvg-ycen) (size (- s
x (+ pad thk)) (- w (* 2 thk)) infinity)
>??? ??? ??? ??? (material air)
>??? ??? ??? )
>??? ??? ??? (make block (center? wvg-xcen (- (* -0.5 pad) (* 0.5 th
k))) (size (- w (* 2 thk)) (- sx (+ pad thk)) infinity)
>??? ??? ??? ??? (material air)
>??? ??? ??? )
>??? ??? ???
>??? )
>)
>
>(set! pml-layers
>??? (
>??? ??? list
>??? ??? ??? (make pml (thickness 1.0))
>??? )
>)
>
>(set! resolution 10)
>
>(set! sources
> ??? (???
>??? ??? ??? list
>??? ??? ??? (make source
>??? ??? ??? ??? (src (make continuous-src (wavelength lam)))
>??? ??? ??? ??? (component Ez)
>??? ??? ??? ??? (center (+ 1 (* -0.5 sx)) wvg-ycen)
>??? ??? ??? ??? (size 0 w)
>??? ??? ??? )
>??? ??? ???
>??? ??? )
>)
>
>(run-until ctime
>??? (at-beginning output-epsilon)
>??? (to-appended "ez" (at-every 0.6
> output-efield-z))
>)
>
>---------------------------snap---------------------------------------
>
>Thanks
>
>
>Sendy Phang
>Engineering Physics
>Computational Electromagnetics Resarch Group ?
>Bandung Institute of Technology
>Indonesia
<<attachment: metal-black-ez.png>>
_______________________________________________ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
