-
高功率光纤激光器在热管理、转换效率以及成本方面相比固体激光器具有优势,光纤自身柔性传输的特性给系统集成和应用也带来了便利,在先进工业制造方面具有极其重要的应用[1-2]。目前光纤激光器领域的领头羊IPG公司基于掺镱(Yb)有源光纤的1 μm光纤激光器产品准单模单纤输出功率已超过10 kW[3],国内很多机构也开展高功率光纤激光器方面的工作[4-8],其中中国工程物理研究院和中国科学院上海光学精密机械研究所先后报道了单纤激光输出功率超过10 kW的实验结果[6-7]。受光纤光栅等器件的损伤阈值所限,高功率光纤激光器多采用主振荡功率放大器(MOPA)结构,即由振荡器产生一定功率的种子光,再经过一级或多级高功率放大器放大后得到万瓦级激光输出。由于高功率放大器泵浦功率及其提供的激光增益极高,一旦加工对象表面具有较高的反射率,其反射的激光返回放大器链路的纤芯中反向传输,将得到放大形成逐渐积累增强反向的信号光。反向信号光一方面消耗反转粒子数、引起激光器输出功率的波动和下降;另一方面其功率的积累将导致放大级输入端的泵浦吸收和受激辐射过程增强,加剧单位长度上的热负载,过高的反向光功率本身也极易造成光纤器件和泵浦源的损伤,严重时甚至导致整个放大链路的烧毁。由于光纤激光系统本身很高的输出功率,事实上很难在激光器输出端设置隔离等针对反向光的防护措施。因此,建立反向光在高功率放大器中放大过程的理论模型,分析反向光对系统运转状态的影响,对于高功率光纤激光器的设计优化具有重要的意义。
文中,笔者采用速率方程模型对高功率光纤激光系统中反向回光放大过程分别进行数值模拟,得到反向光功率放大特性分布及其对正向激光输出功率的影响,结合石英光纤的损伤阈值讨论了可能导致系统损伤的因素。
-
首先确定高功率光纤激光器的基本结构模型。表1给出了典型的万瓦级光纤激光系统相关文献报道参数,激光器的基本结构均为一级振荡器+一级放大器的MOPA结构,振荡器谐振腔由光纤光栅构成,输出功率170 W~1 kW不等,经放大器放大后达到10 kW, 振荡器和放大器均采用双包层有源光纤,包层泵浦。泵浦波长方面,IPG公司曾使用掺Yb光纤激光器的1018 nm输出对放大器进行级联泵浦实现单纤10 kW输出,与常用的976 nm半导体激光器(LD)泵浦方式相比,1018 nm级联泵浦具有更低的量子亏损,热效应相对较轻,被认为是实现单纤高功率输出的关键,但随着高功率光纤激光器研究的不断深入,研究人员发现1018 nm泵浦主要解决的是泵浦源亮度问题,即光束质量较高的1018 nm激光提高了能够耦合进入有源光纤的最大泵浦功率,随着LD以及合束器等无源器件工艺的进步,直接使用976 nm LD作为泵浦源也能够为系统提供足够的泵浦功率,其相对1018 nm级联泵浦方式略高的热负载对系统运行的稳定性也不会产生影响,后续的高功率光纤激光器相关报道中多采用976 nm LD直接对激光器进行泵浦,因此文中仿真计算中也采用976 nm LD半导体激光器作为泵浦源,端面耦合,正向包层泵浦,图1给出了掺镱光纤放大器的示意图。
表 1 典型万瓦光纤激光器系统参数
Table 1. Parameters of typical 10-kW-level fiber laser systems
Parameters IPG Photonics[3] CAEP report [7] SIOM report[7] Output power/kW 10 10.6 10.14 Structure Oscillator and one-stage amplifier Laser wavelength/nm 1 070 1 080 1 070 Seed power 1 kW 1 kW 170 W Fiber core size and lengths 30 μm,15 m 30/900 μm,25 m 30/900 μm,18 m Pump scheme 1018 nm,12690 W,
backward pumping976 nm,11500 W,
forward pumping976 nm,11359 W,
bidirectional pumping图1所示的掺镱光纤放大器中,其上能级反转粒子数密度、泵浦光、信号光随时间t和增益光纤位置z的变化规律可以用如下速率方程进行描述[8-11],此处z以有源光纤的正向起点为原点:
$$\begin{split} \dfrac{{\partial {N_2}(z,t)}}{{\partial t}} =& \dfrac{{{\lambda _{\rm{p}}}{\varGamma _{\rm{p}}}}}{{hcA}}\left[ {{\sigma _{{\rm{pa}}}}{N_1}(z,t) - {\sigma _{{\rm{pe}}}}{N_2}(z,t)} \right] \cdot \left[ {P_{\rm{p}}^ + (z,t) + P_{\rm{p}}^ - (z,t)} \right] - \dfrac{{{N_2}(z,t)}}{\tau } +\\ & \dfrac{{{\varGamma _{\rm{s}}}{\lambda _{\rm{s}}}}}{{hcA}}\left[ {{\sigma _{{\rm{sa}}}}({\lambda _{\rm{s}}}){N_1}(z,t) - {\sigma _{{\rm{se}}}}({\lambda _{\rm{s}}}){N_2}(z,t)} \right] \cdot \left[ {P_{\rm{s}}^ + (z,t,{\lambda _{\rm{s}}}) + P_{\rm{s}}^ - (z,t,{\lambda _{\rm{s}}})} \right] \\ \end{split} $$ (1) $$N = {N_1}(z,t) + {N_2}(z,t)$$ (2) $$\begin{split} \pm \dfrac{{\partial {P_{\rm{s}}}^ \pm (z,t,{\lambda _{\rm{s}}})}}{{\partial z}} + \dfrac{1}{v}\dfrac{{\partial {P_{\rm{s}}}^ \pm (z,t,{\lambda _{\rm{s}}})}}{{\partial t}} =& {\varGamma _{\rm{s}}}\left[ {{\sigma _{{\rm{se}}}}({\lambda _{\rm{s}}}){N_2}(z,t) - {\sigma _{{\rm{sa}}}}({\lambda _{\rm{s}}}){N_1}(z,t)} \right] \cdot {P_{\rm{s}}}^ \pm (z,t,{\lambda _{\rm{s}}}) - {\alpha _{\rm{s}}}{P_{\rm{s}}}^ \pm (z,t,{\lambda _{\rm{s}}}) \\ \end{split} $$ (3) $$ \pm \dfrac{{\partial {P_{\rm{p}}}^ \pm (z,t)}}{{\partial z}} + \dfrac{1}{{{v_{\rm{p}}}}}\frac{{\partial {P_{\rm{p}}}^ \pm (z,t)}}{{\partial t}} = - {\varGamma _{\rm{p}}}\left[ {{\sigma _{{\rm{pa}}}}{N_1}(z,t) - {\sigma _{{\rm{pe}}}}{N_2}(z,t)} \right] \cdot {P_{\rm{p}}}^ \pm (z,t) - {\alpha _{\rm{p}}}{P_{\rm{p}}}^ \pm (z,t)$$ (4) 式中:N为YDFA所用增益光纤的掺杂浓度;N1(z,t)为基能态粒子密度随时间t和增益光纤位置z的变化函数;N2(z,t)为上能级反转的粒子数密度随时间t、增益光纤位置z的变化函数;Γs和Γp为ASE光和泵浦光的重叠因子;A为增益光纤的模场面积;σpa和σpe分别为泵浦光的吸收和发射截面;σsa
和σse分别表示信号光的吸收和发射截面;c为真空中的光速;τ为荧光寿命;αs和αp分别为泵浦光的损耗系数;h为普朗克常数;“+”和 “−”分别为正向和反向传输;各物理量数值见表2。对连续波反向回光进行分析时,系统处于稳态,公式(1)~(4)中各物理量对时间的导数项均为零,结合放大器边界条件(5)~(6)对方程进行求解,即可获得光纤激光器中上能级粒子数密度、泵浦光功率以及正向和反向信号光功率随增益光纤位置的变化规律。边界条件(5)~(6)中,P1和P2分别表示正向和反向注入的泵浦光功率;P3和P4分别表示正向和反向注入的、中心波长为λs的信号光功率,其余波长对应功率为零。 表 2 模拟仿真所用参数
Table 2. Parameters adopted in the calculation
Parameters Value Parameters Value λp/nm Pump wavelength 976 λs/nm Signal wavelength 1080 Γp Signal overlap efficiency 0.85 Γs Signal overlap efficiency 0.85 Rcore/μm Core diameter 30 Rcladding/μm Clad diameter 900 σa (λp)/m2 Pump absorption cross section 8.5×10−25 σa(λs)/m2 Pump absorption cross section 8.5×10−26 σe(λp)/m2 Pump emission cross section 8.2×10−25 σe(λs)/m2 Pump emission cross section 9.9×10−26 h/J·s−1 Planck constant 6.626×10−34 c/m·s−1 Light speed 3×108 N/m−3 Yb3+doping concentration 1.3×1026 τ/ms Upper laser level lifetime 0.84 $$P_{\rm{p}}^ + (0) = {P_1},P_{\rm{p}}^ - (L) = {P_2}$$ (5) $${P^{\rm{ + }}}(0,{\lambda _{\rm{s}}}) = {P_3},{P^ - }(L,{\lambda _{\rm{s}}}) = {P_4}$$ (6) 参考表1中相关报道,有源光纤芯径取30 μm,另外高功率系统中为避免高泵浦吸收使热负载过于集中,应采用较小的芯包比,取包层直径900 μm。考虑高功率光纤激光系统多工作于1070~1080 nm波长,进行反向光模拟时取正向和反向信号光波长均为1080 nm,计算中代入相应的受激发射截面和受激吸收截面。计算中所用的相关参数数值如表2所示。
-
基于上一节中的速率方程模型计算光纤放大器中泵浦吸收和激光放大过程,结果如图2所示,入射泵浦功率Pp(0) 为11.8 kW,信号光功率Ps+(0) 为1 kW时,20 m长有源光纤对于入射的976 nm泵浦光吸收~13 dB,经放大后的1 080 nm信号光输出功率达到10.24 kW,斜效率83%,与相关文献报道中的仿真和实验结果相符[6-7],验证了模型的合理性,后续采用此模型进一步分析放大器中的反向回光放大过程。
图 2 光纤放大器正向信号、泵浦光功率和反转粒子数随光纤位置z向的变化(Pp(0)=11.8 kW、Ps+(0)=1 kW)
Figure 2. Forward signal, pump power and normalized population inversion of fiber amplifier as functions of fiber position z (Pp(0)=11.8 kW、Ps+(0)=1 kW)
在前述模型中,在z=20 m处引入一反向的Ps-,计算其在有源光纤中的放大情况。反向信号的来源可能是高功率光纤激光器自身输出被反射或其他高功率连续波激光器所提供,因此,其波长、发射截面和吸收截面均取与正向信号一致。图3给出固定输入泵浦功率Pp(0)=11.8 kW,正向种子光功率Ps+(0)=1 kW时,0~500 W范围内不同初始功率的反向信号光在有源光纤中的放大过程。可以看出,反向信号光在有源光纤中反向传输的过程中被逐渐放大,接近入射端处反向信号光被放大的速度更快,即反向信号的放大主要实在接近有源光纤入射端实现的,这是因为入射端吸收泵浦功率较高,反转粒子数较大,能提供更高的增益。
图 3 不同初始功率的连续波反向信号在有源光纤中的放大过程(Pp(0)=11.8 kW、Ps+(0)=1 kW)
Figure 3. Amplification of different incident backward powers in the active fiber (Pp(0)=11.8 kW、Ps+(0)=1 kW)
经过有源光纤放大后到达其输入端的反向信号光功率随初始反向信号光功率的变化关系如图4。可以看出,在反向注入功率500 W时,放大后的反向光功率可接近4 kW,如此高的激光功率经过极易造成信号泵浦合束器和振荡器光纤光栅的损伤,也有可能经过合束器进入LD泵浦源,造成泵浦源的损伤,即使在初始反向功率100 W和10 W时,放大后到达光纤输入端的反向功率也会分别达到1.05 kW和~110 W。反向信号的存在不仅会在经有源光纤放大后形成高功率反向光,成为系统损伤的隐患;另一方面,即使系统没有发生损伤,反向信号对反转粒子数的提取消耗也会导致正向输出功率的下降,图4同时给出正向信号输出功率随注入反向信号功率的变化关系,随着反向注入回光功率从0增加到500 W时,放大器的输出信号光功率从10.24 kW降至~7 kW。
图 4 不同初始功率反向信号在有源光纤中的放大结果和对输出功率的影响(Pp(0)=11.8 kW、Ps+(0)=1 kW)
Figure 4. Signal output power and amplified backward power as functions of incident backward power (Pp(0)=11.8 kW, Ps+(0)=1 kW)
根据速率方程能够进一步求得有源光纤中各处纵向单位长度上的热负载,结果如图5所示。从图中可以看出,有源光纤内部单位长度热负载最高处并不是泵浦功率最高的放大器有源光纤输入端,而是在其后一定位置。这主要是由于输入端的信号光功率并不很高,受激发射过程的速率相对较低,其对上能级反转粒子的提取效率有限,从而限制了泵浦光的吸收,即泵浦吸收饱和。随着信号光在有源光纤中得到放大,受激发射过程的强度增大,热负载也相应上升,其变化趋势与图2中剩余泵浦光曲线的变化率相符。
图 5 有源光纤内部单位长度产热负载及总热负载(Pp(0)=11.8 kW、Ps−(10)=1 kW)
Figure 5. Thermal load per unit length in the gain fiber and the total thermal load(Pp(0)=11.8 kW, Ps−(10)=1 kW)
当存在反向信号时,反向信号经放大到达放大器输出端时具有较高的功率,会增加泵浦功率最高的放大器输出端的受激吸收和受激发射过程的速率,加剧热负载。基于热负载分布的计算结果,结合热方程模型计算不同初始反向信号情况下有源光纤内部纵向温度分布趋势,计算结果如图6所示。可以看出,存在反向信号时,放大器的输入端,也即泵浦合束器与有源光纤熔接点处的温度明显上升,给激光器的稳定运转带来隐患。
图 6 不同反向功率下有源光纤内部纵向温度分布规律(Pp(0)=11.8 kW、Ps−(10)=1 kW)
Figure 6. Temperature distribution with different backward powers in the active fiber(Pp(0)=11.8 kW, Ps−(10)=1 kW)
通过上述分析可知,反向信号在有源光纤中的放大可能引起输出功率下降甚至系统的损伤,因此应尽量避免。直观上,正向的种子光功率越高,对于反转粒子数的提取越充分,在竞争中优势更明显,则反向光得到放大增益越低。为此,计算正向种子光功率Ps+(0)在500~1500 W范围内功率10 W的反向信号的放大过程,如图7。可以看出,放大器输入端处的反向回光功率会随着正向种子功率的增加而降低,而正向种子光功率降至500 W时,放大后的反向光功率超过200 W,相比正向种子光1 kW时增加了一倍以上。因此,在保证系统稳定性的和效率的前提下,选择较高的振荡器输出种子光功率有助于抑制反向回光的放大。实际工作中,振荡器输出的种子光功率的选择还需综合考虑振荡器稳定性、光纤光栅的功率承受能力、放大器有源光纤输入端热负载以及受激拉曼散射(SRS)等因素。
图 7 不同正向种子光功率时有源光纤中反向光功率的放大过程(Pp(0)=11.8 kW, Ps−(10)=10 W)
Figure 7. Amplification of backward signal in the active fiber with different forward signal powers(Pp(0)=11.8 kW, Ps−(10)=10 W)
连续波激光作用下石英光纤的实际损伤阈值鲜有实验报道,国防科技大学的研究人员通过理论计算认为30 μm芯径的光纤所能承受的最高功率为15.9 kW,进一步增加激光的功率应选择更大芯径的光纤[12]。根据前文的计算结果,连续波反向回光放大并不会使放大器有源光纤中产生相对没有反向信号时更高的功率密度,正向和方向的总功率之和仍低于损伤阈值。但是,被放大的反向光会明显增加振荡器光纤光栅处的激光功率密度,而光纤光栅的损伤阈值远低于石英光纤,目前商用高功率光纤光栅一般标称承受功率能力为2-3 kW,如光纤光栅功率处理能力的余量较小,反向光放大导致振荡器光纤光栅发生损伤的风险也较大。根据图7的计算结果可知,一旦反向回光放大导致振荡器所提供的正向种子光信号发生波动下降,反向信号在有源光纤中将得到更高的增益,到达放大器输入端的反向功率会进一步上升,对系统带来更大的损伤风险。
-
通过速率方程模型对MOPA结构的10 kW级高功率光纤激光器中存在反向回光信号时正向和方向光的功率放大特性进行数值分析,结果显示连续波反向信号会被放大器所明显放大,反向信号放大对反转粒子数的消耗会影响激光器的输出功率。考虑光纤和器件的损伤阈值,存在反向回光放大时可能导致MOPA系统中振荡器光纤光栅损伤。反向信号在放大器有源光纤中得到的放大增益受正向种子信号对反转粒子数的饱和影响,提高正向信号光功率有助于对反向信号放大的抑制,而一旦反向信号造成振荡器损伤或引起功率波动,反向信号的放大过程会被一进步加强,带来更严重的损伤风险。
Numerical analysis on backward light amplification and damage in high-power fiber laser
-
摘要: 利用速率方程模型对主振荡−功率放大器结构的1 μm波段掺镱(Yb)高功率光纤激光器中存在连续波反向信号光时的功率特性进行了理论分析,结果显示反向信号光功率会被高功率激光放大器所明显放大,10 kW级的光纤激光器中,100 W的反向信号经过放大器后功率会被放大至kW量级;与此同时,反向信号放大过程对反转粒子数的消耗会导致激光器的正向输出功率的严重下降。另外,反向信号放大也会导致放大器输出端的激光功率过强,加剧泵浦吸收和受激发射过程,增加该处的热负载、导致温度大幅上升100 ℃以上,对稳定性产生潜在影响。反向信号导致振荡器提供的正向种子光功率波动和下降时,正向信号不能充分饱和有源光纤中的增益,会进一步加强反向信号在主放大级中的放大作用,进而对系统造成更严重的影响。提高正向种子光功率、增强正向信号对激光增益的饱和作用,有助于抑制反向信号的放大过程,但需综合考虑种子源稳定性、热负载、热致模式不稳定和受激拉曼散射等因素合理选择种子光功率。Abstract: The amplification of continuous-wave backward signal in 1 μm high-power master-oscillator-power-amplifier based Yb-doped fiber laser was investigated using rate equation model. The results show that the backward light power would be amplified significantly by the high-power amplifier. The 100 W backward signal from the output end of the fiber amplifier can be amplified to up to kW level. Meanwhile, the amplification of backward signal can consume the population inversion, saturate laser gain and thus decrease the laser output power seriously. Furthermore, the backward signal amplification would result in a much higher laser intensity at the incident end of the amplifier gain fiber, where the highest pump power existed. The temperature at the incident end of the fiber can be 100 ℃ higher than that without backward signal. The higher laser intensity at the incident end could break the pump absorption saturation, and enhance the rates of pump absorption and stimulated emission a lot, hence increased the thermal load and the temperature significantly. Since the backward signal gain was determined by the saturation of population inversion by the forward seed, the power fluctuation of oscillator caused by the amplified backward signal may aggravate the backward signal amplification, and further increase the risk of damage. Higher forward seed power resulted in stronger saturation of the laser gain in the active fiber, which could suppress the backward signal amplification effectively. However, higher seed power put forwards much complex requirements to laser oscillator, and the thermal load in the active fiber of the laser amplifier would be more concentrated, which made the thermal management more difficult. Furthermore, with higher seed power, the stimulated Raman scattering and thermal induced transverse mode instability are more likely to occur. Therefore, it is important to optimize the seed laser power based on a comprehensive consideration of the above issues, and to prevent the backward light from coupling into the fiber amplifier.
-
Key words:
- fiber laser /
- fiber amplifier /
- rate equation /
- backward signal amplification /
- damage threshold
-
表 1 典型万瓦光纤激光器系统参数
Table 1. Parameters of typical 10-kW-level fiber laser systems
Parameters IPG Photonics[3] CAEP report [7] SIOM report[7] Output power/kW 10 10.6 10.14 Structure Oscillator and one-stage amplifier Laser wavelength/nm 1 070 1 080 1 070 Seed power 1 kW 1 kW 170 W Fiber core size and lengths 30 μm,15 m 30/900 μm,25 m 30/900 μm,18 m Pump scheme 1018 nm,12690 W,
backward pumping976 nm,11500 W,
forward pumping976 nm,11359 W,
bidirectional pumping表 2 模拟仿真所用参数
Table 2. Parameters adopted in the calculation
Parameters Value Parameters Value λp/nm Pump wavelength 976 λs/nm Signal wavelength 1080 Γp Signal overlap efficiency 0.85 Γs Signal overlap efficiency 0.85 Rcore/μm Core diameter 30 Rcladding/μm Clad diameter 900 σa (λp)/m2 Pump absorption cross section 8.5×10−25 σa(λs)/m2 Pump absorption cross section 8.5×10−26 σe(λp)/m2 Pump emission cross section 8.2×10−25 σe(λs)/m2 Pump emission cross section 9.9×10−26 h/J·s−1 Planck constant 6.626×10−34 c/m·s−1 Light speed 3×108 N/m−3 Yb3+doping concentration 1.3×1026 τ/ms Upper laser level lifetime 0.84 -
[1] Richardson D J, Nilsson J, Clarkson W A. High power fiber lasers: current status and future perspectives [J]. J Opt Soc Am B, 2010, 27: B63-92. doi: 10.1364/JOSAB.27.000B63 [2] Shi W, Fang Q, Zhu X, et al. Fiber lasers and their applications [J]. Appl Opt, 2014, 53: 6554-68. doi: 10.1364/AO.53.006554 [3] Stiles E. New developments in IPG fiber laser technology[C]// Proc 5th Int Workshop Fiber Lasers, 2009. [4] Fang Q, Li J, Wei Shi, et al. 5 kW near-diffraction-limited and 8 kW high-brightness monolithic continuous wave fiber lasers directly pumped by laser diodes [J]. IEEE Photon J, 2017, 9: 1506107. [5] 王标, 庞璐, 衣永青, 等. 国产 25/400 μm 掺镱光纤实现 3.2 kW 激光输出[J]. 红外与激光工程, 2019, 48: 0706009. Wang Biao, Pang Lu, Yi Yongqing, et al. 3.2 kW laser output by domestic 25/400 μm Yb-doped fiber [J]. Infrared & Laser Engineering, 2019, 48: 0706009. (in Chinese) [6] Lin H, Xu L, Li C, et al. 10.6 kW high-brightness cascade-end-pumped monolithic fiber lasers directly pumped by laser diodes in step-index large mode area double cladding fiber [J]. Results in Physics, 2019, 14: 102479. doi: 10.1016/j.rinp.2019.102479 [7] Chen X, Lou F, He Y, et al. Home-made 10 kW fiber laser with high efficiency [J]. Acta Optica Sinica, 2019, 39(3): 423-425. [8] 马建立, 姜诗琦, 于淼, 等. 基于波长锁定泵浦单振荡级千瓦光纤激光器[J]. 红外与激光工程, 2016, 45(11): 1105002. Ma Jianli, Jiang Shiqi, Yu Miao, et al. kW fiber laser oscillator based on wavelength locking [J]. Infrared & Laser Engineering, 2016, 45(11): 1105002. (in Chinese) [9] Nilsson J, Paschotta R, Caplen J E, et al . Yb3+-ring doped fiber for high-energy pulse amplification [J]. Opt Lett, 1997, 22: 1092-1094. doi: 10.1364/OL.22.001092 [10] Paschotta R, Nilsson J, Tropper A C, et al. Ytterbium-doped fiber amplifiers [J]. IEEE J Quantum Electron, 1997, 33: 1049-1056. doi: 10.1109/3.594865 [11] Wang Y, Po H. Dynamic characteristics of double-clad fiber amplifiers for high-power pulse amplification [J]. Journal of Lightwave Technology, 2003, 21(10): 2262-2270. doi: 10.1109/JLT.2003.818166 [12] Zhu J, Zhou P, Ma Y, et al. Power scaling analysis of tandem-pumped Yb-doped fiber lasers and amplifiers [J]. Optics Express, 2011, 19(19): 18645-18654. doi: 10.1364/OE.19.018645 -