📚 课程目录第3章光电协同仿真中的多物理场基础

3. 光电协同仿真中的多物理场基础

CPO（Co-Packaged Optics，共封装光学）把光引擎与交换 ASIC 放在同一封装基板上，光、电、热、力四个物理场在毫米甚至微米尺度内同时作用。任何一次“光电协同仿真”，本质上都是在求解这四个场之间的双向或单向耦合问题。理解每个场用什么方程描述、场与场之间通过什么物理量传递影响，是决定仿真能否收敛、结果是否可信的前提。

四个物理场的基本方程

光学场在 CPO 中主要涉及波导、微环、光栅耦合器和调制器。其基本描述是麦克斯韦方程组，工程上通常化为频域亥姆霍兹方程 ∇×(1/μᵣ ∇×E) − k₀²εᵣE = 0，用 FDTD 或 FDE/FEM 求解。对硅波导，关键输出量是有效折射率 n_eff、群折射率 n_g、传播损耗 α 和模式重叠积分。温度与应力会通过热光系数和弹光效应改变 εᵣ，这正是光学场与热、力场耦合的入口。

电学场包含两部分：一是器件级载流子输运，用漂移-扩散方程描述，包括泊松方程 ∇·(ε∇φ) = −q(p−n+N_D⁺−N_A⁻)、电子与空穴连续性方程；二是封装级互连，用传输线或 RLC 网络描述，关注插入损耗、回波损耗和阻抗匹配。对硅基调制器，载流子浓度变化同时改变折射率（等离子体色散效应）和吸收，因此电学解直接驱动光学参数。

热学场由傅里叶导热方程 ρc_p ∂T/∂t = ∇·(k∇T) + Q 描述，稳态下退化为 ∇·(k∇T) = −Q。CPO 的热源包括 ASIC 功耗（数十至数百瓦）、激光器、调制器和金属互连焦耳热。关键输出是温度分布 T(x,y,z) 与热梯度 ∇T。硅的热光系数约 1.8×10⁻⁴ /K，意味着 1 K 温升会让微环谐振峰漂移约 0.06–0.1 nm，这是热-光耦合最直接的体现。

力学场用弹性力学平衡方程 ∇·σ + f = 0 与几何方程、本构方程 σ = C:ε 描述。CPO 中力学问题的来源是热膨胀失配：硅（CTE≈2.6×10⁻⁶/K）、二氧化硅（≈0.5×10⁻⁶/K）、铜（≈17×10⁻⁶/K）与封装基板（≈10–17×10⁻⁶/K）在回流焊降温或工作温升时产生应力。应力通过弹光效应改变折射率，也可能导致光纤阵列与光栅耦合器错位，造成耦合损耗上升。

耦合机制：谁驱动谁

四个场之间的耦合并非对称，实际仿真中要区分“强耦合”与“弱耦合”。电-热耦合是最强的一环：调制器与 ASIC 的功耗直接作为热源 Q 注入热学方程，而温度又反过来改变载流子迁移率和器件电阻，形成双向迭代。在 CPO 中，ASIC 功耗可达 50–300 W，热流密度远超传统可插拔光模块，因此电-热耦合通常必须双向求解。

热-光耦合是 CPO 最敏感的单向链路。硅的热光系数 dn/dT ≈ 1.8×10⁻⁴ /K，微环谐振波长漂移约 0.06–0.1 nm/K。若微环 Q 值约 10⁴、3 dB 带宽约 0.1 nm，则 1 K 温升就足以让工作点偏离谐振峰，导致调制效率骤降。因此热-光耦合通常按“热场解 → 光学参数更新”的单向方式处理，但若涉及热调谐环，则需加入加热器功率与温度的闭环。

热-力耦合来自 CTE 失配。以硅光芯片键合到有机基板为例，假设回流焊峰值温度 260 ℃、室温 25 ℃，温差 ΔT = 235 K，硅与基板 CTE 差约 10×10⁻⁶/K，则等效热应变约 2.35×10⁻³，对应热应力可达百 MPa 量级。这个应力通过弹光效应改变硅的折射率（典型弹光系数约 −1.5×10⁻⁴ /MPa 量级），同时使光纤阵列相对光栅耦合器产生亚微米级位移，直接抬高耦合损耗。

