COMSOL加热保温过程

在激光加热过程中往往需要考虑保温过程。举个简单的例子，某次加热过程，在高斯热源的作用下，前90s的时间热源一直在作用，现在保温30秒后我们需要查看他的温度，COMSOL是否可以实现这一功能，答案是可以的。

首先我们先来研究没有保温过程的加热，元件半径为15mm厚度为5mm材料为熔石英，材料参数由文献取得这里不过多赘述。取高斯激光功率为15 kW,光斑直径为2 mm，考虑前后表面吸收且吸收率相等的情况（增透光学元件），具体仿真参数由图1所示。

图1 仿真参数

给定相关对流等边界条件后，设置好求解器进行仿真分析。在不考虑体吸收的作用下，膜层吸收是影响温度分部的主要过程，温度由表面向内部扩散，造成局部温升的过程。

图2 表面温度分布图

图3 截面温度分布

如图2、图3为一组计算的温升分布三维图。从图也可以看出，因为考虑的是前后表面的吸收，所以样品前后表面的温度相对高于样品中间的温度。这也与熔石英样品的热导率较小有关，样品中温度的高低主要取决于热源带来的能量多少，而由于热传导导致的能量再分布是次要因素。

图4 高斯激光辐照熔石英样品某次计算的样品中心温升随时间的变化关系

图4为高斯激光辐照熔石英样品某次计算的样品中心温升随时间的变化关系可以看到在高斯激光辐照熔石英样品时，在前面提到的温度边界条件下，温度很容易就能达到稳态。在样品热导率等参数确定的前提下，达到稳态的时间取决于温度边界条件中的换热系数大小。换热系数越大，温度达到稳态所用的时间越短，反之，换热系数越小，温度达到稳态的时间越长。

90s的温度有了，我们如何实现保温过程呢？很简单，所谓保温也就是自然散热，也就是把热源去掉后的过程。我们都知道COMSOL软件也是MATLAB软件的分支，所以我们再添加热源上增加一个if函数来实现这一功能如if(t<=90,3*Q1*c/(pi*R1^2)*exp(-3*(x^2+y^2)/R1^2),0)，首先if函数是被COMSOL软件识别的，所以我们就不需要调用分段函数来实现，通过一个简单的操作来实现这一功能。

图5 保温30s后表面的温度变化

如图5所示，没有热源的作用后，表面温度降得很快，在120秒的时候最大温度只有56.6K左右，温度分布也降不在集中在光斑内部，而是随着时间向周围扩散。

图6 高斯激光辐照熔石英样品保温后样品中心温升随时间的变化关系

如图6为保温30秒，表面温度大概只有50K左右，90秒停止加热的一瞬间温度下降的最快，和加热温升过程类似，并且温度的下降也不是没有边界的，最后趋于一种平衡状态，或许受外界室温的影响。