Here in the following you can find MY code for a time changing material.
You can find different solutions in the web.
import meep as mp
import numpy as np
resolution = 50
cell_size = mp.Vector3(14,10,0)
pml_layers = [mp.PML(thickness=2,direction=mp.X)]
rot_angle = np.radians(0)
fsrc = 1.0
eps1=1
eps2=10
default_material1 = mp.Medium(epsilon=eps1)
k_point1 = mp.Vector3(fsrc*np.sqrt(eps1)).rotate(mp.Vector3(z=1),rot_angle)
sources1= [mp.EigenModeSource(src=mp.ContinuousSource(fsrc),
center=mp.Vector3(-5,0),
size=mp.Vector3(y=10),
direction=mp.AUTOMATIC if rot_angle == 0 else
mp.NO_DIRECTION,
eig_kpoint=k_point1,
eig_band=1,
eig_parity=mp.EVEN_Y+mp.ODD_Z if rot_angle ==
0 else mp.ODD_Z,
eig_match_freq=True)]
sim1 = mp.Simulation(cell_size=cell_size,
resolution=resolution,
boundary_layers=pml_layers,
sources=sources1,
k_point=k_point1,
default_material=default_material1,
symmetries=[mp.Mirror(mp.Y)] if rot_angle == 0 else []
)
def cng_mat(sim):
mt=sim.meep_time()
if sim.meep_time() > 15:
sim.set_materials(default_material=mp.Medium(epsilon=eps2))
else :
sim.set_materials(default_material=mp.Medium(epsilon=eps1))
sim1.use_output_directory()
sim1.run(mp.at_every(0.1,cng_mat,
mp.in_volume(mp.Volume(mp.Vector3(0,0),size=mp.Vector3(0,0)),mp.to_appended("epsT",mp.output_epsilon)),
mp.in_volume(mp.Volume(mp.Vector3(0,0),size=mp.Vector3(10,4)),mp.to_appended("ezT",mp.output_efield_z))),
until=30)
Enjoy...
Massimo
Il giorno sab 9 nov 2019 alle ore 16:58 Cícero Julião <
juliaojun...@gmail.com> ha scritto:
> Anyone have a working code with epsilon changing during simulation?
>
> Thank you.
>
> --
>
> Cícero Julião da Silva Junior
> ____________________________________
> Instituto Federal de Alagoas - IFAL
> Campus Murici
> http://www.ifal.edu.br/
> _______________________________________________
> meep-discuss mailing list
> meep-discuss@ab-initio.mit.edu
> http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
--
--
Massimo Moccia
Fields & Waves Lab
Department of Engineering - University of Sannio
Palazzo Dell'Aquila Bosco Lucarelli
Corso Garibaldi 107, I-82100 Benevento, Italy
Voice: +39 0824 305809
Email: massimo.moc...@unisannio.it
Skype: max.moccia
URL: http://fw-lab.org
_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
