Dear meep user,
I am beginner in meep.  I want to simulate the thin film.
I want to learn transmission and the reflection when light is put.
I made the model while seeing the tutorial.  
I am suitable of this model made for myself or not confident. 
It suitable in this as the model?
My best regards. 
 
Input light  ramda = 500nm = 1.6678204
The glass substrate  n = 1.52 = epsilon = 2.3104
Thin film            n = 1.36 = epsilon = 1.8496
image
 |----------------|<- refl(reflection)
 |                |
 |      Air       |
 |                |
 |  Light sourse  |
 |                |
 |----------------|
 |   Thin film    |
 |----------------|
 |                |
 |   The glass    |
 |   substrate    |
 |                |
 |----------------|<- trans(transmission)

-----------------------------------------------------------------------------
(define-param sx 16) ; size of cell in X direction                              
(define-param sy 32) ; size of cell in Y direction                              
(set! geometry-lattice (make lattice (size 16 32 no-size)))
(define-param w 1) ;
(define-param no-bend? true) ;    

(set! geometry
          (list
           (make block
             (center 0 -13)
             (size infinity 6 infinity)
             (material (make dielectric (epsilon 2.3104)))) ; n=1.52
           (make block
             (center 0 -9.975)
             (size infinity 0.05 infinity)                  ; thin film 50nm
             (material (make dielectric (epsilon 1.8496)))))) ;n=1.36

(use-output-directory "trre_27_010")

(define-param fcen 1.6678204) ; pulse center frequency  ramda=500nm             
                  
(define-param df 0.5)    ; pulse width (in frequency)                           
  
(set! sources (list
               (make source
                 (src (make gaussian-src (frequency fcen) (fwidth df)))
                 (component Ez)
                 (center 0 0)
                 (size 0 w))))

(define-param reso 10); resolution

(set! pml-layers (list (make pml (thickness 3.0))))
(set-param! resolution reso)

(define-param nfreq 1001) ; number of frequencies at which to compute flux      
       
(define trans ; transmitted flux                                                
      (add-flux fcen df nfreq
                    (make flux-region
                     (center 0 -14) (size 32 0))))
(define refl ; reflected flux                                                   
      (add-flux fcen df nfreq
                 (make flux-region 
                   (center 0 14) (size 32 0))))

(if (not no-bend?) (load-minus-flux "refl-flux" refl))
(run-sources+ 500 (at-beginning output-epsilon))
(if no-bend? (save-flux "refl-flux" refl))

(run-sources+ 
 (stop-when-fields-decayed 50 Ez
                               (vector3 0 0))
                           1e-3)

(display-fluxes trans refl)
 
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss

Reply via email to