{"id":383257,"date":"2024-09-20T03:24:55","date_gmt":"2024-09-20T03:24:55","guid":{"rendered":"https:\/\/innovationspace.ansys.com\/forum\/forums\/reply\/383257\/"},"modified":"2024-09-20T05:27:49","modified_gmt":"2024-09-20T05:27:49","slug":"383257","status":"publish","type":"reply","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/reply\/383257\/","title":{"rendered":"Reply To: Simulation error in port source"},"content":{"rendered":"

<p> <\/p><p> <\/p><p> <\/p><div>Thank you for your response.<\/div><div>Sorry for asking abstract questions without sending details.<\/div><div>The python script we used for the setup is as follows (I inserted 2d photonic crystals structure as holes in the part that imports the gds files)<\/div><div> <\/div><div><div>import numpy as np<\/div><div>import lumapi<\/div><div>fdtd = lumapi.FDTD()<\/div><div> <\/div><div>#—————set constant——————————————-<\/div><div>a=0.46e-6<\/div><div>b=a*1e6<\/div><div>nx=14<\/div><div>ny=100<\/div><div>G=11<\/div><div> <\/div><div>#—————-set FDTD————————————————-<\/div><div> <\/div><div>fdtd.addfdtd(dimension=”3D”, x = 0,x_span=4*b*nx*1e-6, y = 0, y_min=-20*b\/0.5*1e-6, y_max=ny*b*1e-6 , z = 0.0, z_span=220e-9)<\/div><div>fdtdObj= fdtd.getObjectById(“::model::FDTD”)<\/div><div>fdtdObj.simulation_time=3000e-15<\/div><div>fdtdObj.simulation_temperature=300<\/div><div>fdtdObj.background_material=”<object>”<div>fdtdObj.index=1.0<\/div><div>fdtdObj.x_max_bc=”PML”<\/div><div>fdtdObj.x_min_bc=”PML”<\/div><div>fdtdObj.y_max_bc=”PML”<\/div><div>fdtdObj.y_min_bc=”PML”<\/div><div>fdtdObj.z_max_bc=”PML”<\/div><div>fdtdObj.z_min_bc=”Symmetric”<\/div><div>fdtdObj.same_settings_on_all_boundaries=False<\/div><div> <\/div><div>#——————set mesh, movie monitor, Frequency-domain field and power monitor, SiSlab(rectangle)—————————<\/div><div> <\/div><div>fdtd.addmesh(x = 0,x_span=4*b*nx*1e-6,  y_min=-20*b\/0.5*1e-6, y_max=ny*b*1e-6 , z = 0.0, z_span=220e-9),<\/div><div>fdtd.addmovie(name=”movie_z0″,monitor_type=3,x = 0,x_span=4*b*nx*1e-6,y_min=-20*b\/0.5*1e-6, y_max=ny*b*1e-6 , z = 0.0, lock_aspect_ratio=1, horizontal_resolution=240 ),<\/div><div>fdtd.addpower(name=”monitor_z2″,monitor_type=7,x=0,x_span=2*b*(G-1)*1e-6,<\/div><div>                y=b*(ny-18)*1e-6, y_span=24*b*1e-6, z=1450e-9)<\/div><div>r= fdtd.addrect(name=”SiSlab”, x = 0,x_span=4*b*nx*1e-6, y_min=-20*b\/0.5*1e-6, y_max=ny*b*1e-6 , z = 0.0, z_span=220e-9)<\/div><div>r.material = “Si (Silicon) – Palik”<\/div><div> <\/div><div>#——————–set port——————————————–<\/div><div> <\/div><div>fdtd.addport(),<\/div><div>port1=fdtd.getObjectById(“::model::FDTD::ports::port 1”)<\/div><div>port1.x=0<\/div><div>port1.x_span=5*a\/0.5<\/div><div>port1.y=-8*a\/0.5<\/div><div>port1.z=0<\/div><div>port1.z_span=5<object>*a\/0.5<\/object><\/div><div>port1.direction=”Forward”<\/div><div>port1.amplitude=100<\/div><div>port1.phase=0<\/div><div>port1.mode_selection=”fundamental mode”<\/div><div>port1.number_of_field_profile_samples=1<\/div><div>port1.injection_axis=”y-axis”<\/div><div> <\/div><div># ————set global monitor————————<\/div><div> <\/div><div>fdtd.setglobalmonitor(“use source limits”, False)<\/div><div>fdtd.setglobalmonitor(“minimum wavelength”, 1.25e-6)<\/div><div>fdtd.setglobalmonitor(“maximum wavelength”, 1.45e-6)<\/div><div>fdtd.setglobalmonitor(“frequency points”,41)<\/div><div>fdtd.setglobalmonitor(“use wavelength spacing”,True)<\/div><div> <\/div><div># ————–set global source————————-<\/div><div> <\/div><div>fdtd.setglobalsource(“set wavelength”, True)<\/div><div>fdtd.setglobalsource(“wavelength start”, 1.25e-6)<\/div><div>fdtd.setglobalsource(“wavelength stop”, 1.45e-6)<\/div><div>fdtd.setglobalsource(“eliminate discontinuities”, True)<\/div><div>fdtd.setglobalsource(“optimize for short pulse”, True)<\/div><div>fdtd.setglobalsource(“eliminate dc” , False)<\/div><div> <\/div><div> <\/div><div>#——-import layer from gds files———————————–<\/div><div> <\/div><div>layerlist_1=[0,10,20,30,35]<\/div><div>for layer in layerlist_1:<\/div><div>    fdtd.gdsimport(“gds file 1”, “cell name 1”, layer, “etch”, -110e-9, 110e-9) <\/div><div>layerlist_p2=[41,42]<\/div><div>for layer in layerlist_2:<\/div><div>    fdtd.gdsimport(“gds file 2”, “cell name 2”, layer, “etch”, -110e-9, 110e-9)<\/div><div> <\/div><div> <\/div><div> <\/div><div>fdtd.run()<\/div><div> <\/div><div> <\/div><div> <\/div><div> <\/div><div> <\/div><div> <\/div><div>The setup works fine, but when I run the simulation, nothing is recorded on the monitor.<\/div><div>Note that when I set up and ran the simulation manually, it worked fine. (a=0.5)  I tried to automate it with a python script to change the variable “a” ,but I am facing this problem.<\/div><p> <\/p><br><br>
\n\t\t\t\t\t<\/object><\/div><\/div><p> <\/p><p> <\/p><\/p>\n","protected":false},"template":"","class_list":["post-383257","reply","type-reply","status-publish","hentry"],"aioseo_notices":[],"aioseo_head":"\n\t\t\n\t\n\t\n\t\n\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t