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

Attachment: 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

Reply via email to