在 COMSOL Multiphysics 中模拟激光与材料的相互作用

作者 Walter Frei

2015年 6月 22日

用户经常会问我们:能不能在 COMSOL Multiphysics 中模拟激光与材料的相互作用及加热现象?具体的答案需要取决于您求解的问题类型,因为不同的问题须采用不同的建模技术。今天,我们将探讨应用于激光照射材料的加热问题的多种模拟方法。

简要介绍激光与材料相互作用的模拟

尽管激光光源的类型多种多样,但它们的输出却非常相似。激光接近于单频(单波长),并且具有相干性。通常情况下,输出的激光被聚焦为一束狭窄的准直光。这个相干的单频准直光源在很多领域被用作精准热源,其中包括癌症治疗焊接退火材料研究半导体加工等。

当激光照射到固体材料上时,它的部分能量会被材料吸收,导致局部加热。当然,激光也可以加热液体和气体(以及等离子体),但是流体加热总会引发显著的对流效应。在本篇博客文章中,我们忽略对流,只关注固体材料的加热现象。

对于特定的激光波长,固体材料对激光可能是半透明的,也可能是完全不透明的。因此根据材料透明度的不同,需要选择不同的激光热源模拟方法。此外,我们一定要关注材料与激光波长的相对尺寸。若激光为紧聚焦光束,其适用方法将与宽光束完全不同。若与光束相互作用的材料的几何特征与激光波长尺度接近,那么我们还需要考虑光束具体是怎样与这些小型结构发生相互作用的。

在开始模拟激光与材料的相互作用之前,首先应该确定所模拟材料的光学性质,包括材料在特定的激光波长下以及红外波段范围的特性。您还要了解加热物体的相对尺寸,以及激光波长与光束特性。这些信息将有助于引导您找到正确的建模方法。

表面热源

如果在具有特定波长的激光的照射下,材料为不透明或几乎不透明,那么我们便可以将激光看作表面热源。利用 COMSOL Multiphysics 5.1 版本中“传热模块”的沉积的电子束功率 特征(如下图所示),我们就可以轻易地将激光处理为表面热源。当然,就算仅使用 COMSOL Multiphysics 的核心软件包,我们也能方便地手动创建这个表面热负荷,如链接中的案例所示

表面热源这一说法本身假定了光束与吸收其能量的材料之间的距离是可忽略的(相对于被加热物体的尺寸而言)。有限元网格只要精细到足以解析温度场和激光光斑尺寸即可。我们没有对激光本身进行明确的建模,同时假设被材料反射出去的部分激光不会再次反射回来。使用表面热负荷时,您必须手动添加材料在该波长激光照射下的吸收率,并适当调整沉积电子束功率的范围。

图像展示了使用沉积电子束特征模拟激光与材料的相互作用。
使用“传热模块”中的“沉积电子束功率”特征模拟两束交叉的激光。绘图显示了生成的表面热源。

体热源

在材料为半透明的情况下,激光功率的沉积将发生在域内,而非表面上。根据几何尺寸和波长的相关性,任何不同的方法都可能是合适的。

射线光学

如果被加热物体的尺寸比波长大得多,但是激光在穿过一组光学元件时发生了会聚和发散,并且很可能被镜面反射,那么这种情况下,射线光学模块是您的最佳选择。在该方法中,激光被视为一种射线,软件会追踪它穿过均匀材料、非均匀材料以及有损材料时的路径。

当激光穿过有损材料(例如光学镜片)并照射到表面时,部分功率沉积会对材料产生加热作用。我们可以通过复折射率来模拟域内的吸收现象;对于材料表面,可以使用反射系数或吸收系数。这些属性都可以随温度变化。对于有兴趣实践该方法的用户,“案例下载”中的热致焦距漂移模型是一个很好的入门教程。

图像展示了穿过两个透镜聚焦的激光束。
激光束穿过两个光学透镜进行聚焦。由于激光强度较高,透镜升温并导致焦点偏移。

Beer-Lambert 定律

如果被加热物体和激光光斑尺寸都远大于波长,此时应当使用 Beer-Lambert 定律来模拟材料内激光的吸收现象。这种方法假设激光束为完全平行的单向光束。

当使用 Beer-Lambert 定律方法时,材料的吸收系数与材料表面的反射率必须是已知的。这两种材料属性都可能随温度而变化。至于建立此类模型的正确方法,您可以参阅我们之前发布的博客《借助 Beer-Lambert 定律模拟激光与材料相互作用》。

若您已经知悉入射激光的强度,并且材料内或边界处未发生激光反射,可以采用 Beer-Lambert 定律方法。

图像展示了使用 Beer-Lambert 定律模拟半透明固体在激光照射下的加热情况。
使用 Beer-Lambert 定律模拟半透明固体在激光照射下的加热情况。

波束包络方法

