Hi every body,
I am trying to find the flux spectrum between two parallel plates caused by
thermal emission. In fact, I am going to reproduce the results of the
following paper (red dashed line on Fig.2 bottom):
Frequency-Selective Near-Field Radiative Heat Transfer between Photonic
Crystal Slabs: A Computational Approach for Arbitrary Geometries and
Materials - By: Alejandro W. Rodriguez et al
At the end my written CTL file is attached. I have few questions regarding
to this:
1- In this paper a Drude material is used for plates using this
formula:
Epsilon(omega)=epsilon_inf-sigma/(omega*(omega+i*gamma))
epsilon_inf=12.5
sigma=0.2533(2*pi*c/a)^2 "a" is unit length (here taken as 1e-6m)
gamma=1.5915(2*pi*c/a)
I used "noisy-drude-susceptibility" for hot plate with:
Frequency=1
Gamma=1.5915*2*pi=10
Sigma=0.2533*(2*pi)^2=10
And for cold plate I used "drude-susceptibility" with same numbers of
frequency, gamma and sigma. Is this definition correct?
2-in order to set Bloch periodicity along x-direction with the phase
difference of exp(i*kx*a), I set the k-point as (vector3 1 0 0). Does it
mean the same concept?
3-After doing these, I set a point "gaussian-src" source at the
center of the hot plate with Ex component and I add flux region along the
center line of the cold plate to calculate the flux (which is Phi in the
paper on figure 2, bottom).
Doing all these, I got wrong results. Actually, my flux spectrum is very
sensitive to resolution and calculation time. I am really confused. Can
anybody help me?
Thanks
CTL file ------------------------------------------------------------
(define-param pml_t 0.1) ;thickness of the PML layer
(define-param h 0.2) ;thickness of the plates
(define-param d 0.2) ;gap distance
(define-param sx 1) ; length of the plates along x-direction
(define-param SX (+ sx (* 2 pml_t)))
(define-param sy (+ h h d))
(define-param SY (+ sy (* 2 pml_t)))
(define-param res 85) ;resolution
(define-param time 180) ;calculation time
(define-param amp (/ 1 (sqrt time))) ;calculating the noise-amp
(set! geometry-lattice (make lattice (size sx SY no-size)))
(set! pml-layers (list (make pml (thickness pml_t) (direction Y))))
(set-param! resolution res)
(set! k-point (vector3 1 0 0)) ;periodicity along x. Bloch-periodicity
(set! ensure-periodicity true)
;setting the material of the hot plate
(define hot (make dielectric (epsilon 12.5) ;reference length is 1e-6 m
(E-susceptibilities
(make noisy-drude-susceptibility
(noise-amp amp) (frequency 1) (gamma 10) (sigma 10))
)))
;setting the material of the cold plate
(define cold (make dielectric (epsilon 12.5) ;reference length is 1e-6 m
(E-susceptibilities
(make drude-susceptibility
(frequency 1) (gamma 10) (sigma 10))
)))
(define-param fcen 0.095)
(define-param df 0.191)
(define-param fnum 500)
(set! geometry (list
(make block (center 0 (* -1 (+ (/ d 2) (/ h 2))) 0) (size sx h infinity)
(material hot ))
(make block (center 0 (+ (/ d 2) (/ h 2)) 0) (size sx h infinity)
(material cold )) ))
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df))) ;
continuous-src
(component Ex) ;(size sx h)
(center 0 (* -1 (+ (/ d 2) (/ h 2))) 0)) ))
(define trans ; transmitted flux
(add-flux fcen df fnum
(make flux-region (center 0 (* 1 (+ (/ d 2) (/ h 2))) 0) (size sx 0 0))
))
;(run-until 100 (at-every 0.1 (output-png Ez "-Zc bluered")))
;(run-until 1000 (at-beginning output-epsilon))
(run-until time)
(display-fluxes trans)
_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
