Hi meep users,
I am trying to change the nonlinear index by control wave, but I can not.
Flux does not change even though kerr coefficient or control power is
increase.
I can see the change in flux at this frequency, if we increase the index
manually.
My ctl code is the following.
If someone shows me where I'm wrong I'll be very pleased.
Best regards…
Onder
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
(define-param npoly 1.46) ;; base index
(define-param epspoly (* npoly npoly))
(define-param nwg 1.488) ;; index of line defect waveguide
(define-param epswg (* nwg nwg))
(define-param epsair 1)
(define-param a 1)
(define-param r 0.2)
(define-param ax (* a 0.5))
(define-param ay (* a (/ (sqrt 3) 2)))
(define-param wgwidth (/ (* 1 a (sqrt 3)) 2)) ;; width of line defect
waveguide
(define-param sy 75)
(define-param sx 54)
(set! geometry-lattice (make lattice (size sx sy no-size)))
(set! default-material (make dielectric (epsilon epspoly)))
(define-param axx (* 2 a 0.5))
(define-param ayy (* (/ (sqrt 3) 2) 2 a))
(set! geometry
(append
(geometric-objects-lattice-duplicates
(list
(make cylinder (center 0 0 0) (radius r) (height infinity) (material
(make dielectric (epsilon epsair)))
)
(make cylinder (center ax (* -1 ay) 0) (radius r) (height
infinity) (material (make dielectric (epsilon epsair)))
)
)
axx ayy)
))
(set! geometry (append geometry
(list
(make block (center 0 0) (size sx 1 infinity)
(material (make dielectric (epsilon epspoly)
) ) ) ) ))
(set! geometry (append geometry
(list
(make block (center 0 0) (size sx wgwidth infinity)
(material (make dielectric (epsilon epswg)
) ) ) ) ) )
(define-param epsblok epswg) ; index of nonlinear block
(define-param kerr 0.59) ; kerr coefficient chi3
;; append nonlinear block in waveguide with Kerr
(set! geometry (append geometry
(list
(make block (center 0 0) (size 0.5 wgwidth infinity)
(material (make dielectric (epsilon epsblok) (chi3
kerr)
) ) ) ) ) )
(set! sources (list
(make source
(src (make continuous-src (frequency 0.78) )) ;;
control source
(component Hz)(center -22 0) (size 0 wgwidth) (amplitude 10)
)
(make source
(src (make continuous-src (frequency 0.60) (start-time 100) )) ;;
input source
(component Ez)(center -20 0) (size 0 wgwidth) (amplitude 1)
)
))
(set! pml-layers (list (make pml (thickness 2.0))))
(set! resolution 32)
(define-param flxdout 23) ; output flux region location
(define-param flxsize wgwidth)
(define transout
(add-flux 0.6 1 500
(make flux-region
(center flxdout 0) (size 0 flxsize)
) ) )
(run-until 250
(at-beginning output-epsilon)
(to-appended "ez" (at-every 5 output-efield-z)))
(display-fluxes transout)
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
