test_absorbing_media.py

import numpy as np
from simulator import Simulation1D, um, fs


K = 200  #Number of elements
zpos_vertex = np.linspace(0,10*um,K+1) # z position of the vertices of the elements
Np = 4 #polynomial order

#medium parameters
ε = np.ones(K)
𝜇 = np.ones(K) 
σ = np.ones(K)*0.01
sim = Simulation1D(zpos_vertex, Np ,  ε, 𝜇 , σ,
                   boundaries = 'perfect_conductor_boundaries')





#Set initial conditions
def intial_fields(z):
    sig    = 0.5*um
    z0 = 5.0*um

    Ex = np.exp(-1/sig**2*((z -z0))**2)
    Hy = Ex*0
    return Ex, Hy


sim.set_initial_conditions(intial_fields)


#plot fields at t = 0 
sim.plot_fields(title = "t = 0 fs")



# run the simulation 
total_time   = 3 * fs
nt = int(500 * total_time)
dt = total_time/nt
sim.run(total_time, dt)
sim.plot_fields(title = "t = 3 fs")



# run the simulation 
total_time   = 4 * fs
nt = int(500 * total_time)
dt = total_time/nt
sim.run(total_time, dt)
sim.plot_fields(title = "t = 6 fs")