指南1:如何计算一个简单激光腔中端面泵浦棒的热透镜 i_LF`JhEQT 目录 |"$uRV=qm 1. 启动LASCAD和定义一个简单激光腔 1 68'>Zbelb 2. 定义和分析端泵浦的棒 2 do>"[RO 2.1. 选择晶体的类型和泵浦结构 3 W;Jx<-#1 2.2. 定义泵浦光的分布 3 .`Zf}[5[ 2.3. 定义棒的冷却 6 =KX<_;E 2.4. 定义工作物质参数 6 cQZ652F9 2.5. 定义合成物质 7 n1:v HBM@\ 2.6. 定义控制FEA计算程序的选项 7 g#K'6VK{ 2.7. FEA结果的可视化 9 1(GHCxA8G 2.7.1. 3D观察器 9 F X1ZG! 2.7.2. 2D数据分布和抛物线拟合 10 B7-RU<n 2.8. 高斯模的计算 10 Q\J,}1<`6 2.9. 在模式图中插入晶体 11 Y&r]lD 3. 调整腔参数 12 E)t 4. 激光腔特性分析工具 14 m^\TUj 4.1. 激光腔稳定性分析 14 {- 63/z 4.2. 高斯横模的分布以及与泵浦光分布的重叠显示 14 IdRdW{o 4.3. 计算激光功率输出 15 K,}w]b 5. 光束传输代码(BPM) 16 Y5!b)vke EZ(^~k=I 1. 启动LASCAD和定义一个简单激光腔 fRg=!<#%
?Zyok]s 选择 Start/Programs/LASCAD/lascad启动LASCAD,定义工作路径,单击“OK”打开LASCAD 主窗口,单击工具栏最左端的“New Project”按钮或者主菜单上的“File”选项,将“Number of Face Elements”设置为4,输入适当的波长,其余选项保持默认设置,单击“OK”。 dyFKxn`,
a~^Srj!}x
{rkn q_;0 图1
U06o;s(
c! H 9yk
现在可以看到LADCAD 的主菜单出现在窗口的上方,另外两个窗口在下方,其中一个窗口的标题是“Standing Wave Resonator”,另外一个的标题是“Parameter Field”,如图1所示。 3On
JWuVfZ
R+$8w2#
上方的窗口是有4个元件的简单腔的模式图,下方窗口显示这个腔的参数。在每一个元件的下面的一列里显示这个元件的参数,比如在标有“Type-Param”的一例里显示每个镜的曲率半径。可以用标有元件序号的下拉框按钮来改变相应元件的类型,可以选择Mirror、Dielectric Interface或者Lens类型。每一元件序号下面的列显示相应的两元件的空间属性,比如折射率、“Refractive Parameter”(表示相应的抛物线型折射率分布的二阶导数)。关于这个窗口的其他功能的描述可以在第3节、快速入门或者使用手册中找到。 5y_"
s"0b%0?A
2. 定义和分析端泵浦的棒 ^4`&EF
l@d
gJ %XXkVK` 图2
f@:CyB GQ
{B yn{?w
2.1. 选择晶体的类型和泵浦结构 {.#zHL
;
3BMS_,P
单击LASCAD 主窗口中的“FEA/Parameter Input & Start of FEA Code”打开标题为“Crystal, Pump Beam and Materials Parameters”的窗口,如图2。注意有六个标签页来定义不同的参数组。 n'64;J5
在“Models”标签的对话框中列举了不同的晶体和泵浦结构供选择。默认的设置是“(Dual) end pumped cylindrical rod”。这个窗口的底部可以定义棒的尺寸,比如我们输入棒的长度为6mm 和直径为2mm,如图2所示。 (@ea|Fd#4
J/4y|8T/y
2.2. 定义泵浦光的分布 Plz-7fy33
5&4F,v[zp
选择“Pump Light”标签,打开如图3所示的窗口。 ~2N-k1'-'
.~dEUt/|) 图3
\kwe51MQ
5(}H
?
在对话框中可以选择超高斯函数近似泵浦光的分布。超高斯函数是指光的分布具有高斯分布的数学形式,但其指数比2大。增加超高斯函数的指数可以形成平顶分布。关于这些内容的详细描述请参阅LASCAD手册中。 hBCR]=']
iAz UaF
在“Incident pump power, W”标签的对话框中可以分别定义晶体左右端面的泵浦功率。对端泵而言,通常泵浦光束聚焦到棒中。在"Spot size at beam waist", "Distance of waist from left end face" 和"Beam divergence"的对话框中可以定义相关量的大小。 iiwpSGFl]
单击“From the right end”按钮就可以输入从右端输入泵浦光的参数。 PF ;YE6
~{
.,8jE
通过下列操作可以看到超高斯指数效应: Rx.
rj~
PbUI!Xqe`
选定“x-y-symmetry”选项,打开“Super-gaussian exponent”选项y平面的下拉框,选择其参数比如10,单击底部右端的标有“Show Pump Light Distribution”的按钮,其泵浦光分布窗口将打开,如图4所示。 S.fXHtSx
c57b f A. Nz_! 图4
9c5!\m1 图5
(/q}mB 输入超高斯指数,光分布沿x轴是通常的高斯分布,而y轴呈平顶分布。从左至右图形拖动窗口下方的滑动条,能看到由于吸收而导致的吸收泵浦能量强度的减少。单击图形可以得到最大高度的图形。
;2W2MZ!TF
g@(30{ 5~yb
~0 2.3. 定义棒的冷却 0e9A+&r 选定“Boundaries”标签打开如图5所示的窗口。
T^.{9F]*S 在对话框中可以对棒的每个表面分别定义冷却条件,可以选择用固体冷却也可以选择用液体。对于液体冷却需要选定 “Fluid Cooling”选项。
Z)v)\l9d 对于固体冷却情况,其保持常温的表面温度可以在“Temperature, K”对话框中定义。而对于液体冷却情况,定义液体的体积温度,另外,在底部的膜系数(它描述与界面处的固体或液体的热传导)对话框必须定义。更详细的描述请参阅LASCAD手册6.10.3节。
zlX!xqHj 对话框“Reference temperature”可以用来计算畸变和相应的晶体在受热前的温度。当边界温度是用开氏温度定义时,注意输入正确的温度值。
WRMz]|+}4 lIy/;hIc 图6
HU3Vv<lz 7dhn'TW 2.4. 定义工作物质参数选定“Material Param”标签打开如图6所示的窗口。
L9'- “absorption coefficient”描述由于泵浦光的吸收,泵浦光按 呈指数衰减,它也与掺杂浓度有关系。各向异性的晶体对泵浦光的吸收还与泵浦光的偏振态有关,这种情况下的吸收因子可以定义为两个独立的偏振方向。更详细的描述请参阅LASCAD 手册的附录部分。对Nd:YAG晶体,只需要输入相等的值到相应的对话框。另外,更复杂的内容(比如,温度对物质参数的依赖,由两种不同物质组成的固体物质的情况)请参阅使用手册。
>#i $Tw 还可以将物质参数保存为一个文件便于以后在新项目中使用这些参数。2.5. 定义合成物质
kOM- 选择“Doping &Mats”标签打开如图7所示的窗口。在现在这个例子中,我们用输入框定义端盖未掺杂的棒的泵浦区域范围[mm]from z=1.0 to z=5.0意思是说,泵浦光只有在1 ≤ z ≤ 5 mm这个掺杂区域才被吸收。由于我们例子里使用的晶体是6mm长,我们定义两非掺杂端盖为1mm长。
y84XoDQ 图7
?lG;,,jc,W bG1 ofsU 这个窗口显示的参数与在“Material Param”窗口的选择有关。
;G$)MS'nB vcD'~)G(* 2.6. 定义控制FEA计算程序的选项 &.*T\3UO e6es0D[>5 选定“FEA Options”标签打开如图8所示的窗口。
fsb=8>}63} $UmE 有限元代码(FEA code)使用晶体中一个规则长方形网格来计算,这个长方形与它所相邻的那些小的不规则的网格有关。
1(F'~i|5 “Mesh size in x- and y-direction”是指与晶轴垂直的体积元的边界长。
Pb=rFas*C “Mesh size in z-direction”指与晶轴平行的体积元的边长。
pgfu+K7?w 单击“Estimated number of elements”按钮可以得到进行计算所需的随机存储内存(RAM)的大小。可以按照下面的规则计算所需RAM:如果有n个元素,则至少需要2n/1000MB RAM,否则计算机会因为内存与硬盘的数据交换速度太慢而使得计算速度很慢。(你可以这样在任务管理窗口来确定当前空RAM的大小:同时按下WinNt/2K/XP操作
系统的Ctrl-Alt-Del 三键,从弹出的窗口中选定任务管理按钮,打开“System Performance”窗口,在物理内存窗口(Physical Memory)显示了以K计算的当前RAM。
H.;yLL= 对话框“Convergence limits”用来控制迭代计算程序的收敛性。对热分析情况,默认值为1.0E-7时如果最高温度的前7位没有变化,代码将停止运行。对结构分析的限制是最大节点位移的绝对温度值。如果最大迭代次数超过“Maximum number of iterations”对话窗口中设置的次数,不论其收敛性如何,迭代程序将停止。在对话框“Directory for output of FEA results”中可以定义生成的FEA代码文件的保存路径,默认的路径是工作路径的子路径LascadFEA。
R mgxf/ 在对话框“Position of cutting plane perpendicular to z-axis”中可以定义放置一个垂直于z轴的切平面并以3D可视化方式显示晶体中物理量的分布,下一节会具体描述3D可视化的问题。如果输入的值为0,切平面将沿着棒的轴线和垂直于y轴方向。如果输入的值大于0但小于棒长,切平面垂直于棒的轴线。如果输入值大于棒的长度,整个晶体都会显示出来。
9w$7VW; 按下“Apply and run FEA”按钮可以将输入的数据传送到内部变量后运行FEA代码。弹出的“Finite Element Analysis”窗口显示当前迭代计算的次数。此外,它还在热分析时显示最高温度,在结构分析时显示最大节点位移的值。计算完成后将出现“FEA finished successfully”消息框,单击“ok”按钮关闭对话框。
;op'V6iG 需要注意的是,对于元素个数较多的情况,FEA的初始化和网格的生成需要一段时间。
V?WMj
$l< 6A@Lj*:2m zrTY1Asw;4 图8
|<2JQ[] ^8U6"O6|X 图9
kOzt"t&