(define-param d (/ (sqrt 3) 2))
(define-param N 17)
(define-param sy (* N d)) ; the (odd) number of lateral supercell periods
(define-param dpml 1)
(define-param pad 1)
(define-param sx (+ (* (+ pad dpml) 2) N))
(define-param eps 12) ; the dielectric constant of the rods
(define-param r 0.45) ; the rod radius in the bulk crystal
(define-param rc 0.27)
(define-param s 1)
(define-param chirp 2)
(set! geometry-lattice (make lattice (size sx sy no-size)))
(set! geometry
(append
(list (make block (center 0 (* 0.5 d))
(size sx (* 2 d) infinity)
(material (make dielectric (epsilon eps)))))
(geometric-object-duplicates (vector3 1 0) 0 (/ (- N 1) 2)
(make cylinder (center 0 0) (radius r) (height infinity)
(material air)))
(geometric-object-duplicates (vector3 -1 0) 0 (/ (- N 1) 2)
(make cylinder (center 0 0) (radius r) (height infinity)
(material air)))
(geometric-object-duplicates (vector3 1 0) 0 (/ (- N 1) 2)
(make cylinder (center 0.5 (* 1 d)) (radius r) (height infinity)
(material air)))
(geometric-object-duplicates (vector3 -1 0) 0 (/ (- N 1) 2)
(make cylinder (center -0.5 (* 1 d)) (radius r) (height infinity)
(material air)))))
(set! geometry (geometric-objects-lattice-duplicates geometry infinity (* 2 d) infinity))
(set! geometry (append geometry
(geometric-object-duplicates (vector3 1 0) 0 (/ (- N 1) 2)
(make cylinder (center 0 0) (radius r) (height infinity)
(material (make dielectric (epsilon eps)))))
(geometric-object-duplicates (vector3 -1 0) 0 (/ (- N 1) 2)
(make cylinder (center 0 0) (radius r) (height infinity)
(material (make dielectric (epsilon eps)))))
(geometric-object-duplicates (vector3 1 0) 0 (- s 1)
(make cylinder (center 1 0) (radius rc) (height infinity)
(material air)))
(geometric-object-duplicates (vector3 -1 0) 0 (- s 1)
(make cylinder (center -1 0) (radius rc) (height infinity)
(material air)))
(geometric-object-duplicates (vector3 2 0) 0 (- chirp 1)
(make cylinder (center 0 (* 2 d)) (radius r) (height infinity)
(material (make dielectric (epsilon eps)))))
(geometric-object-duplicates (vector3 -2 0) 0 (- chirp 1)
(make cylinder (center 0 (* 2 d)) (radius r) (height infinity)
(material (make dielectric (epsilon eps)))))
(geometric-object-duplicates (vector3 2 0) 0 (- chirp 1)
(make cylinder (center 0 (* -2 d)) (radius r) (height infinity)
(material (make dielectric (epsilon eps)))))
(geometric-object-duplicates (vector3 -2 0) 0 (- chirp 1)
(make cylinder (center 0 (* -2 d)) (radius r) (height infinity)
(material (make dielectric (epsilon eps)))))
))
(set! pml-layers (list (make pml (thickness dpml))))
(set-param! resolution 20)
(define-param fcen 0.3629) ; pulse center frequency
(define-param df 0.2) ; pulse width (in frequency)
(define-param nfreq 500) ; number of frequencies at which to compute flux
; false = transmission spectrum, true = resonant modes:
(define-param compute-mode? false)
(if compute-mode?
(begin
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Hz) (center 0 0))))
(run-sources+ 400
(at-beginning output-epsilon)
(after-sources (harminv Hz (vector3 0) fcen df)))
(run-until (/ 1 fcen) (at-every (/ 1 fcen 20) output-hfield-z))
)
(begin
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ey)
(center (+ (* -0.5 sx) dpml pad) 0)
(size 0 (* 1 d)))))
(define trans ; transmitted flux
(add-flux fcen df nfreq
(make flux-region
(center (- (* 0.5 sx) dpml pad) 0) (size 0 (* 2 d)))))
(run-sources+ 500
(at-beginning output-epsilon))
(display-fluxes trans) ; print out the flux spectrum
))
fdtd.rar_FDTD透射率_meep_光子晶体_微腔_波导
版权申诉
5星 · 超过95%的资源 179 浏览量
2022-09-24
07:14:02
上传
评论 1
收藏 1KB RAR 举报
小贝德罗
- 粉丝: 69
- 资源: 1万+
最新资源
- 基于matlab实现多车辆车辆路径问题,用遗传算法编程,保证可用.rar
- 基于matlab实现多层极限学习机实现手写体识别,准确率超过99%.rar
- 基于matlab实现电力系统最优潮流程序,可以应用于电力市场下的最优潮流计算 适合电力系统专业的同仁使用.rar
- 串行通信的通信协议串口协议串行通信的通信协议串口协议
- Vue构建交互式的单页面应用程序Vue构建交互式的单页面应用程序
- Spring Boot快速搭建和部署应用程序
- Redis开源的高性键值存储系统广泛应用于缓存、实时消息传递
- Qt跨平台的应用程序开发框架Qt跨平台的应用程序开发框架
- 单片机是集成处理器、存储器和输入输出设备
- 深度学习机器习方法模仿人脑神经网络深度学习机器习方法模仿人脑神经网络
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
评论1