HTML
-
高功率固体激光器是当前激光技术研究中的重要发展方向。为提升固体激光器的平均功率,大量泵光被注入到增益介质中,造成了有害热量积累与不良热光效应,如,热致应力双折射、热退偏等,使得光束质量、波前、功率等指标的严重下降,甚至造成晶体断裂。因此,在高功率固体激光技术中,冷却技术的研究是一项极其重要的工作[1]。高效的激光介质冷却技术能够保证其正常的工作温度和理想的温度分布,从而获得更高功率、更高光束质量的激光输出。另一方面,大功率激光介质尺寸增大与有害热累积间的矛盾,研究人员开发了多种几何形状的增益介质与冷却方案[2]。不同于常规轴对称的棒状激光晶体,板条形的增益介质具有较高的面积体积比[3-8]。
冷却换热方法包括传导换热、对流换热与辐射换热[9]。以水为基础的传导、对流冷却技术以其长期可靠、效率高、稳定换热的特点而被广泛采用。Jingchun Min[10]提出一种基于微通道热沉的新型高效传导冷却方案,用于高功率固体板条激光器,其主要参数为Re=2 000,平均对流换热系数2×104 W·m−2·K−1,然而该方案由于界面热阻的原因使换热效率无法进一步提升,此外流体方向也使得温度分布的不均匀性出现在了垂直于光束传输的方向;J. R. Wang[11]采用类薄板高速流动方案对板条晶体进行直接冷却,通过控制狭缝的高度与雷诺数,成倍增大了局部对流换热系数,然而流体方向使得晶体出现了不均匀热应变,影响光束放大的效果;Xing Fu[12]对30 kW级多板条阵列进行了数值仿真,通过对面放置的相同板条晶体组成振荡器,利用相反的流体运动方向控制晶体热应变,在信号光通过多块晶体后实现整体热补偿。
信号光在板条形激光介质中以zig-zag形式传播连续渡越前后反射面,这使得厚度方向上热梯度所带来的影响被消除。结合大面积换热,板条晶体能够在高功率注入下实现高效、高质量的激光输出。然而晶体壁面上流体温度的非定常特性与流场的方向性,传统冷却方案一般会在板条的宽度方向上引入有害热梯度,鉴于此文中基于冷却流场与目标温度匹配控制的思路,提出了一种用于双大面侧泵的大功率板条激光器的冷却结构,实现沿板条晶体长度方向的均匀流场。利用非定常边界条件的流−固耦合有限元方法对其冷却过程中的流动与传热状态进行了详细分析,对比其它传统冷却结构方案,研究了该冷却结构的工作特性及参数特性,并利用实验对仿真结果进行了验证。
-
统计模拟计算结果,如表1所示。从表1明显看出:纵向强制对流冷却方案要优于全腔直接冷却方案,且优于实际可行的微通道冷却方案,板条晶体的最高温度、平均温度、温度均匀性等各项指标均优于其他冷却方案。
Cooling configuration Flow rate/ L·min−1 Max static
pressure/ MPaMax flow
velocity/ m·s−1Max
temperature/ ℃Average
temperature/℃Temperature
varianceCavity forced convection 40.00 0.31 13.72 32.23 23.55 3.78 Micro-channel conduction 0.60 (one side) 0.87 5.95 32.30 29.00 2.05 Longitudinal forced convection 20.00 (one side) 0.22 11.33 26.1 22.69 0.95 Table 1. Simulation results of several cooling configurations
全腔直接冷却是一种较为简单的流体冷却方案,没有考虑腔体内发热元件的热分布情况,其温度分布极为不均匀,在水流入口端由于水流温度最低,此处板条晶体温度较低,在水流出口附近的水流已经被一定程度的加热,这附近的板条晶体温度相对较高,纵观水流流向,可以看出:其速度矢量方向大体为入口端指向出口端,是一个斜穿过板条晶体发热区域的方向,而在方向矢量的两侧,则因为结构而出现大量低速漩涡,因此,在此方向矢量的两侧出现了不同程度的热积累,最终形成了如图3(a)的温度场分布。这种温度分布下的板条晶体,热应力作用非常不均匀,板条晶体热畸变及热致应力折射效应较大,对波前会产生极不良的影响。
Figure 3. Temperature contour on the x-z center cross section of laser slab. (a) Cavity forced convective configuration; (b) Micro-channel conductive configuration; (c) Longitudinal forced convective configuration
在微通道传导冷却方案中,基于文中所用板条晶体,其模拟仿真结果如图3(b)所示,热积累区域集中且分布较均匀,当较小的种子光束通过板条晶体时,即通过图中红色区域,其波前在z轴方向上能量放大相对均匀。但也应该看到,如此特征的温度分布仍不理想,温度场沿着水流方向向上稍偏移了一些,这是由于流体流动的入口效应以及流动过程中流体温度升高带来的合理现象[10]。分析流道,宏通道水流在宏通道与微通道界面处会因为水阻的骤然增大而产生大量涡流,一方面消耗了水流动能,另一方面还产生了大量有害动压作用,使水流的冷却效能进一步下降,这也是微通道冷却方案必然要面对的问题。
在纵向强制对流冷却方案中,如图3(c)所示,泵浦区温度分布非常理想,其温度分布是沿x-z平面、x-y平面基本中心对称的分布状态,沿z方向的温度梯度非常小,仅在板条晶体的边缘及非泵浦区边缘存在热梯度,对于较大种子光斑,该红色热积累区域都能够完全覆盖。该方案能够很好地消除z方向上的一阶热应力问题,而在x轴方向上同样具有良好的一致性,结合zig-zag传输方式,在板条晶体中传输的激光能够被均匀放大。
-
图4为纵向强制对流冷却方案中各参数与入口流量的关系,可以看出:当流量为20 L/min时,板条晶体已经处于一个较理想的温度水平以及静压条件。继续增加水流量,流道内流速近似线性增加,水压的增加趋势也较为平缓,说明流道通路并没有达到其最大可承载的流量,同时也说明该结构还具备较大的冷却潜力;从插入的曲线来看,随着水流量增加,温度均方差随之下降,当水流量超过20 L/min时,均方差接近于1,此后曲线稳定波动于0.95附近,可知该流量已能够获得较理想的温度场分布均匀性。当然,过高的流速与静压容易引发密封失效、流体泄漏、设备短路等问题。
-
在研究充分发展流与非充分发展流两种不同流场状态时发现,两者冷却性能有一定差距,文中用非充分发展流温度场参数曲线减去充分发展流参数曲线,得到图5的非充分发展流温度场参数相对充分发展流温度场参数在同等入口流量情况下的差值曲线。红色斜线部分代表充分发展流的板条晶体平均温度与最高温度均较非充分发展流大,而绿色斜线则代表充分发展流的温度场均匀性要好于非充分发展流。
Figure 5. Relationships of the different parameter differences vs inlet flow about fully developed and developing flow
基于流体壁面阻力系数与表面对流换热系数的共性方程,推导得到公式(8):
式中:τw为壁面剪切力。
考虑到充分发展流场的特征,流动道内流场沿流动方向均匀不变,沿水流方向压力损失呈线性,则公式(8)、(9)中τw在整个流场范围内应保持不变,则对流换热系数hx也保持不变,说明在整个流场范围内流体的冷却效率具有良好的一致性;另一方面,由于充分发展流的壁面附近速度要小于非充分发展流,因此,充分发展流在壁面附近的湍流强度、湍流脉动速度等相对较小,也就使得流场整体的对流换热能力略小。
取30 L/min入口流量充分发展流与非充分发展流的模拟结果作比对,图6为两种流场在x-y中心截面上板条晶体固−液界面对流换热系数随x轴的变化曲线,黑色曲线与红色曲线为对流换热系数沿x轴分布,绿色与蓝色曲线分别为黑色与红色曲线的高阶多项式拟合结果。观察两条拟合曲线,能够看到充分发展流场的对流换热系数沿x轴变化较小,曲线较平缓。仿真结果与前述分析一致,说明充分发展流流场的换热能力在整个流场范围内较为稳定。如图6所示,在入口段的初期,由于流体速度在沿板条晶体大面方向上的速度分量较小,此时流体对板条晶体的正压力作用大于剪切作用,根据公式(8),此时对流换热系数较小,进入入口段的中后期由于水流沿板条表面速度分量基本达到最大值,加上由来流扰动造成的大量湍流,使入口段中后期的对流换热系数大大增加,并出现较大的跳跃;在距板条晶体中心±50 mm冷却区域内,对流换热系数在±30 mm的范围内分布均匀,考虑到波前畸变的影响,在冷却结构设计过程中,应增加入口段的长度,使得在冷却液换热工作区的流场处于充分发展流的状态。
-
为研究流道壁面质量在纵向强制对流冷却方案中所起的作用,以入口流量为30 L/min为基础,设流道壁面粗糙度水平(Rz:轮廓最大高度)由3e-4、1e-3、0.01、0.05、0.1等变化至30 μm,模拟结果如图7所示。当在表面粗糙度小于3.2 μm时,各参数曲线较平缓,当表面粗糙度大于3.2 μm后,板条晶体温度水平明显下降。另一方面看,随着表面粗糙度的增大,晶体对流换热表面的平均对流换热系数逐渐增大,并在大于20 μm后几乎提升了一个数量级。表面粗糙度的增大能够提升表面对流换热系数水平,一方面是由于表面粗糙度增大,使流体边界层情况复杂化,在高速主流的带动下,流体边界层内的湍流时均动能增加,边界层换热系数增加;另一方面,表面粗糙度的增大也意味着固−液接触面积的增加,降低了固−液系统的宏观界面热阻,增强了固−液界面的热传导。当然一般光学元件表面粗糙度极低,因此为提高系统对流换热系数,可在不影响泵浦光照射的情况下,在板条晶体表面加工一定形状的粗糙结构,这与Nd:YAG螺纹晶体棒的做法是一致的。
-
根据模拟结果设计装配了一款大功率板条激光放大器,水冷机出口总流量设置约为40 L/min。在短腔的储能测试过程中,待系统稳定运行,利用红外热像仪拍摄板条晶体端面温度,列举仿真模拟结果,截取晶体端面温度云图并图例归一化。
结合图8与表2数据能够发现:系统仿真结果与实际模块工作时的结果基本吻合,达到了预期冷却效能;通过测量,纵向强制对流冷却方案下的板条晶体放大器其热焦距能够达到数十米,在大功率全固态激光器系统中能够提供低热畸变与低热透镜效应的放大效能。
Figure 8. Crystal end faces temperature in investigation with laser amplifier working in longitudinal forced convection- infrared thermal imager (up) and simulation results (down)
Parameter Value Thermal power/W 850 900 960 1 020 Measurement/℃ 34.6 36.3 38.6 40.1 Simulation/ ℃ 33.51 34.38 36.44 38.17 Table 2. Crystal end faces temperature in investigation with laser amplifier working in longitudinal forced convection