Folks, I just started using MPB to study various rectangular dielectric waveguide geometries. I’m interested in looking at the portion of the field energy in a rectangular core and a rectangular cladding. When I use compute-energy-in-objects, it curiously returns more than 100%.
I’m using:
mpbi 1.7.0, Copyright (C) 1999-2012, MIT
Using libctl 4.2.0 and Guile 3.0.2.
On an Ubuntu installation.
When I’m looking at just the first mode:
(get-dfield 1) ; get d-field for first band
(compute-field-energy)
(print
"Gamma_guide: "
(compute-energy-in-objects (make block (size infinity w h) (center 0 0 0)
(material air)))
"\n")
(print
"Gamma_clad: "
(compute-energy-in-objects
(make block (size infinity (+ w (* 2 tClad)) (+ h (* 2 tClad)))
(center 0 0 0) (material air)))
"\n")
(print
"Gamma_cell: "
(compute-energy-in-objects (make block (size infinity sc-z sc-y) (center 0 0
0) (material air)))
"\n")
I received the below output:
D-energy-components:, 102, 1, 0.0919373, 0.907724, 0.000338274
$35 = (0.9999999999997949 3.843388742030928e-9 0.09193731941226466
0.9077244029679912 9.523365114586835e-11 3.382689149641591e-4
4.766081590605888e-9)
Gamma_guide: 1.085048755617456
Gamma_clad: 1.0923722948863863
Gamma_cell: 1.0923723008598665
What am I missing?
-JW
smime.p7s
Description: S/MIME cryptographic signature
_______________________________________________ mpb-discuss mailing list mpb-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
