/
chirped_pulse.py
36 lines (27 loc) · 1.1 KB
/
chirped_pulse.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
## linear-chirped pulse planewave with higher frequencies at the front (down-chirp)
import meep as mp
import numpy as np
resolution = 40
dpml = 2
pml_layers = [mp.PML(thickness=dpml,direction=mp.X)]
sx = 40
sy = 6
cell_size = mp.Vector3(sx+2*dpml,sy)
v0 = 1.0 # pulse center frequency
a = 0.2 # Gaussian envelope half-width
b = -0.5 # linear chirp rate (positive: up-chirp, negative: down-chirp)
t0 = 15 # peak time
chirp = lambda t: np.exp(1j*2*np.pi*v0*(t-t0)) * np.exp(-a*(t-t0)**2+1j*b*(t-t0)**2)
sources = [mp.Source(src=mp.CustomSource(src_func=chirp),
center=mp.Vector3(-0.5*sx),
size=mp.Vector3(y=sy),
component=mp.Ez)]
sim = mp.Simulation(cell_size=cell_size,
boundary_layers=pml_layers,
resolution=resolution,
k_point=mp.Vector3(),
sources=sources,
symmetries=[mp.Mirror(mp.Y)])
sim.run(mp.in_volume(mp.Volume(center=mp.Vector3(), size=mp.Vector3(sx,sy)),
mp.at_every(2.7, mp.output_efield_z)),
until=t0+50)