切换到宽版
  • 广告投放
  • 稿件投递
  • 繁體中文
    • 2646阅读
    • 1回复

    [原创]RP Fiber Power仿真设计掺铥光纤激光器代码详解 [复制链接]

    上一主题 下一主题
    离线小火龙果
     
    发帖
    932
    光币
    2176
    光券
    0
    只看楼主 倒序阅读 楼主  发表于: 2020-05-28
    (* B;?)   
    Demo for program"RP Fiber Power": thulium-doped fiber laser, )3PQ|r'  
    pumped at 790 nm. Across-relaxation process allows for efficient YpMQY-n  
    population of theupper laser level. Q.Uyl:^PxU  
    *)            !(*  *)注释语句 A$.woE@  
    '>-  C!\t  
    diagram shown: 1,2,3,4,5  !指定输出图表 5fuOl-M0W  
    ; 1: "Powersvs. Position"     !分号是注释;光纤长度对功率的影响 ;*ebq'D([  
    ; 2:"Variation of the Pump Power"  !泵浦光功率变化对信号输出功率的影响 ?3jOE4~aHr  
    ; 3:"Variation of the Fiber Length"!信号输出功率vs 光纤长度的变化,仿真最佳光纤长度 Q`}1 B   
    ; 4:"Transverse Profiles"             !横向分布,横坐标为半径位置 v8p-<N)  
    ; 5:"Transition Cross-sections"    !不同波长的跃迁横截面,横坐标波长,纵坐标为横截面 .q#2 op  
    YFgQ!\&59  
    include"Units.inc"         !读取“Units.inc”文件中内容 R}7>*&S:  
    ]@_M)[ x  
    include"Tm-silicate.inc"    !读取光谱数据 WL(u'%5  
    jrT5Rw_}q  
    ; Basic fiberparameters:    !定义基本光纤参数 }8J77[>/  
    L_f := 4 { fiberlength }      !光纤长度 +O>1 Ed  
    No_z_steps := 50 {no steps along the fiber } !光纤步长,大括号{ }是注释,相当于备注 '^ "6EF.R  
    r_co := 6 um { coreradius }                !纤芯半径 &g;4;)p*8  
    N_Tm := 100e24 { Tmdoping concentration }  !纤芯Tm离子掺杂浓度 +ew2+2  
    _UKH1qUd4  
    ; Parameters of thechannels:                !定义光信道 6n 37R#(  
    l_p := 790 nm {pump wavelength }                !泵浦光波长790nm uCuXY#R+  
    dir_p := forward {pump direction (forward or backward) }   !前向泵浦 g &~T X  
    P_pump_in := 5 {input pump power }                    !输入泵浦功率5W -_BS!T%r  
    w_p := 50 um {radius of pump cladding }               !包层泵浦相应的半径 50um E&r*[;$  
    I_p(r) := (r <=w_p) { pump intensity profile }          !泵浦光强度分布 Ld6j;ZJ';  
    loss_p := 0 {parasitic losses of pump wave }           !泵浦光寄生损耗为0 `P}T{!P+6  
    <Okk;rj2  
    l_s := 1940 nm {signal wavelength }                   !信号光波长1940nm +~mBo+ ,  
    w_s := 7 um                          !信号光的半径 J'B6l#N  
    I_s(r) := exp(-2 *(r / w_s)^2)            !信号光的高斯强度分布 4SSq5Ve<  
    loss_s := 0                            !信号光寄生损耗为0 j;@7V4'  
    Iu`eQG  
    R_oc := 0.70 {output coupler reflectivity (right side) }      !输出耦合反射率 )>08{7  
    i)=!U>B_0  
    ; Function for defining themodel:   !定义模型函数,一定要有calc命令,否则函数只会被定义,但不会被执行 #|`/K[.xd%  
    calc `~Nd4EA)2  
      begin !Ziq^o.  
        global allow all;                   !声明全局变量 Z[:fqvXQ  
        set_fiber(L_f, No_z_steps, '');        !光纤参数 9NvV{WI-1  
        add_ring(r_co, N_Tm); 2^N 4(  
        def_ionsystem();              !光谱数据函数 qg?O+-+  
        pump := addinputchannel(P_pump_in, l_p,'I_p', loss_p, dir_p);  !泵浦光信道 d54(6N%  
        signal_fw := addinputchannel(0, l_s, 'I_s',loss_s, forward);      !前向信号光信道 zn|~{9>y  
        signal_bw := addinputchannel(0, l_s, 'I_s',loss_s, backward);    !后向信号光信道 v`8dRVN  
        set_R(signal_fw, 1, R_oc);                                 !设置反射率函数 Uv'.]#H<  
        finish_fiber();                                   u1 d{|fF  
      end; PW)XDo7  
    sxcpWSGA^  
    ; Display someoutputs in the Output window (on the right side): !在Output aera区域显示输出 Cn4o^6?"  
    show "Outputpowers:"                                   !输出字符串Output powers: qYLOq `<f  
    show"pump:     ", P_out(pump):d3:"W"  !输出字符串pump:和计算值(格式为3个有效数字,单位W) AG6tt  
    show"signal:   ",P_out(signal_fw):d3:"W" !输出字符串signal:和计算值(格式为3个有效数字,单位W) KA1Z{7UK%  
    vG9A'R'P  
    <EKDP>,~  
    ; ------------- ]5b%r;_  
    diagram 1:                   !输出图表1 ]v96Q/a  
    diN5*CF'~  
    "Powers vs.Position"          !图表名称 Mo`7YS-Y  
    zMasA  
    x: 0, L_f                      !命令x: 定义x坐标范围 {'P7D4w  
    "position infiber (m)", @x      !x轴标签;@x 指示这些字符串沿坐标轴放置 %Z?2 .)  
    y: 0, 15                      !命令y: 定义y坐标范围 >(+g:p  
    y2: 0, 100                    !命令y2: 定义第二个y坐标范围 _ Js & _d  
    frame          !frame改变坐标系的设置 }ybveZxv5A  
    legpos 600, 500  !图行在图表窗口中的位置(相对于左上角而言) xe5|pBT  
    hx             !平行于x方向网格 FaO1?.  
    hy              !平行于y方向网格 qXO@FW]  
    HH/ bBM!  
    f: P(pump, x),    !命令f: 定义函数图;P(pump, x)函数是计算x位置处的泵浦光功率 zTb!$8D"g  
      color = red,  !图形颜色 ++ !BSQ e  
      width = 3,   !width线条宽度 ((L=1]w  
      "pump"       !相应的文本字符串标签 m/l#hp+  
    f: P(signal_fw, x),  !P(signal_fw ,x) 函数是计算x位置处的前向信号光功率 7kapa59  
      color = blue,     EJ&[I%jU  
      width = 3, jeM %XI  
      "fw signal"  J5 PXmL  
    f: P(signal_bw, x),   !P(signal_bw ,x) 函数是计算x位置处的后向信号光功率 3D>syf  
      color = blue, F.ml]k&(m  
      style = fdashed, =y.!Ny5A  
      width = 3, m }I@:s2  
      "bw signal" tpp. 9  
    td{M%D,R"  
    f: 100 * n(x, 2),    !n(x ,2) 函数是计算x位置处激活粒子数在能级2上的占比 _r0[ z  
      yscale = 2,            !第二个y轴的缩放比例 y\7 -!  
      color = magenta, kx=.K'd5H  
      width = 3, 3x2*K_A5:Q  
      style = fdashed, J*s!(J |Q  
      "n2 (%, right scale)" Y(QLlJ*)/  
    U6V+jD}L]  
    f: 100 * n(x, 3),          !n(x ,3) 函数是计算x位置处激活粒子数在能级3上的占比 lrg3n[y-l  
      yscale = 2, CC,_I>t  
      color = red, $OMTk  
      width = 3, H(bR@Qok  
      style = fdashed, ^6Y4=  
      "n3 (%, right scale)" t3%[C;@wB  
    & yFS  
    ,YzrqVY  
    ; ------------- YSi[s*.G  
    diagram 2:                    !输出图表2 Z%_"-ENT  
    r}ZL{uWMW  
    "Variation ofthe Pump Power" --*Jv"/0  
    PScq-*^  
    x: 0, 10 \d~sU,L;]  
    "pump inputpower (W)", @x .9X,)^D  
    y: 0, 10 LlAMtw"  
    y2: 0, 100 .*+?]  
    frame aq/'2U 7  
    hx 0c*y~hUVZ  
    hy $|~YXH~O  
    legpos 150, 150 r9[{0y!4  
    5&V0(LT]C  
    f: (set_P_in(pump, x);P_out(signal_fw)), !set_P_in(pump,x)改变泵浦信道功率;P_out(signal_fw)输出前向信号光 3 D<s #  
      step = 5, FJ]BB4 K  
      color = blue, _ZUtQ49  
      width = 3, Qu4Bd|`(k  
      "signal output power (W, leftscale)",     !相应的文本字符串标签 U1zcJ l^  
      finish set_P_in(pump, P_pump_in) !Cse,6/Z  
    := OdjfhY  
    f: (set_P_in(pump,x); 100 * n_av(2)),   !改变泵浦信号功率对能级2上激活粒子占比的影响 ~Y=v@] 2/  
      yscale = 2, HPM ggRs  
      step = 5, w7d(|`  
      color = magenta, @&!`.Y oy  
      width = 3, Ak(_![Q:q\  
      "population of level 2 (%, rightscale)", Wp!#OY1?  
      finish set_P_in(pump, P_pump_in) CjW`cHd  
    @ 63Uk2{W>  
    f: (set_P_in(pump,x); 100 * n_av(3)),   !改变泵浦信号功率对能级3上激活粒子占比的影响 0< i]ph  
      yscale = 2, Jo%5NXts4  
      step = 5, uLok0"}  
      color = red, AC*> f&  
      width = 3, a "*DJ&  
      "population of level 3 (%, rightscale)", 09FHE/L  
      finish set_P_in(pump, P_pump_in) 1u4)  
    'F3cvpc`  
    Z x&gr|)}  
    ; ------------- A UCk]  
    diagram 3:                         !输出图表3 J=QuZwt  
    r3.A!*!  
    "Variation ofthe Fiber Length" xw2dEvjgp%  
    9='a9\((mH  
    x: 0.1, 5 ~loJYq'y  
    "fiber length(m)", @x ~fL:pVp  
    y: 0, 10 cWG>w6FI  
    "opticalpowers (W)", @y e)Q{yO  
    frame .9r+LA{  
    hx (sX=#<B%  
    hy .Z[Bz7  
    3=r#=u5z  
    f: (set_L(x);P_out(signal_fw)),     !改变光纤长度对信号光输出功率的影响 WO=X*O ne  
      step = 20,             ~M@'=Q*~  
      color = blue, $F> #1:=v<  
      width = 3, z@WuKRsi  
      "signal output" v]`}T/n  
    :)/%*<vq,  
    ;f: (set_L(x);P_out(pump)),                     !改变光纤长度对泵浦信号输出功率的影响 Vn:BasS%  
       step = 20, color = red, width = 3,"residual pump"  U~%V;*|4  
    879x(JII  
    ! set_L(L_f) {restore the original fiber length } 5v1f?btc  
    lHg&|S&J  
    )C \ %R  
    ; ------------- cM9z b6m  
    diagram 4:                                  !输出图表4 !Db 0r/_:G  
    J$Huzs#  
    "TransverseProfiles" JPeZZ13sS  
    9wR-0E )  
    I_max :=maxr(I(pump, -1, 0, 0), I(signal_fw, -1, 0, 0)) 3_%lN4sz  
    EVovx7dr  
    x: 0, 1.4 * r_co /um j Wa%vA  
    "radialposition (µm)", @x &hciv\YT2W  
    y: 0, 1.2 * I_max *cm^2 g~zz[F 8U  
    "intensity (W/ cm&sup2;)", @y qx#k()E.U  
    y2: 0, 1.3 * N_Tm >FrF"u:kM  
    frame  ;'^5$q  
    hx WD"3W)!  
    hy <p_r{  
    l| / tKW  
    f: N_dop(1, x * um,0),      !掺杂浓度的径向分布 d[h2Y/AR  
      yscale = 2, l`#XB:#U  
      color = gray, #;4afj:2g  
      width = 3, ;4E.Yr*  
      maxconnect = 1, NF/@'QRT  
      "N_dop (right scale)" f>g< :.k*  
    YX||\  
    f: I(pump, -1, x *um, 0) * cm^2,    !泵浦光沿光纤径向的强度分布 +4[L_  
      color = red, >.P/fnvJ  
      maxconnect = 1,           !限制图形区域高度,修正为100%的高度 iw`,\V&  
      width = 3, -}2'P)Xp  
      "pump" z]P =>w  
    i5oV,fiZo  
    f: I(signal_fw, -1,x * um, 0) * cm^2,  !信号光沿光纤径向的强度分布 MS`wd  
      color = blue, cmw2EHTT<  
      maxconnect = 1, H@!\?5I  
      width = 3, 43i@5F]  
      "signal" x3X^\ Ig  
    9*AH&/EXth  
    xwSi}.  
    ; ------------- v0~'`*|&  
    diagram 5:                                  !输出图表5 Y[oNg>Rz  
    RR*eq.;  
    "TransitionCross-sections" jzRfD3_s  
    fc}G6P;3{  
    I_max :=maxr(I(pump, -1, 0, 0), I(signal_fw, -1, 0, 0)) D1Yc_  
    u}zCcWP|L  
    x: 1450, 2050 +/">]QJ  
    "wavelength(nm)", @x ]_8bX}_n  
    y: 0, 0.6 :>@6\    
    "cross-sections(1e-24 m&sup2;)", @y $d<vPpJ3  
    frame 80i-)a\n  
    hx Idy{(Q  
    hy SGuR-$U`)  
    5"x=kp>!d  
    f: s12_Tm(x * nm) /1e-24,      !Tm3+吸收截面与波长的关系 iL^bf*  
      color = red, z_)`='&n  
      width = 3, XkG:1H;Q%  
      "absorption" =F!_ivV  
    f: s21_Tm(x * nm) /1e-24,  !Tm3+发射截面与波长的关系 \v7->Sy8  
      color = blue, z#$>f*b  
      width = 3, %Tc P[<  
      "emission" B%~hVpm,eM  
    r&^xg`i[z>  
     
    分享到
    离线lileisgsz
    发帖
    14
    光币
    69
    光券
    0
    只看该作者 1楼 发表于: 2021-09-28
    感谢,视频上有点看不清楚