4. 温度均匀性指标：定义与计算方法

在MZI结构的热光调谐与热沉布局优化中，温度场的均匀性直接决定了多通道干涉相位的控制精度。本章建立一套量化的温度均匀性评价体系，包含温度梯度、最大温差（ΔT_max）和均方根误差（RMSE）三个核心指标，用于后续仿真结果的横向对比与布局方案筛选。

4.1 指标定义与物理意义

4.2 温度梯度：空间热流的方向标

温度梯度是矢量量，在COMSOL后处理中通过 d(T,x)、d(T,y)、d(T,z) 算子提取。对于MZI结构，重点关注波导传播方向（x方向）和垂直于波导的方向（y方向）的梯度分量。

• 纵向梯度（∂T/∂x）：影响光波沿传播方向的相位累积均匀性，梯度过大会导致两臂等效折射率沿程变化不一致。
• 横向梯度（∂T/∂y）：反映MZI两臂之间的温度差异，是造成两臂相位差的主要来源。
• 法向梯度（∂T/∂z）：在薄层SOI结构中通常较小，但在多层热沉界面处不可忽略。

4.3 最大温差（ΔT_max）：极差控制红线

ΔT_max 是工程中最直观的均匀性指标。在COMSOL中，通过 maxop1(T) - minop1(T) 积分算子计算，其中 maxop1 和 minop1 定义在MZI波导区域边界上。

// COMSOL 后处理表达式示例（定义在MZI波导域）
// 最大温差
DeltaT_max = maxop1(T) - minop1(T)

// 归一化最大温差（用于不同加热功率下的对比）
DeltaT_max_norm = (maxop1(T) - minop1(T)) / (maxop1(T) - T_ambient)

ΔT_max 的局限性在于它只反映两个极端点，无法描述中间区域的温度分布形态。因此，需要结合RMSE进行综合评价。

4.4 均方根误差（RMSE）：全局离散度度量

RMSE 是温度均匀性评价中最具统计意义的指标。其计算流程如下：

1. 定义参考温度 T_ref：通常取MZI波导区域的平均温度 T_avg，或目标调谐温度 T_target。
2. 采样点选取：在MZI两臂波导中心线上等间距取 N 个点（建议 N ≥ 50）。
3. 计算残差平方和：Σ(T_i − T_ref)²。
4. 开方取平均：RMSE = √(残差平方和 / N)。

# Python 后处理脚本：从COMSOL导出数据后计算RMSE
import numpy as np

def compute_rmse(temperatures, ref_temp=None):
    """
    计算温度均匀性RMSE
    
    参数:
        temperatures: ndarray, 采样点温度数组
        ref_temp: float, 参考温度，默认取平均值
    
    返回:
        rmse: float, 均方根误差
    """
    if ref_temp is None:
        ref_temp = np.mean(temperatures)
    residuals = temperatures - ref_temp
    rmse = np.sqrt(np.mean(residuals ** 2))
    return rmse

# 示例：假设从COMSOL导出了50个采样点的温度数据
sample_temps = np.array([348.2, 348.5, 348.1, 347.9, 348.3, ...])  # 实际数据
rmse_value = compute_rmse(sample_temps)
print(f"RMSE = {rmse_value:.3f} K")

4.5 三指标联合评价流程

下图展示了温度均匀性指标的计算与评价流程，从COMSOL仿真结果到最终方案决策的完整链路：