如果被加热的域较大,但激光却集中聚焦于域内,这时射线光学与 Beer-Lambert 定律建模法都无法准确地求解焦点附近的场与损耗,这是因为两种方法都将激光看作射线,因而无法直接求解 Maxwell 方程组。在这种情况下,“波动光学模块”中的波束包络方法是您的最佳选择。

当光场包络变化缓慢时,波束包络方法可以求解完整的 Maxwell 方程组。如果波矢在整个域中几乎是已知的,我们也清楚光的大致传播路线,这种方法就非常适合。模拟激光聚焦和波导结构(例如 Mach-Zehnder 调制器环形谐振器)的案例就符合这类情况。因为光束方向是已知的,因此传播方向上的有限元网格可以非常粗化,从而有效减少了计算量。

图像显示了在圆柱形材料域内聚焦的激光束。
在圆柱形材料域内聚焦的激光束。绘制了入射处与材料内的激光强度以及网格。

通过电磁热源 多物理场耦合,我们可以将波束包络方法与固体传热 接口相结合。在添加物理场 下点击添加激光加热 接口后,系统将自动建立上述耦合。

截图展示了在COMSOL Multiphysics 中的激光加热接口。
激光加热接口添加了 波束包络 固体传热接口,以及它们的多物理场耦合。

全波

最后,如果被加热结构的尺寸与波长相当,便有必要求解完整的 Maxwell 方程组,而无需假定激光在建模空间内的任何传播方向。这时,我们需要使用“波动光学模块”或 RF 模块中的电磁波,频域 接口。此外,“RF 模块”中还内置有微波加热 接口(类似于上文的激光加热 接口),能将电磁波,频域 接口与固体传热 接口耦合在一起。虽然这样命名,但实际上“RF 模块”和微波加热 接口仍适用于较宽的频带范围

当执行全波方法时,为了解析激光波长,有限元网格需要足够精细。由于光束可能在所有方向上散射,网格大小必须尽量地统一。模拟平面波照射下金纳米球中的损耗(如下图所示)是一个合理使用电磁波,频域 接口的绝佳案例。

图像展示了被激光加热的金纳米球。
激光加热金纳米球。图像沿网格绘制了球体内的损耗与周围的电场大小。

模拟材料内部与周围的传热、对流和再辐射

您可以从上述五种方法中选择任意一种,来模拟固体材料中由激光光源引起的功率沉积。不过若要模拟材料内部与周围的温升和热通量,还需要用到固体传热 接口。COMSOL Multiphysics 的核心软件包提供了这一接口,它能够模拟固体传热,并包含固定温度、绝缘和热通量边界条件。固体传热 接口还包含其他多种边界条件,可用于模拟被加热体与周围空气或流体的对流传热,以及趋于已知环境温度的辐射冷却过程。

有时候,您可能希望有这样一种流体:它可以对问题中的材料进行大幅加热或冷却,但不能被当做边界条件近似处理。为此,您需要使用传热模块CFD 模块对流体流动进行建模,因为两个模块都可以同时对温度和流场进行求解。此外,两个模块也都可用于解决层流和湍流流动问题。不过“CFD 模块”还包含了一些额外的湍流建模功能,详情请参阅之前发表的一篇博客

如果您认为在不同温度下,被加热物体与周围任何物体之间可能会发生显著的辐射,“传热模块”中的一个附加功能便可以计算灰体辐射角系数和辐射传热。快速热退火教程模型对这一功能进行了详细的演示。如果您希望温度变化较为显著,那么还应在仿真中加入随波长变化的表面发射率

如果您研究的材料对激光透明,那么它们很可能对热(红外波段)辐射是呈半透明的。红外光既不是相干光,也不是准直光,因此我们不能用上述任何一种方法来描述半透明介质中的再辐射。这时我们可以使用另一种方法,即参与介质中的辐射。当辐射会在材料内导致十分显著的热通量时,这种方法便适用于模拟其传热现象。“案例下载”中的玻璃板辐射冷却案例便使用了这种方法。

总结

在本篇博客文章中,我们探讨了多种可在 COMSOL Multiphysics 的软件环境中用于模拟激光加热固体材料现象的建模技术,其中包括表面加热和体加热等方法,同时我们还对传热建模功能进行了简述。目前为止,我们只考虑了无相变固体材料的加热现象。对于液体和气体这类会发生相变的材料,我们将会在后续博客中进行探讨。敬请关注!


评论 (27)

正在加载...
强 赵
强 赵
2018-09-15

我用proe建了尺寸是微米级的模型,导入comsol中,comsol总是报错“Internal error code: 900”。所以就想问问comsol可以就进行微米,甚至是纳米的仿真吗?如果可以如何解决这个问题?

爱春 樊
爱春 樊
2019-09-18

不可以直接在comsol的几何模型里面创建几何体么

