COMSOL模拟大功率激光系统中由于热变形导致的焦点移动
热透镜的原理
模拟激光聚焦系统
模拟吸收介质中的射线传播
-
几何光学 — 计算射线轨迹。 -
固体传热 — 计算透镜中的温度。 -
固体力学 — 模拟透镜的热膨胀。 -
动网格 — 使与透镜相邻的域中的有限元网格变形。
计算射线轨迹
式中, 是光线位置, 是波矢量, 是角频率。波矢量和角频率的关系为
(2)
式中, 是真空中的光速。在吸收介质中,折射率可以表示为 ,其中 和 是实数。
当光线进入和离开透镜时,根据斯涅尔定律,它们会发生折射,
(3)
式中, 和 分别是入射角和折射角。
计算每条射线传输的功率
每条光线传输的总功率 ,在非吸收域中保持不变。在均匀的吸收域中,功率呈指数衰减
(4)
式中, 是射线的自由空间波数。
为了将由光线的功率损失用作固体传热 接口的源项,有必要在吸收域中添加一个沉积射线功率节点。对于由于所选域中的射线衰减引起的体积热源,这个节点定义了一个变量 (国际单位:)。当光线通过透镜传播时,它们对 值有贡献,
(5)
式中,(SI 单位:W)是某条光线的传输的功率, 是光线的总数, 是狄拉克函数。实际上,每条射线不能在其精确位置产生热源项,因为射线在空间中占据无限小的点,而底层网格单元具有有限大小,因此每条射线损失的功率均匀分布在射线所在的网格单元上。
计算温度
可以通过求解传热方程来计算透镜的温度
(6)
式中, 是介质的热导率。使用热通量节点,在暴露于周围空气的所有边界上施加对流冷却
(7)
温度的变化对透镜中的总非弹性应变贡献了一个热应变项 。该热应变由下式定义
(8)
式中, 是热膨胀系数, 是介质的温度, 是参考温度。产生的位移场 由固体力学接口计算。
获得合理的解
如果光束传输的功率非常低,则光线损耗到周围环境的能量不会明显改变介质的温度。但是如外力和热源等其他现象,仍然有可能改变透镜的形状或温度。
在这个案例中,需要先计算域内的位移场和温度,然后再计算射线轨迹。这是一个单向或单向耦合,因为温度变化和结构变形会影响射线轨迹,但反过来却不是这样。
如果光束传输的功率足够大,吸收介质中的能量损耗可能会产生足够的热量,从而显著改变区域的形状或介质中的折射率。在这个示例中,光线轨迹会影响在周围域上定义的变量,例如温度,而这些变量又会影响光线轨迹,这是一个双向或双向耦合。
在这个示例中,我们假设激光器以恒定功率运行,因此最好使用稳态研究步骤计算温度和位移场。但是,射线轨迹是在时域中计算的。
为了在射线轨迹与温度场和位移场之间建立双向耦合,我们首先创建一个稳态研究步骤来模拟透镜的加热和变形,然后添加一个射线追踪研究步骤来计算射线轨迹。最后,使用 For 和 End For 节点将相应的求解器封闭在一个循环中。下图显示了在射线轨迹与温度场和位移场之间建立双向耦合的求解器序列。
For 和 End For 节点之间的节点将重复 For 节点的设置窗口中指定的次数。此外,每次运行求解器时,它都会使用前一个求解器的解。通过这种方式,可以在两个研究之间建立双向耦合,并在它们之间进行迭代,直到获得合理的解。
热透镜模拟结果
3000w 光束(左)和 1w 光束(右)沉积射线功率的比较。 为了进行比较和可视化,沉积功率的颜色表达式已被归一化并用对数尺度表示。
我们已经看到,大功率激光系统中的温度变化和由此产生的热膨胀会明显改变光束的焦点。您在设计此类系统时可以使用射线光学模块考虑热变形引起的焦移。
文章来源:COMSOL