Hello!
I am trying to calculate field energy in air before the field passes through
a medium and after it passes through. The medium has a refractive index of
2.1609.
I have used two cases: with periodic boundary conditions (using k-point) and
without it.
Also I have put metal nano cylinder and this leads us to a total of four
cases as given below.
nc- nano cylinder
kp -k-point
One - field -energy after the medium
two - field - energy before the medium
the field energy computations were done at the last time step (after third
time step in this case)
1. *no nc and yes kp* : One-3.34630186473115e-5 two-3.5681340908767e-5
2. *yes nc and no kp* : One-6.37494514942021e-8 two-3.80564486417271e-5
3. *yes nc and yes kp* : One-1.75693509082845e-6 two-6.01227187135615e-5
4. *no nc and no kp* : One-3.05763975245008e-5 two-2.52388151348767e-5
Question:
How are the results getting affected by k-pojnt condition?
Especially if we compare 'yes nc' case then there is a difference of 3
orders of magnitude!! My question is that is my k-point condition correct in
the program?
-----------------------------------------------------------------------------------------------------------
CTL FILE:
(set! geometry-lattice (make lattice (size .5 .09 no-size)))
(set! geometry (list
(make block (center 0 0) (size .3 .09 0.4)
(material(make dielectric (epsilon 2.1609))))))
;(make cylinder (center 0 0) (radius .035) (height .110)
;(material(make dielectric (epsilon 2.16091)
; (polarizations
;(make polarizability
;(omega 1e-30) (gamma .1561792)
;(delta-epsilon 41.9731539e+60))))))))
(set! k-point (vector3 1 1 0))
(set! sources (list
(make source
(src (make continuous-src (frequency 2.1777)))
(direction X)
(component Hz)
(center -.2 0)
(size 0 .09 0))))
(set! pml-layers (list (make pml (thickness 0.025)(direction X))))
(set! resolution 1000)
(run-until 2
(at-beginning output-epsilon))
(run-until 1
(at-every .5 output-efield-x)
(at-every .5 output-efield-y))
(print "One" (field-energy-in-box (volume (center 0.15 0) (size 0.01 0.01
0))) " "
"two" (field-energy-in-box (volume (center -0.15 0) (size 0.01 0.01 0))))
----------------------------------------------------------------------------------------------------------------
Thanks!
regards
--
Ajay Ram S K S M
MSc(Physics)
Sri Sathya Sai University
Prasanthi Nilayam
India.
_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
