Hi meep users

 

I want to get the energy flux in dielectric medium about different
wavelength

 

So I use ‘flux-in-box’ to get the poynting vector integral results. 

 

Other conditions are same.  I just change only the wavelength.

 

For example, if I put the wavelength 0.44, the integral results have ‘0.01
- 0.02’ range

 

But, if I put the wavelength 0.52, the integral results have ‘2 - 2.5’
range.

 

It is a very big difference

 

I can’t understand this results 

 

I hope your advises.

 

My ctl file is below. 

 

Thanks in advance

 

Cho.

 

 

(define sio (make dielectric (index 1.5)))

(define-param dpml 1) ; PML thickness; The cell dimensions

 

(set! geometry-lattice (make lattice (size 12 12 12)))

(set! geometry

(list (make block (center 0 1.5 0) (size infinity 3 infinity)

              (material (make dielectric (index 1.5))))

             (make block (center 0 3.1 0) (size infinity 0.2 infinity)

              (material (make dielectric (index 1.8))))

             (make block (center 0 3.29 0) (size infinity 0.18 infinity)

              (material (make dielectric (index 1.75))))

             (make block (center 0 3.43 0) (size infinity 0.1 infinity)

              (material (make dielectric (epsilon -1e20))))

             (make block (center 5.5 0.25 0)  (size 1 6.5 infinity)

              (material (make dielectric (epsilon -1e20))))

             (make block (center -5.5 0.25 0) (size 1 6.5 infinity)

              (material (make dielectric (epsilon -1e20))))

             (make block (center 0 0.25 5.5) (size infinity 6.5 1)

              (material (make dielectric (epsilon -1e20))))

             (make block (center 0 0.25 -5.5) (size infinity 6.5 1)

              (material (make dielectric (epsilon -1e20))))))

 

(set! pml-layers (list (make pml (thickness dpml))))

(set-param! resolution 10)

(set! Sources

(list

          (make source

                (src (make continuous-src (wavelength 0.44)))

                (component Ez) (center 0 3.3 0) (size 0 0 0))))

 

(define (flux) (print "flux: "

   (flux-in-box Y (volume (center 0 -3 0) (size 10 0 10))) "\n"))

 

(run-until 200

(at-every 0.4 flux))

 

 

 

 

 

 

_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss

Reply via email to