Ansys Learning Forum Maintenance

NOTICE: We will be performing backend maintenance of our Learning Forum from April 5 to April 12, 2021. The result will be a new infrastructure but with little impact to user experience and design. Currently the forum is accessible in read-only mode as we complete our final migration. Thank you for your patience. For urgent issues please visit HERE.

Thermally tuned Waveguide tutorial trouble — Ansys Learning Forum

Thermally tuned Waveguide tutorial trouble

mattsarg55mattsarg55 Member Posts: 14

Hello fellow lumerical users,

I am new to the software and was attempting to do the thermally tuned waveguide tutorial found here:

I was able to complete the tutorial succesfully but I am attempting to apply the logic found here to my own waveguide. I create the entire structure in HEAT, run a sweep, download the temperature grid, and upload that info to MODE to get this (Note: I only included the waveguide and a bit of the surrounding material in MODE which is why my temp grid structure is much longer):

However, once I begin running the script provided in the tutorial:


My temp grid shrinks back to the origin and I see this:

an I still get results for my sweeps but they are all zero. Why is this happening? Do I need to change something in the provided script? Also, where does the equation, Phase = 2*pi*(neff-neff0)*L/lambda come from? Does it have a name?

Hopefully someone can help me with this as lumerical seems to be an amazing program and I want to know how to apply it to the work at my job


  • gsungsun Ansys Employee Posts: 227

    This script file does not set simulation region, nor the attribute size. so you will need to check your FDE settings. Before doing sweep, try to calculate modes at one temperature, by assigning temperature index, and see of the simulation region is, you add a "break" after this section:




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


    and see if it works as expected.

    "where does the equation, Phase = 2*pi*(neff-neff0)*L/lambda come from? Does it have a name?"

    As you can see, the script file gives neff0 in the above script, which is a reference point. At other temperature, the neff will be different from it, and thus creates a phase change.

    Please test step by step.

  • mattsarg555mattsarg555 Member Posts: 7

    Hello! Thank you for the response! I went through it and it is making much more sense!

    I am still wondering about the equation, I understand the effective indexes and how we are getting them but what about L? I know we assign its value but do I need to change this for my own experiments? What does the variable 'L' in the script represent?

  • gsungsun Ansys Employee Posts: 227

    L is the length of the device, so 2*pi*(neff-neff0)*L/lambda gives a phase difference to the reference temperature. Think that a waveguide mode with neff,and neff0 will create a phase difference through a propagation length of L.

  • mattsarg555mattsarg555 Member Posts: 7

    Gotca! Thank you gsun! Finally, is the phase difference in degrees or radians?

  • gsungsun Ansys Employee Posts: 227

    it uses pi, so it is radian. You can change it to degree if you want.

  • mattsarg555mattsarg555 Member Posts: 7

    Okay! So i may have said that I understood your explanation of the variable 'L' without actually looking at the model. If we look at the model used in HEAT in the example given:

    the length of this object in the x is 100um and in the y direction is 50um. When we say the propagation length, are we refering to the 'z-direction'?

  • gsungsun Ansys Employee Posts: 227

    yes, you are right. it is the propagation length. Phase can only be accumulated through propagation for regular waveguide.

  • mattsarg555mattsarg555 Member Posts: 7

    Great! Glad I understood that, my final question is one that you may or may not be able to answer... in the example given we find a phase shift of 128 degrees(roughly 2.2 radians)and that seems to be quite large! I am new to photonics so having a shift be that large at such small input powers seems destructive and I would expect it instead to be on the order of only a few degrees. Why is this? If you cannot answer this then I totally understand as it clearly deals more with optical theory and not the software.

  • gsungsun Ansys Employee Posts: 227
    edited April 5

    Since this is a thermally modulated device, so the result will depend on the temperature. In addition, your simulation can be different from real experiment, since simulation region is limited, and you may use a constant power source and insulate boundaries. With time, if the convection is less than the source power, the temperature will constantly increase. So you may need to check the max temperature and see if it is reasonable. if it is too high , it might mean that your boundary conditions are improper.

    In addition, the phase change can be purposely to be as large as 180 degree, depending on the device its self. This will be thermal modulation.

    I hope this helps.

Sign In or Register to comment.