-
Notifications
You must be signed in to change notification settings - Fork 613
/
ring-mode-overlap.py
67 lines (51 loc) · 1.83 KB
/
ring-mode-overlap.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
# Calculating 2d ring-resonator modes, from the Meep tutorial.
import meep as mp
n = 3.4 # index of waveguide
w = 1 # width of waveguide
r = 1 # inner radius of ring
pad = 4 # padding between waveguide and edge of PML
dpml = 2 # thickness of PML
sxy = 2 * (r + w + pad + dpml) # cell size
cell = mp.Vector3(sxy, sxy)
# Create a ring waveguide by two overlapping cylinders - later objects
# take precedence over earlier objects, so we put the outer cylinder first.
# and the inner (air) cylinder second.
geometry = [
mp.Cylinder(radius=r + w, height=mp.inf, material=mp.Medium(index=n)),
mp.Cylinder(radius=r, height=mp.inf, material=mp.air),
]
pml_layers = [mp.PML(dpml)]
resolution = 20
# If we don't want to excite a specific mode symmetry, we can just
# put a single point source at some arbitrary place, pointing in some
# arbitrary direction. We will only look for Ez-polarized modes.
fcen = 0.118 # pulse center frequency
df = 0.010 # pulse width (in frequency)
sources = [
mp.Source(
src=mp.GaussianSource(fcen, fwidth=df),
component=mp.Ez,
center=mp.Vector3(r + 0.1),
)
]
# exploit the mirror symmetry in structure+source:
symmetries = [mp.Mirror(mp.Y)]
sim = mp.Simulation(
cell_size=cell,
resolution=resolution,
geometry=geometry,
boundary_layers=pml_layers,
sources=sources,
symmetries=symmetries,
)
h1 = mp.Harminv(mp.Ez, mp.Vector3(r + 0.1), fcen, df)
sim.run(mp.after_sources(h1), until_after_sources=300)
fields2 = sim.fields
sim.reset_meep()
fcen = 0.236
h2 = mp.Harminv(mp.Ez, mp.Vector3(r + 0.1), fcen, df)
sim.run(mp.after_sources(h2), until_after_sources=300)
def overlap_integral(r, ez1, ez2):
return ez1.conjugate() * ez2
res = sim.integrate2_field_function(fields2, [mp.Ez], [mp.Ez], overlap_integral)
print(f"overlap integral of mode at w and 2w: {abs(res)}")