王 刚
王 刚
2018-10-19

通常可以做这种微/纳米级别的建模的,您提到的这种错误有较大可能性是几何结构在导入时出现了错误,建议发给 Support 帮你检查一下。

旋 周
旋 周
2018-11-30

有人知道如何使用沉积的电子束功率来模拟激光加热吗?我找不到激光光斑大小输入的地方。

飞 赵
飞 赵
2019-01-03

请问沉积电子束功率 激光的波长在何处设置

欣 李
欣 李
2019-04-18

你好,请问在5.3中就没有沉积的电子束功率了吗

wei bao
wei bao
2019-04-25

激光的波长信息直接反应在功率参数中了,不需要单独设置

泽昊 李
泽昊 李
2019-07-10

HUA ZJ
HUA ZJ
2019-07-15

你好,请问热源是纳秒级的激光还能否模拟到它的传热过程

振江 陈
振江 陈
2019-08-05

在激光热疗中会涉及一些吸收光谱的问题(输入特定波长的激光,能量的利用率比较高。对周围的细胞的损伤相对来说比较小一点),那么吸收光谱(例如金纳米颗粒)该如何模拟呢?COMSOL能解决吗?
注意(吸收光谱并非散射光谱,之前我问过COMASOL工作人员,他告诉我贵公司案例库有,但是并不能解决我的问题)

changjaing sun
changjaing sun
2020-12-03

可以用comsol模拟螺旋波产生的等离子体与电磁波的相互作用(加热等离子体)问题吗

孝通 苑
孝通 苑
2021-09-22

能模拟激光穿过多层材料反射会来的材料光谱信息吗

Qingbin Yuan
Qingbin Yuan
2021-09-28 COMSOL 员工

可以的,您可以看一下这个链接中对分布式布拉格反射镜的模拟:http://cn.comsol.com/model/distributed-bragg-reflector-19275
其中DBR是一种由交替电介质层构成的周期性结构,可以通过在射线光学中使用薄层来等效处理。

淼 刘
淼 刘
2021-10-20

激光加热融化金属这种相变固体的博客出了吗?

洋洋 张
洋洋 张
2021-10-21 COMSOL 员工

固体熔化相变界面计算可以参考案例:http://cn.comsol.com/model/tin-melting-front-6234

婷婷 钟
婷婷 钟
2021-11-22

可以研究激光器内的激光晶体受激光或者说是泵浦光的作用吗,如果可以,如果想要模拟激光晶体的损伤过程,怎么设置

Liwen Yang
Liwen Yang
2021-11-26 COMSOL 员工

如果激光晶体的损伤是指结构损伤,那么可以在结构力学模块下考虑损伤建模,可以参考案例 http://cn.comsol.com/model/brittle-damage-in-uniaxial-tension-62761

君贤 张
君贤 张
2022-06-29

激光加热金纳米球 模型的温度该如何计算呢?

Qihang Lin
Qihang Lin
2022-07-04 COMSOL 员工

该类问题属于电磁热分析,建议您学习微波加热案例如何仿真:http://cn.comsol.com/model/microwave-oven-1424

含 韦
含 韦
2023-04-03

激光加热液体这种会发生相变的材料的博客出了吗?

越 赵
越 赵
2023-04-04 COMSOL 员工

您好,目前还没有激光加热液体相变的博客。

啊喵喵
啊喵喵
2023-04-13

请问如何模拟激光辐照材料,导致材料损伤啊

Liwen Yang
Liwen Yang
2023-04-19 COMSOL 员工

一般激光的入射热通量可以作为结构表面分布的热源来处理,从而通过固体传热计算温度分布,再耦合固体力学计算热应力。在力学分析中考虑材料损伤,可以引入“损伤”本构来计算。关于激光加热的温度计算,您可以参考案例 https://cn.comsol.com/model/laser-heating-of-a-silicon-wafer-13835 ;关于力学损伤建模,您可以参考案例 https://cn.comsol.com/model/brittle-damage-in-uniaxial-tension-62761

Yang Liu
Yang Liu
2023-04-19

矩形激光 矩形光斑怎么设置

Qihang Lin
Qihang Lin
2023-04-25 COMSOL 员工

需要您自己写入矩形的表达式,可以考虑通过组合两个与位置相关的step函数实现。

Yang Liu
Yang Liu
2023-05-15

矩形激光的扫描路径的公式需要用到分段函数来设置吗

Yuhong Zhang
Yuhong Zhang
2023-05-18 COMSOL 员工

您好,comsol里面有很多的函数类型,分段函数、波形函数、三角函数等,您可以根据您想要实现的具体路径形式进行相应选择,关于激光扫描的路径设置,您可以参考案例:https://cn.comsol.com/model/laser-heating-of-a-silicon-wafer-13835.

浏览 COMSOL 博客