基于MATLAB的Nd:CNGG激光器的熱效應研究

李奇楠，丁狀，向萬貴，羅小潔

基于MATLAB的Nd:CNGG激光器的熱效應研究

李奇楠，丁狀，向萬貴，羅小潔

（齊齊哈爾大學 理學院，黑龍江 齊齊哈爾 161006）

基于MATLAB的有限差分法對激光二極管（LD）端面泵浦Nd:CNGG激光器的熱透鏡效應進行了仿真分析。在LD端面抽運的條件下，設定泵浦光波長808nm和885nm，分別改變了泵浦光功率、晶體的尺寸和泵浦光束腰半徑的位置，得到了不同條件下的晶體內部溫度分布圖。分析得到，以上因素都會對晶體內部溫度分布造成影響。對采用808nm泵浦光，束腰位置不同時晶體的溫差，最高約68K，而使用平均溫升較低的885nm泵浦，溫差約10K。

激光器；Nd:CNGG；熱效應；MATLAB仿真

自上世紀60年代激光器誕生以來，眾多科研人員紛紛投入研究，為激光器及其相關技術的發展注入了極大的活力。此后，固體激光器的變化可謂日新月異，其中激光晶體是固體激光器不可或缺的一環。激光晶體種類繁多，但以摻雜Nd3+的一系列晶體更受科研以及相關人員的喜愛[1]。Ca3(Nb,Ga)2-Ga3O12(CNGG)是一種具有無序結構的晶體，它的吸收帶較寬，使其有利于進行激光二極管抽運[2]。因此摻雜Nd3+的Nd:CNGG晶體應運而生。其發射光譜中的935nm波長是水蒸氣的吸收峰，可以用作探測光源，將之倍頻也可以用作藍光光源。其熔點低，還可在空氣中生長，所以是一種有前景的激光材料。

在固體激光器中，最影響激光器穩定性的是激光晶體工作時的熱效應。它會影響激光器的輸出波長、輸出功率、模式穩定等，進而導致激光器的輸出光束質量變差，甚至不再輸出光。造成熱效應的廢熱來源主要有量子虧損、能級之間的能量差轉化、晶體的再吸收等[3-5]。由于激光晶體內部受熱膨脹不均，導致激光晶體內部密度不均，產生類似光線通過透鏡一樣的折射效果，即熱透鏡效應[6-7]。由于熱效應對固體激光器的種種負面影響，研究人員對晶體熱效應做了大量研究[8-11]。為了研究計算晶體的熱效應，向來都是研究端面泵浦或者雙端泵浦、選擇波長不同的泵浦、泵浦光分布對球差的影響等方面進行研究，少有對激光晶體本身相關的因素進行研究。本文通過建立激光晶體的熱模型，改變晶體的尺寸、泵浦光在晶體中的束腰位置和泵浦光的輸入功率，利用MATLAB進行理論模擬，尋找相對影響晶體內部溫度分布的因素。

1 理論分析

以簡單的端面泵浦平行平面腔為例，如圖1所示。

圖1 激光器實驗裝置示意圖

本文采用波長為808nm和885nm的LD作為泵浦源。泵浦光束從光纖通過耦合系統經過腔鏡M1照射在激光晶體Nd:CNGG的端面中心。一般采用銦箔包裹晶體，再將其置入紫銅夾中，紫銅夾注冷卻水循環散熱降溫，控制晶體周邊溫度恒定。散熱裝置示意圖如圖2所示。相較于808nm使用885nm的LD泵浦，激光晶體擁有更低的溫差，熱焦距更大[12]。

根據圖1, 2，建立沿晶體端面幾何中心抽運、邊緣恒溫的長方體晶體的熱模型，由此分析端面抽運的激光晶體內溫場分布狀態[13-15]。

圖2 熱模型示意圖

假定耦合到激光晶體端面的抽運光具有理想的高斯分布[16-17]，并且抽運光沿晶體中心平行于軸方向傳播。在端面（0）處的光強0的分布為[18]

式(1)中，ω為高斯光束的光斑半徑；為泵浦光束的功率。

泵浦光束在晶體中傳播時，光強會因晶體的吸收而逐步減弱，衰減程度呈e指數形式。那么在處的光強(,,)有：

式(2)中，為晶體的吸收系數。對于Nd:CNGG，對波長為808nm和885nm的光吸收系數分別為264m-1和50m-1。

在泵浦光工作時，晶體內部的光束可作為熱源，其溫度分布應遵循泊松方程：

其中，(,,)是晶體內部的熱功率密度；為各向異性的熱導率（熱傳導系數）。

