Skip to content

Conversation

FilipeFcp
Copy link
Contributor

I believe the results look like the paper, although I had to make some guesses on the material properties and heat power.
@momchil-flex, feel free to take a look if you are interested.

Copy link
Contributor

github-actions bot commented Sep 18, 2025

Spell Check Report

CPOHeat.ipynb:

Cell 1, Line 1: 'co-packaged'
  > # Steady-state and transient thermal analysis for co-packaged optics
Cell 29, Line 16: 'lightgreen'
  > Ax[0].plot(np.log10(time_steps), temp3.values.squeeze(), color="lightgreen", label="ASIC left edge")
Cell 29, Line 28: 'lightgreen'
  > Ax[1].scatter(100, 0, color="lightgreen", s=50, marker="x")
Cell 45, Line 3: 'microrings'
  > This analysis is important for understanding thermal crosstalk between adjacent devices and for defining the requirements for stabilization and control loops that keep the microrings properly tuned.

Checked 1 notebook(s). Found spelling errors in 1 file(s).
Generated by GitHub Action run: https://github.com/flexcompute/tidy3d-notebooks/actions/runs/18017479746

Copy link
Contributor

@tomflexcompute tomflexcompute left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Great result @FilipeFcp and showcase of transient thermal analysis!

  1. I suggest a more formal title such as "Steady-state and transient thermal analysis for co-packaged optics"

  2. Since there are quite a lot of warnings, explain why they are and why it's fine to ignore in this case?

  3. In the temperature plot, plot the structure boundaries as well?

@tomflexcompute
Copy link
Contributor

Also tagging @prashkh since he is interested in this model as well.

Copy link
Collaborator

@momchil-flex momchil-flex left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Looks great and seems to match the paper reasonably well! I wonder if the difference could be coming from the insulating BCs along x and y, in reality there's probably some heat going out of there too?

Just one note: the schematic appears on the top of the first cell but at the end of that cell it says "A schematic of the simulated structure can be seen in the image below:"

Just out of curiosity I noted that the meshing and solver seem to take about the same time ~25 minutes each. Hopefully we'll improve on the meshing side!

image

:maxdepth: 1

../../HeatSolver
../../CPOHeat
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Isn't the usual policy for one notebook to go in one place only?

I'm not necessarily opposed to linking in multiple places just wondering.

Copy link
Contributor

@alec-flexcompute alec-flexcompute left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Excellent work as always @FilipeFcp! A couple suggestions:

  1. You might want to define the unstructured mesh in terms of one of the feature variables. Raw numbers might be more unclear as to why they were chosen.

  2. For the two time_specs, you might want to deifne the time_step and total-time_steps variables at the beginning of cell 71 to have it match cell 79

Copy link
Contributor

@marc-flex marc-flex left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thanks @FilipeFcp. This is an awesome notebook!

Is it possible to use MultiPhysicsMedium instead of Medium to run this notebook? I feel like it would make sense to use it.

@prashkh
Copy link
Contributor

prashkh commented Sep 19, 2025

Very nice notebook! I am sure a lot of people will be interested in this work. Thanks @FilipeFcp !!

@FilipeFcp
Copy link
Contributor Author

Thanks everyone for the comments. Sorry, it took me a while to implement; I was running some mesh tests, and the 2D transient simulation takes quite a long time to run.

With a coarser mesh, the 2D transient simulation shows some strange dips, so I kept the first one.

@momchil-flex, these are the boundaries that they report in the paper. I think the discrepancies are in the source definition. They won't report in volumetric power, and I didn't find it in the paper any information to directly translate it.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

6 participants