Hello ANSYS Support Team,
I am running my simulation using the C: drive, which has plenty of free space,also i checked power and temperature everything. However, during my inverse design simulation in ANSYS, the computer shuts down automatically after about 10–12 steps.
Previously, the same code ran successfully up to 220-250 steps without any issues.
Also I am attaching my code
base diamond Code
switchtolayout;
selectall;
delete;
## SIMULATION PARAMETERS
# Waveguide and Geometry Parameters
num_wg = 2;
wg_width = 2.0e-6;
wg_height = 1.0e-6;
out_wg_dist = 1e-6;
total_wg_h_span = num_wg*wg_width + (num_wg-1)*out_wg_dist;
# Optimization and Simulation Region Size
opt_size_x = 4e-6;
opt_size_y = 5e-6;
opt_size_z = wg_height;
size_x = opt_size_x + 1e-6;
size_y = opt_size_y + 1e-6;
size_z = 2.2e-6;
# Material and Mesh Parameters
wg_index = 2.417;
bg_index = 1.46;
fine_mesh_size = 20e-9;
## GEOMETRY SETUP
# Input Waveguide
addrect;
set(‘name’,’input wg’);
set(‘x min’,-size_x);
set(‘x max’,-opt_size_x/2 + 5e-7);
set(‘y’,0);
set(‘y span’,wg_width);
set(‘z’, 0);
set(‘z span’, wg_height);
set(‘index’,wg_index);
# Output Waveguides
for( i=1:num_wg ) {
addrect;
set(‘name’,’output wg ‘+num2str(i));
set(‘x min’,opt_size_x/2 – 5e-7);
set(‘x max’,size_x);
set(‘y’,-total_wg_h_span/2 + (wg_width/2) + (i-1)*(out_wg_dist+wg_width));
set(‘y span’,wg_width);
set(‘z’, 0);
set(‘z span’, wg_height);
set(‘index’,wg_index);
}
## SOURCE SETUP
addmode;
set(‘direction’,’Forward’);
set(‘injection axis’,’x-axis’);
set(‘x’,-size_x/2+0.2e-6);
set(‘y’,0);
set(‘y span’,1.5*wg_width);
set(‘z’,0);
set(‘z span’,size_z – 0.2e-6);
set(‘center wavelength’,500e-9);
set(‘wavelength span’,0);
set(‘mode selection’,’fundamental TE mode’);
## FDTD SIMULATION REGION (FDTD)
addfdtd;
set(‘dimension’,’3D’);
set(‘index’,bg_index);
set(‘mesh accuracy’,3);
set(‘x min’,-size_x/2);
set(‘x max’,size_x/2);
set(‘y min’,-size_y/2);
set(‘y max’,size_y/2);
set(‘z min’,-size_z/2);
set(‘z max’,size_z/2);
set(‘z min bc’,’PML’);
set(‘z max bc’,’PML’);
set(‘auto shutoff min’,1e-6);
set(‘simulation time’,5000e-15);
## MESH REFINEMENT FOR 100nm FEATURES
addmesh;
set(‘name’,’opt_mesh’);
set(‘x’,0);
set(‘x span’,opt_size_x);
set(‘y min’,-opt_size_y/2);
set(‘y max’,opt_size_y/2);
set(‘z’, 0);
set(‘z span’, opt_size_z);
set(‘override x mesh’,true);
set(‘dx’,fine_mesh_size);
set(‘override y mesh’,true);
set(‘dy’,fine_mesh_size);
set(‘override z mesh’,true);
set(‘dz’,fine_mesh_size);
## MONITORS
# Optimization Fields Monitor
adddftmonitor;
set(‘name’,’opt_fields’);
set(‘monitor type’,’3D’);
set(‘x’,0);
set(‘x span’,opt_size_x);
set(‘y min’,-opt_size_y/2);
set(‘y max’,opt_size_y/2);
set(‘z’, 0);
set(‘z span’, opt_size_z);
# FOM (Figure of Merit) Monitors
for( i=1:num_wg ) {
y_center = -total_wg_h_span/2 + (wg_width/2) + (i-1)*(out_wg_dist+wg_width);
adddftmonitor;
set(‘name’,’fom_’+num2str(i));
set(‘monitor type’,’2D X-normal’);
set(‘x’, size_x/2-0.2e-6);
set(‘y’, y_center);
set(‘y span’, 1.5*wg_width);
set(‘z min’, -size_z/2);
set(‘z max’, size_z/2);
}
# Global Index Monitor (for visualization)
addindex;
set(‘name’,’global_index’);
set(‘monitor type’,’3D’);
set(‘x min’,-size_x/2);
set(‘x max’,size_x/2);
set(‘y min’,-size_y/2);
set(‘y max’,size_y/2);
set(‘z min’,-size_z/2);
set(‘z max’,size_z/2);
#INITIAL GUESS
addstructuregroup;
set(“name”,”initial_guess”);
for( i=1:num_wg ) {
y_out = -total_wg_h_span/2 + (wg_width/2) + (i-1)*(out_wg_dist+wg_width);
addwaveguide;
set(“name”,”bend”+num2str(i));
set(“base width”, wg_width);
set(“base height”, wg_height);
set(“base angle”,90);
poles = [-opt_size_x/2, 0;
0, 0;
0, y_out;
opt_size_x/2, y_out];
set(“poles”,poles);
set(“index”,wg_index);
addtogroup(“initial_guess”);
}
python code
######## IMPORTS ########
# General purpose imports
import os
import sys
import numpy as np
# Optimization specific imports
from lumopt import CONFIG
from lumopt.geometries.topology import TopologyOptimization2D,TopologyOptimization3DLayered
from lumopt.utilities.load_lumerical_scripts import load_from_lsf
from lumopt.figures_of_merit.modematch import ModeMatch
from lumopt.optimization import Optimization
from lumopt.optimizers.generic_optimizers import ScipyOptimizers
from lumopt.utilities.wavelengths import Wavelengths
######## DEFINE BASE SIMULATION ########
def runSim(params, eps_bg, eps_wg, x_pos, y_pos,z_pos, size_x, size_y,size_z, filter_R, min_feature_size, working_dir, beta):
######## DEFINE GEOMETRY ########
geometry =TopologyOptimization3DLayered(params=params, eps_min=eps_bg, eps_max=eps_wg, x=x_pos, y=y_pos,z=z_pos, filter_R=filter_R, min_feature_size=min_feature_size, beta=beta)
######## DEFINE FIGURE OF MERIT FOR EACH OUTPUT WAVEGUIDE ########
fom1 = ModeMatch(monitor_name = ‘fom_1’, mode_number = ‘Fundamental TE mode’, direction = ‘Forward’, norm_p = 2, target_fom=1)
fom2 = ModeMatch(monitor_name = ‘fom_2’, mode_number = ‘Fundamental TE mode’, direction = ‘Forward’, norm_p = 2, target_fom=1)
######## DEFINE OPTIMIZATION ALGORITHM ########
optimizer = ScipyOptimizers(max_iter=200, method=’L-BFGS-B’, pgtol=1e-6, ftol=1e-5, scale_initial_gradient_to=0.25)
######## DEFINE SETUP SCRIPT AND INDIVIDUAL OPTIMIZERS ########
script = load_from_lsf(‘C:/Users/lumer/Desktop/Lumerical/base_diamond.lsf’)
script = script.replace(‘opt_size_x=4e-6′,’opt_size_x={:1.6g}’.format(size_x))
script = script.replace(‘opt_size_y=5e-6′,’opt_size_y={:1.6g}’.format(size_y))
wavelengths1 = Wavelengths(start = 520e-9, stop = 540e-9, points = 11)
opt1 = Optimization(base_script=script, wavelengths = wavelengths1, fom=fom1, geometry=geometry, optimizer=optimizer, use_deps=False, hide_fdtd_cad=True, plot_history=False, store_all_simulations=False, save_global_index=True)
wavelengths2 = Wavelengths(start = 600e-9, stop = 800e-9, points = 11)
opt2 = Optimization(base_script=script, wavelengths = wavelengths2, fom=fom2, geometry=geometry, optimizer=optimizer, use_deps=False, hide_fdtd_cad=True, plot_history=False, store_all_simulations=False, save_global_index=False)
######## PUT EVERYTHING TOGETHER AND RUN ########
opt = opt1+opt2
opt.run(working_dir = working_dir)
if __name__ == ‘__main__’:
size_x = 4000
size_y = 5000
size_z = 1000
filter_R = 75e-9
min_feature_size=100e-9
x_points=int(size_x/20)+1
y_points=int(size_y/20)+1
z_points=int(size_z/20)+1
eps_wg = 2.417**2
eps_bg = 1.46**2
x_pos = np.linspace(-size_x/2*1e-9,size_x/2*1e-9,x_points)
y_pos = np.linspace(-size_y/2*1e-9,size_y/2*1e-9,y_points)
z_pos = np.linspace(-size_z/2*1e-9,size_z/2*1e-9,z_points)
## We need to pick an initial condition. Many different options:
params = 0.5*np.ones((x_points,y_points)) #< Start with the domain filled with (eps_wg+eps_bg)/2
#params = np.ones((x_points,y_points)) #< Start with the domain filled with eps_wg
#params = np.zeros((x_points,y_points)) #< Start with the domain filled with eps_bg
#params = None #< Use the structure defined in the project file as initial condition
working_dir = ‘diamond_3D_TE_topology_x{:04d}_y{:04d}_f{:04d}’.format(size_x,size_y,int(filter_R*1e9))
runSim(params, eps_bg, eps_wg, x_pos, y_pos,z_pos, size_x*1e-9, size_y*1e-9,size_z*1e-9,filter_R, min_feature_size,working_dir=working_dir, beta=1)
Could you please advise on what might be causing this and how I can resolve it?
You are navigating away from the AIS Discovery experience