晶體內部升溫的主要因素是晶體的熒光量子效應和內損耗吸收泵浦光，在只考慮二者的情況下，在處晶體的熱功率密度，可表示為

其中，為上述兩個主要因素決定的熱轉換系數。在Nd3+低摻雜濃度的晶體中，可以表示為

其中，λ為泵浦光束的波長808nm或885nm，λ為震蕩光束的波長935nm。在外加紫銅水冷的情況下，采用LD端面泵浦，晶體的側面熱擴散遠遠大于軸向，因此可以認為晶體的端面絕熱，熱量為徑向傳導[19]。設定水冷溫度為u，晶體側面溫度應與其相同，那么邊界條件可認定為

在LD端面泵浦的全固態激光器中，激光晶體的有效熱透鏡焦距可以表達為[20]

式(7)中，熱光系數dd=9.2×10-6/K[21]，0為高斯光束的束腰半徑。

2 仿真結果

上文已給出泵浦光波長、震蕩光波長和晶體的吸收系數，其余參數設定如下：Nd:CNGG晶體的熱導率=3.43W/m/K[16]，設定熱沉溫度u=278K，泵浦光功率為20W。

根據激光原理知識，高斯光束的光斑半徑與束腰半徑之間的關系有：

其中，0為泵浦光束束腰半徑的位置。設定高斯光束束腰半徑0=0.2mm。使用MATLAB，利用有限差分法即可計算出晶體內部的溫度分布。

當晶體長度=12mm，808nm泵浦光束腰位置0=2mm時，晶體內部的溫度分布如圖3所示。

圖3 L=12mm，z0=2mm，808nm晶體內部溫度分布

其它變量保持不變，改變泵浦光束的束腰位置。泵浦光波長808nm和885nm的晶體中軸線上溫度分布如圖4所示。不同的束腰位置，晶體內部溫度不同。束腰位置對應著溫度的拐點，泵浦光波長為885nm時，束腰位置對應的溫度為晶體內部的最高溫度。束腰位置越遠，晶體內部溫度越低。在采用波長808nm的光泵浦時，晶體內溫度最高與最低之間的溫差最大約68K。在泵浦光波長885nm，其余條件保持與808nm一致，晶體內部平均溫度低于808nm泵浦時的溫度，這與文獻[12]中的觀點一致，波長885nm泵浦Nd:CNGG晶體的熱效應更好。此外，還仿真了晶長=12mm，束腰位置0=6mm時，不同功率808nm和885nm束腰處方向的溫度分布，如圖5所示。

圖4 不同束腰位置的晶體內部溫度分布

由圖5可以看出兩種波長下，隨著功率的增加，晶體內部溫度上升，溫度分布逐漸不均勻。

考慮溫度與晶體長度的影響，設定束腰位置在0=2mm處，其它變量保持不變，泵浦光波長為808nm時，晶體內部溫度分布如圖6所示?？梢钥闯鼍w長度越長，晶體內溫度越低。但通過增加晶體長度實現降溫的效果不顯著，晶體長度為8mm時與晶體長度為16mm時的溫差僅為10K左右。在晶長為8mm時，根據式(7)，可得晶體的有效熱焦距與泵浦功率之間的關系，如圖7所示，可以看出熱焦距隨著泵浦功率的增大而減小，在23W后808nm的功率曲線趨于穩定。

圖5 不同功率下的晶體溫度分布

圖6 不同晶長的晶體內溫度分布

圖7 熱焦距隨泵浦功率變化

3 結論

利用MATLAB模擬計算得出，LD端面抽運，泵浦光的輸入功率、改變晶體的尺寸與泵浦光束腰半徑的位置都會對晶體內部溫度造成不同的影響。具體表現為：晶體長度越長，晶體內部溫度越低；泵浦光束腰位置離端面越遠，晶體內部溫度越低；泵浦功率越高，晶體內部溫度越高。對采用808nm泵浦光，束腰位置不同的溫差，最高約68K，而對平均溫升較低的885nm泵浦，溫差約10K。本文對激光實驗的晶體參數選擇有一定的參考價值。

[1] FAN T Y, BYER R L. Diode laser-pumped solid-state lasers[J]. IEEE Journal of Quantum Electronics: A Publication of the IEEE Quantum Electronics and Applications Society, 1988, 24(6): 895-912.

[2] 張輝榮，徐觀峰，李斌. Nd:CNGG晶體的生長及其性能研究[J]. 激光技術，2005(06): 599-600, 625.

