1. 热仿真基础
在微纳光子器件（如微环谐振器、马赫-曾德尔干涉仪、光子晶体腔）中，温度变化会通过热光效应和热膨胀效应显著改变器件的折射率与几何尺寸，进而导致谐振波长漂移。理解三种基本传热机制——热传导、热对流与热辐射——是建立精确热仿真模型的前提。下表总结了三种机制的核心特征与在微纳尺度下的适用性。
| 传热机制 | 能量载体 | 驱动势 | 典型尺度 | 微纳光子器件中的影响 |
|---|---|---|---|---|
| 热传导 | 声子（晶格振动）/ 电子 | 温度梯度 ∇T | 器件内部（μm ~ mm） | 主导波导、衬底内的热扩散；决定热时间常数 |
| 热对流 | 流体分子（空气、水） | 流体速度场 & 温差 | 器件表面（mm ~ cm） | 封装级散热；微流控芯片中显著；真空封装下可忽略 |
| 热辐射 | 光子（电磁波） | T⁴ 温差 | 远场 / 近场（nm ~ m） | 高温（>500 K）或近场辐射时不可忽略；室温下通常次要 |
1.1 热传导：傅里叶定律与微纳尺度修正
热传导是微纳光子器件内部最主要的散热途径。其宏观控制方程为傅里叶定律：
q = -k ∇T其中 q 为热流密度（W/m²），k 为导热系数（W/(m·K)），∇T 为温度梯度。在稳态无内热源条件下，温度场满足拉普拉斯方程 ∇²T = 0。
微纳尺度下的关键修正：
- 尺寸效应：当特征尺寸接近或小于声子平均自由程（硅中约 40 nm @ 300 K）时，边界散射导致等效导热系数显著下降，需使用弹道-扩散模型或分子动力学模拟。
- 界面热阻：不同材料界面（如 Si-SiO₂、Si-空气）存在 Kapitza 热阻，在多层膜波导结构中不可忽略。
- 热各向异性：SOI（Silicon-on-Insulator）衬底中埋氧层导热系数极低（~1.4 W/(m·K)），形成天然热隔离。
仿真提示：在 COMSOL 中设置硅波导热传导时，建议使用“固体传热”物理场，并手动输入温度相关的导热系数 k(T)。对于亚微米波导，可启用“尺寸效应”修正或使用等效各向异性导热系数。
1.2 热对流：牛顿冷却定律与微流控散热
热对流发生在器件表面与周围流体（空气、水、制冷剂）之间。牛顿冷却定律给出：
q = h (T_surface - T_fluid)其中 h 为对流传热系数（W/(m²·K)），其值取决于流体性质、流速和表面几何。在微纳光子器件中，对流通常分为两类：
- 自然对流：h ≈ 5–25 W/(m²·K)，适用于无主动冷却的芯片级封装；散热能力有限。
- 强制对流：h ≈ 50–500 W/(m²·K)，通过微流道、风扇或泵实现；常用于高功率激光器或片上热管理。
注意：在真空封装的光子集成芯片（如 SiPh 收发模块）中，热对流可忽略不计，热传导成为唯一散热路径。此时必须依赖衬底热沉或热电冷却器（TEC）进行温控。
1.3 热辐射：斯特藩-玻尔兹曼定律与近场辐射
热辐射通过电磁波传递能量，不依赖介质。宏观上由斯特藩-玻尔兹曼定律描述：
q = ε σ (T₁⁴ - T₂⁴)其中 ε 为表面发射率（0~1），σ = 5.67×10⁻⁸ W/(m²·K⁴)。在室温（~300 K）下，辐射热流密度通常小于 500 W/m²，远低于传导和对流，因此在大多数微纳光子器件中可忽略。
例外情况：
- 高温器件：如微加热器、热光调制器（工作温度 > 500 K），辐射占比上升。
- 近场辐射：当两个表面间距小于热波长（~10 μm @ 300 K）时，辐射传热可超过黑体极限数个数量级，在近场热光伏、热辅助磁记录中有重要应用。
1.4 三种机制在微环谐振器中的综合影响
下表对比了三种传热机制在典型硅基微环谐振器（半径 10 μm，波导截面 500 nm × 220 nm，SOI 衬底）中的相对贡献：
| 传热路径 | 机制 | 热阻量级 (K/W) | 时间常数 | 对谐振漂移的影响 |
|---|---|---|---|---|
| 波导 → 埋氧层 → 硅衬底 | 传导 | 10³ ~ 10⁴ | μs ~ ms | 主要稳态散热路径 |
| 芯片表面 → 空气 | 自然对流 | 10⁴ ~ 10⁵ | ms ~ s | 封装级温度均匀化 |
| 芯片表面 → 环境 | 辐射 | > 10⁶ | s | 可忽略（室温） |
关键结论：在微环谐振器的热仿真中，热传导是必须精确建模的核心机制；热对流仅在非真空封装或微流控场景下纳入；热辐射在室温标准器件中通常忽略。三种机制的时间常数跨越 6 个数量级，决定了瞬态热仿真的时间步长选择策略。
1.5 热仿真中的多物理场耦合
在微环谐振器的实际仿真中，传热并非孤立物理过程，而是与电磁场、结构力学耦合：
- 热-光耦合：温度变化通过热光系数（dn/dT）改变折射率，导致谐振波长漂移。硅的 dn/dT ≈ 1.86×10⁻⁴ K⁻¹。
- 热-力耦合：热膨胀（CTE ≈ 2.6×10⁻⁶ K⁻¹）改变波导几何尺寸，产生热应力，进一步影响双折射和偏振特性。
- 电-热耦合：微加热器（TiN、Pt 电阻）通过焦耳热产生局部温升，用于热光调谐。
COMSOL 实现建议：使用“电磁热”多物理场接口自动耦合“固体传热”与“波动光学”；在“热膨胀”节点中勾选“热应变”以激活热-力耦合。瞬态仿真时建议将时间步长设为热时间常数的 1/10 ~ 1/5。