Hi, Steven and Meep users:
I tried to use meep to calculate the decay rate of a single fibre with
absord. First I calculated the energy in a fibre without absord, but each
time the result was not stable, as shown in the following.
I wonder if there is anything wrong with the .ctl file or there is something
else that I should care about.
Another question, as you see in my .ctl file, the volume for calculating the
energy was defined in the 3D dimension, so the cross section is a square,
which is a little bigger than the cross section of the fibre. I also wonder
if there is any way to define a volume with a round cross section.
I also want to know how to define the volume if I set the "dimensions" as
the "CYLINDRICAL" .(volume (center 0 0 0) (r 0.5) (z 0)), just like this? or
in another way?
Thank you.
L.H Zhao
&&&&&&&&&&&&&
Here is the last few values for the energy:
field-energy: 1.15191894175373
field-energy: 1.83168135254997
field-energy: 3.23270907203186
field-energy: 1.13683665930588
field-energy: 1.82739746343853
field-energy: 3.23957118549375
field-energy: 1.14224795599941
field-energy: 1.82999995502438
field-energy: 3.2436454407705
field-energy: 1.14772292007047
field-energy: 1.82536789401917
field-energy: 3.22532329952147
.ctl file:
(set! geometry-lattice (make lattice (size 10 5 5)))
(define Si (make dielectric (index 3.2)))
(set! geometry (list
(make cylinder (center 0 0 0) (radius 0.25) (height infinity)
(axis 1 0 0) (material Si))))
(set! sources (list
(make source
(src (make continuous-src (wavelength 1.5)))
(component Ez)
(center -4 0 0))))
(set! pml-layers (list (make pml (thickness 1.0))))
(set! resolution 20)
(define (print-field-energy)
(print "field-energy: " (field-energy-in-box (volume (center 0 0 0)
(size 0 0.5 0.5))) "\n"))
(run-until 1000
(at-beginning output-epsilon)
(at-end output-efield-z)
(at-every 10 print-field-energy))
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