[3] 李旭東. 879nm激光二極管端泵的高重頻高亮度Nd:GdVO4激光器研究[D]. 哈爾濱：哈爾濱工業大學，2010.

[4] 陳紹坤，黃海濤，陳浩. 熱效應維度對固體激光器輸出性能影響的研究[J]. 激光雜志，2016, 37(10): 14-17.

[5] 蒙裴貝，顏凡江，李旭，等. 熱邊界和泵浦結構對激光晶體熱效應的影響[J]. 紅外與激光工程，2015, 44(11): 3216-3222.

[6] 郭龍成. 二極管泵浦固體激光器中晶體熱效應的研究[D]. 西安：西安電子科技大學，2010.

[7] 高蘭蘭. LD泵浦全固體激光器幾個關鍵技術的研究[D]. 長春：中科院長春光學精密機械與物理研究所，2003.

[8] 李萌萌，楊飛，趙上龍，等. 復合Nd:YAG晶體固體激光器熱效應研究[J]. 激光與紅外，2020, 50(01): 42-48.

[9] 馮祝，萬云芳. LD端面抽運Nd:GGG激光器熱效應研究[J]. 激光技術，2014, 38(03): 360-363.

[10] 潘孫強，劉崇，趙智剛，等. 激光二極管端面抽運固體激光器的熱效應和熱透鏡焦距測量[J]. 中國激光，2010, 37(10): 2445-2450.

[11] 畢學進，周壯，董邵然，等. LD端面泵浦矩形Nd:YVO4晶體的熱效應研究[J]. 光學技術，2012, 38(06): 725-728.

[12] 施玉顯，盧鐵林，馮寶華，等. 885nm和808nm LD抽運Nd:CNGG 935nm激光器熱效應研究[J]. 中國激光，2012, 39(11): 17-21.

[13] 李萍，王垚廷，張磊. 熱透鏡效應的數值計算[J]. 齊魯工業大學學報，2017, 31(05): 62-65.

[14] 郭政鑫，高綺妮，范紫萍，等. 低功率激光熱透鏡效應的模擬及研究[J]. 大學物理，2021, 40(04): 79-85.

[15] 史彭，李隆，甘安生. LD端面抽運矩形截面YAG:Nd:YAG復合晶體熱效應[J]. 光電子激光，2006, 17(12): 1461-1465.

[16] OZYGUS B, ERHARD J. Thermal lens determination of end-pumped solid-state lasers with transverse beat frequencies[J]. Applied Physics Letters, 1995, 67(10): 1361-1362.

[17] INNOCENZI M E, YURA H T, FINCHER C L, et al. Thermal modeling of continuous-wave end-pumped solid-state lasers[J]. Applied Physics Letters, 1990, 56(19): 1831-1833.

[18] 史彭，李隆，甘安生，等. 全固態激光器中摻Nd3+離子激光晶體熱效應的研究[J]. 激光技術，2004(02): 177-180.

[19] YU G Y, WANG J Y, ZHANG H J, et al. Thermal characterization of lowly Nd3+doped disordered Nd:CNGG crystal[J]. Optics Express, 2009, 17(11): 9270-9275.

[20] KOECHNER W. Solid-state laser engineering[M]. Beijing: Science Press, 2002, 409.

[21] 艾慶康，常亮，陳檬，等. 808nm與888nm抽運Nd:YVO4熱效應分析[J]. 中國激光，2011, 38(04): 7-12.

Study on thermal effect of Nd:CNGG laser based on MATLAB

LI Qi-nan，DING Zhuang，XIANG Wan-gui，LUO Xiao-jie

(College of Science, Qiqihar University, Heilongjiang Qiqihar 161006, China)

Based on the finite difference method of MATLAB, the thermal lens effect of laser diode (LD) end pumped Nd:CNGG laser is simulated and analyzed. Under the condition of LD end pumping, the pump power, the crystal size and the position of the beam waist radius of the pump light are changed by setting the pump wavelength at 808nm and 885 nm respectively. The temperature distribution in the crystal under different conditions was obtained. The analysis shows that the above factors will affect the temperature distribution inside the crystal. For 808nm pump beam, the maximum temperature difference of the crystal is about 68K when the beam waist position is different. While for 885nm pump beam with low average temperature rise, the temperature difference is about 10K.

laser；Nd:CNGG；thermal effect；MATLAB simulation

TN248.1

A

1007-984X(2023)03-0086-04

2022-11-24

李奇楠（1975-），男，黑龍江齊齊哈爾人，教授，博士，主要從事激光光譜研究，liqinan@163.com。