力-光耦合在 CPO 中常被低估。波导在应力下产生双折射，TE/TM 模式有效折射率差变化，对偏振敏感器件（如微环、定向耦合器）影响显著。仿真中通常把力学场算出的应力张量映射到光学网格，再通过弹光张量修正 εᵣ。

仿真优先级与求解顺序

面对四个场，不可能也不必要每次都做全耦合。合理的优先级判断依据是“时间尺度”和“影响量级”。热学时间常数在毫秒到秒量级，电学在皮秒到纳秒，光学在飞秒量级，力学在毫秒到秒量级。因此工程上常采用分层求解：

先算热：用有限元求解稳态温度场，得到 T(x,y,z)。这是 CPO 中最慢但影响最广的场，也是其他场的输入。
再算电：在温度场基础上求解载流子分布与 I-V 特性，得到 Δn、Δα 和功耗 Q，必要时把 Q 回代热场迭代 2–3 次。
然后算光：把 Δn、Δα 和应力修正后的 εᵣ 代入光学求解器，计算 n_eff、损耗、谐振漂移和调制效率。
最后算力：若关注耦合对准或长期可靠性，再求解热应力与位移，反馈到光学耦合效率。

优先级不是固定不变的。若研究的是调制器带宽，电-光耦合应优先；若研究的是微环热稳定性，热-光耦合应优先；若研究的是封装翘曲导致耦合损耗，热-力-光链路应优先。判断原则是：先确定目标输出量，再倒推哪条耦合路径对它的灵敏度最高。

一个可执行的耦合仿真参数示例

以下给出一个用于教学演示的微环调制器多物理场耦合参数设置示例（数值为合理假设，非实测结果），可用于搭建协同仿真流程：

# 微环调制器多物理场耦合仿真参数示例（教学假设值）
光学:  λ0 = 1310 nm, 环半径 R = 8 μm, Q ≈ 12000
       波导损耗 α = 2 dB/cm, 群折射率 n_g = 4.2
电学:  掺杂浓度 N = 5e17 cm⁻³, 偏压 V = -1.0 V
       调制效率 VπL ≈ 1.2 V·cm, 结电容 C ≈ 80 fF
热学:  硅 k = 130 W/(m·K), SiO₂ k = 1.4 W/(m·K)
       热光系数 dn/dT = 1.8e-4 /K
       加热器功率 P_heat = 5 mW, 热阻 R_th ≈ 8 K/mW
力学:  硅 CTE = 2.6e-6 /K, 基板 CTE = 15e-6 /K
       回流温差 ΔT = 235 K, 弹光系数 ≈ -1.5e-4 /MPa
耦合:  热-光单向, 电-热双向迭代 3 次, 力-光单向映射

按上述参数，若加热器功率从 0 增至 5 mW，温升约 40 K，谐振波长漂移约 2.4–4.0 nm，足以覆盖整个 C 波段的一个信道间隔。仿真时应先确认热场收敛，再检查光学谐振峰是否按预期漂移；若漂移量明显偏离 dn/dT 估算值，通常说明网格或边界条件设置有问题，而不是物理模型本身出错。

多物理场耦合仿真最常见的失败原因不是方程写错，而是网格不匹配和时间尺度混用。光学网格通常在纳米级，热学网格在微米级，力学网格在十微米级。映射时若不做守恒插值，热源或应力会失真。建议在耦合界面处统一采用能量守恒的投影算法，并记录每次迭代的残差，残差不再下降时停止迭代，而不是盲目增加迭代次数。

总结来看，CPO 多物理场仿真的核心不是把四个场都算一遍，而是识别当前设计目标下哪条耦合路径主导，选择正确的求解顺序和耦合强度，再用合理的示例参数验证结果是否符合物理直觉。热场通常是起点，电-热双向是必须，热-光单向是高频路径，力-光耦合则按可靠性需求决定是否纳入。