Dear Steven and Meep's users,
I'm trying to simulate zero refractive index material with cw source point inside it. It was quite good results with two dimension calculations, there was no diverge at the point of the source. The three dimensional calculation of material block with source at the center don't converge and it is happen only inside the material. I also tried to use some sphere resonator around source point cause it's helps to avoid divergence at the center, but in this case I still can't obtain stationary results after large calculation time. I also can't understand why calculation is go down when I use rotate4-sym. Do somebody have any clue for my problems? File metamaterial.ctl: (set! geometry-lattice (make lattice (size 20 20 20))) (define-param f 0.25) (define-param eps_inf 1.001) (define-param mu_inf 1.001) (define-param w 1e-02) (define-param g 1e-15) (define sig (* (/ eps_inf (* w w) ) (- (* f f) (* w w) ))) (set! geometry (list (make block (center 0 0 0) (size 10 10 10) (material (make dielectric (epsilon eps_inf) (mu mu_inf) (E-polarizations (make polarizability (omega w) (gamma g) (sigma sig))) (H-polarizations (make polarizability (omega w) (gamma g) (sigma sig))) ))) (make sphere (radius 1) (center 0 0 0) (material (make dielectric (epsilon 1)))) )) (set! eps-averaging? true) (set! symmetries (list (make rotate2-sym (direction Z) ))) (set! sources (list (make source (src (make continuous-src (frequency f) (width 20))) (component Ez) (center 0 0 0) ))) (set! pml-layers (list (make pml (thickness 3)))) (set! resolution 10) (define T (/ 1 (* 4 f))) (set! force-complex-fields? false) (use-output-directory "meta-out") (run-until 1000 (after-time 10 ( in-volume (volume (center 0 0 0) (size 10 10 0)) (at-every T output-efield-z) )))
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
