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
|
import meep as mp
cell_size = mp.Vector3(6, 6, 0)
geometry1 = [
mp.Cylinder(center=mp.Vector3(), radius=1.0, material=mp.Medium(index=3.5))
]
sim1 = mp.Simulation(cell_size=cell_size, geometry=geometry1, resolution=20)
sim1.init_sim()
geometry2 = [
mp.Cylinder(center=mp.Vector3(1, 1), radius=1.0, material=mp.Medium(index=3.5))
]
sim2 = mp.Simulation(cell_size=cell_size, geometry=geometry2, resolution=20)
sim2.init_sim()
sim1.fields.phase_in_material(sim2.structure, 10.0)
sim1.run(
mp.at_beginning(mp.output_epsilon), mp.at_every(0.5, mp.output_epsilon), until=10
)
|
