SPEED_LIGHT: float = 299_792_458.0 # [m/s] speed of light
## Simulation
# create FDTD Grid
grid = fdtd.Grid(
(2.5e-5, 1.5e-5, 1),
grid_spacing=0.1 * WAVELENGTH,
permittivity=1.0,
permeability=1.0,
)
# sources
grid[50:55, 70:75, 0] = fdtd.LineSource(
period=WAVELENGTH / SPEED_LIGHT, name="linesource"
)
grid[100, 60, 0] = fdtd.PointSource(
period=WAVELENGTH / SPEED_LIGHT, name="pointsource",
)
# detectors
grid[12e-6, :, 0] = fdtd.LineDetector(name="detector")
# x boundaries
# grid[0, :, :] = fdtd.PeriodicBoundary(name="xbounds")
grid[0:10, :, :] = fdtd.PML(name="pml_xlow")
grid[-10:, :, :] = fdtd.PML(name="pml_xhigh")
# y boundaries
# grid[:, 0, :] = fdtd.PeriodicBoundary(name="ybounds")
plt.plot(np.fft.fft(E.real))
#%%
import fdtd
grid = fdtd.Grid(
shape=(25e-6, 15e-6, 1), # 25um x 15um x 1 (grid_spacing) --> 2D FDTD
)
print(grid)
grid[11:32, 30:84, 0] = fdtd.Object(permittivity=1.7**2, name="object")
grid[13e-6:18e-6, 5e-6:8e-6, 0] = fdtd.Object(permittivity=1.5**2)
grid[7.5e-6:8.0e-6, 11.8e-6:13.0e-6,
0] = fdtd.LineSource(period=1550e-9 / (3e8), name="source")
grid[12e-6, :, 0] = fdtd.LineDetector(name="detector")
# x boundaries
grid[0:10, :, :] = fdtd.PML(name="pml_xlow")
grid[-10:, :, :] = fdtd.PML(name="pml_xhigh")
# y boundaries
grid[:, 0:10, :] = fdtd.PML(name="pml_ylow")
grid[:, -10:, :] = fdtd.PML(name="pml_yhigh")
grid.run(total_time=10)
grid.visualize(z=0)
WAVELENGTH = 1550e-9
SPEED_LIGHT: float = 299_792_458.0 # [m/s] speed of light
## Grid setup
# create FDTD Grid
grid = fdtd.Grid(
(1.5e-5, 1.5e-5, 1), # 2D grid
grid_spacing=0.1 * WAVELENGTH,
permittivity=2.5, # same as object
)
# sources
grid[15, :] = fdtd.LineSource(period=WAVELENGTH / SPEED_LIGHT, name="source")
# detectors
grid[-15, :, 0] = fdtd.LineDetector(name="detector")
# x boundaries
# grid[0, :, :] = fdtd.PeriodicBoundary(name="xbounds")
grid[0:10, :, :] = fdtd.PML(name="pml_xlow")
grid[-10:, :, :] = fdtd.PML(name="pml_xhigh")
# y boundaries
# grid[:, 0, :] = fdtd.PeriodicBoundary(name="ybounds")
grid[:, 0:10, :] = fdtd.PML(name="pml_ylow")
grid[:, -10:, :] = fdtd.PML(name="pml_yhigh")
# The absorbing object fills the whole grid