LD端面泵浦激光晶体热效应分析.doc

LD端面泵浦激光晶体热效应分析李健1 孙尧1 李涛1 张帅一1 于果蕾1 1. 山东师范大学物理与电子科学学院，山东 济南 250014摘要：通过求解泊松热传导方程，得出了矩形截面Nd:YVO4 ,Nd:GdVO4和Nd:GdYVO4激光晶体中心泵浦和偏心泵浦情况下，晶体泵浦端面的温度分布。研究表明，偏心泵浦时，晶体泵浦端面最高温升较中心泵浦分别下下降了11.4%,13.5%,7.7%。并且计算了Nd:YVO4,Nd:GdVO4和Nd:GdYVO4激光晶体端面在达到热平衡时的热形变量，三种晶体的最大热形变量分别为2.6045，0.536，0.8549。关键词：Possion方程，激光晶体，偏心泵浦，热形变一、 引言激光二级管（LD）泵浦的全固态激光器具有体积小、效率高、稳定性好和寿命长等优点，再国防、光电子产业、光通讯和医疗卫生等领域有着重要的应用。激光晶体是全固态激光器中最重要的核心部分，在很大程度上它决定了激光器的输出特性，为了获得高效激光输出，在一定运转方式下选择合适的激光晶体是非常重要的。在全固态激光器中，Nd:YVO4和Nd:GdVO4以其优良的物理、化学及激光特性，被广泛的应用在各类激光产品中12。最近在Nd:YVO4和Nd:GdVO4的基础上，一种新的钒酸盐晶体Nd:GdxY1xVO4逐渐引起了人们的注意。Nd:GdxY1-xVO4是由Y(Gd)部分地代替Nd:GdVO4(Nd:YVO4)中的Gd(Y)形成的新的钒酸盐晶体。三种晶体具有相同的晶体结构，它们的吸收波长中心都在808nm附近。Nd:YVO4晶体具有较高的吸收系数和受激发射截面，Nd:GdVO4晶体则具有较高的热传导率3。 Nd:GdYVO4具有和Nd:YVO4，Nd:GdVO4相同的晶体结构，但由于融合了两种晶体的特点，具有自身的优势4。在大功率抽运的过程中，非辐射跃迁产生的热加上当激光介质吸收泵浦光能量产生激光辐射的同时，有相当一部分泵浦光能量会转变成热能消耗于激光晶体内部56。晶体的外部采用半导体制冷或循环水冷却的方式使热量得以扩散，其结果是在激光晶体内形成非均匀的温度梯度分布，这种分布对振荡的激光作用的效果可等效为一凸透镜7。激光晶体的热透镜效应直接影响激光器的性能。本文通过对LD端面泵浦，周边恒温冷却的矩形截面Nd:YVO4 ，Nd:GdVO4和Nd:GdYVO4晶体进行了分析，利用matlab求解泊松热传导方程，得出了Nd:YVO4，Nd:GdVO4和Nd:GdYVO4激光晶体内部中心泵浦以及偏心泵浦温度场分布，并且计算了3种晶体在达到热平衡时的热形变量。二、 矩形截面激光晶体温度分布的数学模型和计算方法1 LD端面泵浦，周边恒温的矩形截面激光晶体热模型Nd:YVO4（ Nd:GdVO4，Nd:GdYVO4）晶体吸收部分泵浦光能量产生热量，通过紫铜块中循环水向周围散失。当泵浦光通过激光晶体中心时，热模型如图1所示8紫铜块泵浦光光斑晶体体xy图1、激光晶体热模型简图Fig1.schematic diagram of laser crystal thermal model（1）激光晶体的边界条件由于把晶体放置在紫铜块中用循环水进行冷却，所以假定晶体侧面温度保持相对稳定，设为，为了数学处理上方便，可以先预设为0，最后再叠加上冷却水的温度。晶体两个通光端面与空气接触，热交换的热量远远小于通过晶体侧面热交换流出的热量，因此可以假设晶体两个端面绝热，边界条件为 （1） （2） （3） a,b为晶体矩形截面的长和宽，l为晶体长度。（2）晶体内部热传导方程热传导方程的一般形式为91011 (1)其中，u为激光晶体的温度；为晶体的热导系数；为泵浦产生的内热源的发热率；为晶体的密度；c为晶体的比热。假设激光晶体的物理性质不随u而发生变化，在晶体侧面传导冷却的情况下，端面是空气对流冷却，其边界热交换系数远远小于传导冷却的边界热交换系数，因此端面散发的热量远远小于侧面热传导散发的热量，故可以忽略晶体纵向散热，那么方程就可以写成二维热传导方程。 （2）在圆对称高斯光束泵浦时，内热源的发热率为： （3）其中，Q为激光晶体中的总热耗；为晶体的热吸收系数；为泵浦光斑的大小；l为晶体的长度;r为晶体的径向坐标；为在z0时的泵浦光入射的热流密度。稳态时，热传导方程可写成 (4) 假设 (5)可得到 (6)故热传导方程表示为求解方程 （7）2Nd:YVO4，Nd:GdVO4和Nd:GdYVO4晶体内部温度场分布表1、Nd:YVO4 ,Nd:GdVO4和Nd:GdYVO4晶体激光特性参数12book1. laser parameters of Nd:YVO4 ,Nd:GdVO4 and Nd:GdYVO4激光晶体Nd:YVO4Nd:GdVO4Nd:Gd0.5Y0.5VO4吸收系数32.719.94晶体尺寸3310mm33310mm33310mm3掺杂浓度0.50.5%0.5热导系数（）0.050.1170.125热膨胀系数（）11.377.39.33在泵浦功率为12W，泵浦光斑半径时，通过求解Possion方程，得出Nd:YVO4 ,Nd:GdVO4和Nd:GdYVO4晶体内部温度场分布，如图2,3,4所示。3种晶体的最高温升分别为312，285，140，由此看见，Nd:GdYVO4的热性能明显优于Nd:YVO4 和Nd:GdVO4。 图2 Nd:YVO4晶体泵浦端面温度分布图 图3. Nd:GdYVO4晶体泵浦端面温度分布Fig2. temperature field distribution diagram of Nd:YVO4 Fig3 temperature field distribution diagram Nd:GdYVO4图4. Nd:GdVO4晶体泵浦端面温度分布Fig4 temperature field distribution diagram Nd:GdVO4图5泵浦光中线（x0,y=0）上温度分布Fig5.temperature distribution diagram of crystal centre从图5中可以得出： Nd:YVO4 ，Nd:GdVO4和Nd:GdYVO4晶体前部中心温升比后部高很多。由于Nd:GdYVO4晶体具有相对其他两种晶体而言最小的晶体中心最高温升以及最大的热导系数，所以温度沿z轴变化最小。Nd:YVO4 和Nd:GdVO4晶体内部的最高温升差不多，但由于Nd:GdVO4热传导系数比较大，所以沿z轴温度衰减的快。3.偏心泵浦对Nd:GdYVO4和Nd:YVO4晶体内部温度场分布的影响在泵浦条件不变的情况下，调节耦合镜的位置，偏心泵浦Nd:YVO4,Nd:GdVO4和Nd:GdYVO4（x0.5mm；y0.5mm）。Nd:YVO4,Nd:GdVO4和Nd:GdYVO4晶体端面温度分布如图6所示。中心温度分别为280,251和130。图6.偏心泵浦激光晶体端面温度分布 从图6中可以看出，晶体内部温度由于泵浦光的偏心泵浦而下降，Nd:YVO4 ，Nd:GdVO4和Nd:GdYVO4晶体中心温度分别下降32，34，10。当偏心泵浦晶体时，晶体内热源距离冷却周边比较近，所以晶体中心温升较低。偏心泵浦时，晶体泵浦端面最高温升较中心泵浦分别下下降了11.4%,13.5%,7.7%。4Nd:YVO4 ，Nd:GdVO4和Nd:GdYVO4晶体端面热形变当激光晶体内部形成稳定温度场之后，由于晶体受热膨胀会使晶体泵浦端面发生热形变，这种热形变是使激光晶体产生热透镜效应的主要原因之一。设激光晶体内部一点温度为0，当达到热平衡之后，温度变为。；其中为激光晶体热膨胀系数。图7给出了当在泵浦功率为12W，泵浦光斑半径时，三种激光晶体沿泵浦光中心方向（x=a/2;y=b/2）上泵浦端面具有的热形变分布对比图。从图中可以看出，在相同的泵浦条件下，Nd:YVO4晶体有着最大的热形变量，Nd:GdYVO4的热形变量介于Nd:YVO4,Nd:GdVO4之间。三种晶体最大的热形变量分别为2.6045，0.536，0.8549。热形变改变了晶体通过方向上的有效长度，并使晶体产生伴随热效应的高阶球差，在激光器谐振腔设计、调试时应予以考虑。图7.晶体泵浦端面热形变对比图Fig7 Thermal distortion contrast of the pumped faces 三、 结论本文通过对LD端面泵浦激光晶体工作特点的研究，提出了矩形截面激光晶体热模型，利用激光晶体周边恒温两端绝热的边界条件求解泊松热传导方程，得出了泵浦功率为12W，泵浦光斑半径时，Nd:YVO4,Nd:GdVO4以及Nd:GdYVO4激光晶体在中心泵浦和偏心泵浦两种情况下的内部温度场分布。并且计算了在达到热平衡时，三种晶体泵浦端面的热形变，对于改善激光晶体热效应影响，提高激光器的稳定性具有指导作用。参考文献1 赵宏明，赵圣之，杨克建等，LD泵浦Nd:YVO4/KTP内腔倍频声光调Q理论和实验研究J，光子学报，2004，33（11）：129412972 赵致民，李隆，田锋等，高功率端面泵浦腔内倍频瓦级绿光激光器J，激光技术，2003，27（04）：3313333 Jie Liu, Qianqian Peng, Jimin Yang, Qiuxia Jiang, Jingliang He, Lianjie Qin, Xianlin Meng,Diode-pumped Nd:YxGd1-xVO4 crystal continuous-wave laserJ,CHINESE OPTICS LETTERS,2004, 2(1): 29304 卓壮，姜其畅，苏艳丽等，新型钒酸盐混晶Nd:Gd0.42Y0.58VO4/Cr4+:YAG被动调Q锁模特性研究J，中国激光，2005，32（10）：142914325 Fan T Y,Heat generation in Nd:YAG and Nb:YAGJ,IEEE J Quantum Electron,1993,29(6):145714596 张行愚，赵圣之，王青圃等，中国激光，2000，27（9）：7777817 Innocenzi M E,Yura H T,Frucher C L,etc,Thermal modeling of continuous-wave end-pumped solid state lasersJ,Appl Phys Lett,1990,56(19):8318338 李隆，史彭，刘小芳，甘安生，陈浩伟，白晋涛，LD端面抽运长方形激光晶体的热效应J,光电工程，2005，32（4）：35389 史彭，李隆，刘小芳，甘安生，华中文，偏心度对矩形截面Nd:GdVO4晶体热效应的影响J，电子激光，2005，16（10）：1187119210 M.E.Innocenzi,H.T.Yura,C.L.Fincher,andR.A.Fieids,Thermal modeling of continuous-wave end-pumped solid-state lasersJ,Appl.Phys.Lett,1990,56(19):1831183311 Ananada K.Cousins,Temperature and Thermal Stress Scaling In Finite-Length End-Pumped Laser RodsJ,IEEE JOURNAL OF QUANTUM ELECTRONICS,1992,28(4):1057106912 祝俐，张怀金，孟宪林等，Nd:GdVO4晶体生长及其1064nm的激光特性J，人工晶体学报，1999，28（03）：229232Analyse of LD End-pumped Laser Crystal on Thermal Effect Li Jian1 Sun Yao1 Li Tao1 Zhang Shuaiyi1 Yu Guolei11. College of Physics and Electronics, Shan Dong Normal University , 250014, Jinan, ChinaAbstract:Based on Possion equation,the c