Learning Forum

[Leon Neymeyer]

Hi,

I am importing temperature data from Lumerical HEAT to simulate a thermo-optical effect. I used your instructions to do this. The only difference is that I don't sweep over the power, but give the electrode a temperature as a boundary condition. I also sweep over this which works fine. Now I import the data and save the effective refractive index.

I calculate the phase shift as follows
phase = 2*pi*(neff-neff0)*L/lambda;

The first temperature value from the thermal simulation is 293.15K, i.e. room temperature. I would expect neff and neff0 to be equal and the phase shift to be 0. Instead I get a value not equal to 0. I have now changed the "reference temperature" in the material database which did not fixed the problem.

I would expect 0 phaseshift at the first value:

Here is my code:

# Temperaturdaten importieren ########################

temp = linspace(293.15,473.15,7);

neff = matrix(length(temp));

 

addgridattribute("temperature");

matlabload("T_temp_sweep.mat");

select('temperature');

importdataset(T);

 

# Brechungsindizes berechnen

switchtolayout;

setnamed('temperature','enabled',0);

setanalysis('wavelength',1.56e-6);

findmodes;

neff0 = real(getdata('mode1','neff'));

 

for (i=1:length(temp)){

switchtolayout;

setnamed('temperature','enabled',1);

setnamed('temperature','Temperatur_index',i);

findmodes;

neff(i) = real(getdata('mode1','neff'));

}

 

L_long = 6760e-6;

L_short = 4327e-6;

L_avg = 5543e-6;

lambda = 1.56e-6;

phase_long = 2*pi*(neff-neff0)*L_long/lambda;

phase_short = 2*pi*(neff-neff0)*L_short/lambda;

phase_avg = 2*pi*(neff-neff0)*L_avg/lambda;

 

Thanks in advance!

Leon

[Leon Neymeyer]

Pretty simple solution. I just had to set the simulation temperature in FDE.

[Ashish Khemka]

Thanks for sharing the answer.