Hello everyone,
I want to simulate a model in Geophysics . I only care about the conventional
conductivity σ. According to the "NOTE" of
"Conductivity and complex ε" in the "Materials in Meep" . I compute the σD ,
but I didn't get the right answer. And I get confused with σ and σD . And when
I do not set the value of ε (it default to be 1). The result is weird .
Here is my code :
*******************************************************************************************
(reset-meep)
(define-param res 8) ;define the reslution
(define-param dpml 2) ;the thickness fo PML
(set! geometry-lattice (make lattice (size 50 50 no-size)))
(set! geometry (list
(make block (center 0 0) (size 25 25 infinity)
(material (make medium (D-conductivity 7.33) (epsilon 5) )))))
(set! pml-layers (list (make pml (thickness dpml))))
(set! resolution res)
(set! sources (list
(make source
(src (make continuous-src (frequency 1)))
(component Ez)
(center 0 -12) (size 25 0))))
(run-until 100
(at-beginning output-epsilon)
(at-beginning output-mu)
(to-appended "ez" (at-every 0.4 output-efield-z)))
**********************************************************************************8
LErmeng
_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss