激光烧蚀模型

激光辐照在材料上，一个瞬时的高能能量使其从固态转换为气态，主要模拟这个过程中温度场的变化以及材料表面烧蚀时其形貌的演变，不考虑材料的熔化过程。

构造一个固液转化的质量分数，中间量

激光热通量

对流热辐射

升华热

温度场是连续的，温度逐渐升高并达到材料的熔点,材料从固态 液态 转化为气态

材料表面的形变，移动网格（动网格技术（ALE，欧拉拉格朗日方法，有形变但质量守恒），变形几何：体积或者质量不守恒（适合熔覆，打孔技术））

模型选用的材料：铝合金

考虑到温度对材料热物理性能的影响，需要考虑材料固相和液相的密度（rho），比热(Cp)，热导率(k)，

flc2hs(T-Tm,deltaT)[1/K]  以Tm为中心，左右各deltaT，在两个断点的位置分别为0,1，[1/K]消去单位

计算有效导热率时需要考虑潜热的影响：表观热熔法

公式中的Cp:为有效热熔  ，Lm:熔化潜热，Lm的系数为一个高斯分布

模型建立过程中应该注意：一般上表面都为零，这是因为多个公式中都是按照零的初始位置建立的

不考虑辐射项

B中不需要单位归一，自动归一

二、对流，一般是设置对流换热系数，这个模型中还需要考虑蒸发带走的热量，蒸发的热量通过对流换热Hv来包含或者通过辐射来包含，经验参数调试（文献中Hv的值较大，蒸发产生的热量通过对流来实现），一般的文献中表面的温度大于熔点，

材料温度大于熔点使蒸发，定义关于Hv的斜坡函数，较大的斜率，Hv很快的变化，但熔点的变化非常小

二、变形几何

变形几何中：首先选择整个几何体为自由变形，让几何体可以动起来

烧蚀发生在上表面，两个侧壁基本不动，侧壁，指定网格，位移为0

法向网格速度：

蒸发带走的热量用对流热通量表示，实际上的蒸发热率小于沸点，或者高一点点，即热通量是正数

打孔式向下的，即就加一个负号表示

网格化：

收敛：有很多方法可以加速收敛

1. 求解器中设置   

   一般对于变形几何和动网格问题，选择全耦合

变形网格一般选择重划分网格

重划分网格时，会有几个对应的重划分网格的条件，一般选择网格质量，0.2表示当网格质量小于0.2时，重新划分网格，即重新开始计算，但需要注意这个时候重新计算时，已经跑了很多步了，步长已经很大了，再次重新划分网格时步长太大，很容易造成不收敛，所以这时需要重新设置初始步长0.001，这个步长设置应该和最开始设置时的初始步长在一个量级，打开一致初始化

上面的设置一般对于毫秒激光已经足够计算，但对于飞秒和皮秒激光而言，超导

还需要继续设置其瞬态求解器的时间步进

这里的初始步长为最开始计算时的步长，

求解器采用的步长为精确

最大BDF一般设置为2或3就已足够

对于超快而言，有时候还是不收敛，这是因为，飞秒激光的步长非常小，在非常小的时间内给材料表面一个非常大的能量，这时用时间的step函数，（找相关的超快案例）

注：comsol 的求解步长，第一个为初始步长1e-4，第二个步长为1e-4/2，第三个步长1e-4/2/2......

在1/1e-4秒，将一个值加速到其最预定值，若模拟过程为1秒，则其可以在极短的时间内将值加到预定值，对其结果几乎没有影响，若从0-0.5秒才加到1，则影响会非常大，前0.5秒没有加上能量，与实际不符

step函数应该配合初始步长来使用，一般step函数的过度区大小是初始步长的4到5倍，

注：comsol ，初始步长为1e-4，第二步长为2e-4,第三个步长3e-4......,即基本在两个步长时就已经加至预定值

step函数做平滑避免初始值不一致

step函数加到一些突变值中，激光能量突变

eta*Qs*step2(t[1/s]) 

注意：step2函数的书写，关于时间的函数，且要去掉量纲，因为在时间阶跃设置时没有单位，去掉单位

没有融化，直接汽化，蒸发热，斜坡函数，对流热通量来表示蒸发热，调节Hv的标准，表面温度大于熔点，注多个材料的Hv值是通用的

比热：利用表观热熔模型，考虑融化潜热的影响

网格设置渐变网格

对于微妙毫秒量级的，选用重划分网格，求解器，全耦合

对于皮秒纳秒飞秒的超快加工过程，需要设置瞬态求解器，

求解器步长为精确，设置初始步长，同时，结合时间的平滑函数，即时间阶跃函数，step函数，设置，并将step函数加在激光能量项中，表示激光能量的瞬时加载，时间非常短，能量会发生突变，很容易不收敛。

这也适用于表面有很大变形的问题，

飞秒激光烧蚀：双温方程，1. 哪些区域中晶格的温度大于熔点的，利用后处理直接将其屏蔽掉，就可以得到烧蚀的坑道形状，2. 利用变形几何

3. 利用标准的公式，计算蒸发量，法向网格技术，得到坑道的形状

复杂的可能会考虑表面联合效应和其他的效应,熔池效应等，飞秒过程非常复杂

修正：​编辑变形几何底下设置的​编辑法向速度没有负号