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

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

    上一主题 下一主题
    离线小火龙果
     
    发帖
    932
    光币
    2176
    光券
    0
    只看楼主 倒序阅读 楼主  发表于: 2020-05-28
    (* Iq[ d5)M4  
    Demo for program"RP Fiber Power": thulium-doped fiber laser, q-3e^-S*  
    pumped at 790 nm. Across-relaxation process allows for efficient 3d>3f3D8;  
    population of theupper laser level. U WU PY  
    *)            !(*  *)注释语句 mu>L9Z~(L_  
    !&f>,?wlP  
    diagram shown: 1,2,3,4,5  !指定输出图表 O-N@HZC  
    ; 1: "Powersvs. Position"     !分号是注释;光纤长度对功率的影响 &^I2NpT  
    ; 2:"Variation of the Pump Power"  !泵浦光功率变化对信号输出功率的影响 `{B<|W$=  
    ; 3:"Variation of the Fiber Length"!信号输出功率vs 光纤长度的变化,仿真最佳光纤长度 zz7#g U  
    ; 4:"Transverse Profiles"             !横向分布,横坐标为半径位置 7rZE7+%]  
    ; 5:"Transition Cross-sections"    !不同波长的跃迁横截面,横坐标波长,纵坐标为横截面 VGVb3@  
    D-S"?aO-  
    include"Units.inc"         !读取“Units.inc”文件中内容 :&'[#%h8  
    y.6Yl**l  
    include"Tm-silicate.inc"    !读取光谱数据 w(EUe4 w{  
    UWPzRk#s"  
    ; Basic fiberparameters:    !定义基本光纤参数 ?j:g.a+U  
    L_f := 4 { fiberlength }      !光纤长度 q=J8SvSRl  
    No_z_steps := 50 {no steps along the fiber } !光纤步长,大括号{ }是注释,相当于备注  (%\tE  
    r_co := 6 um { coreradius }                !纤芯半径 {Rv0@)P$  
    N_Tm := 100e24 { Tmdoping concentration }  !纤芯Tm离子掺杂浓度 .!^}sp,E  
    's>./Pf  
    ; Parameters of thechannels:                !定义光信道 qJyGr ?  
    l_p := 790 nm {pump wavelength }                !泵浦光波长790nm Q$B\)9`v[  
    dir_p := forward {pump direction (forward or backward) }   !前向泵浦 6$y$ VeW  
    P_pump_in := 5 {input pump power }                    !输入泵浦功率5W b;~?a#Z}  
    w_p := 50 um {radius of pump cladding }               !包层泵浦相应的半径 50um i uGly~  
    I_p(r) := (r <=w_p) { pump intensity profile }          !泵浦光强度分布 .271at#-  
    loss_p := 0 {parasitic losses of pump wave }           !泵浦光寄生损耗为0 Tg;1;XM%  
    g4U`Qf3  
    l_s := 1940 nm {signal wavelength }                   !信号光波长1940nm LV6BSQyQ  
    w_s := 7 um                          !信号光的半径 d"#& VlKcv  
    I_s(r) := exp(-2 *(r / w_s)^2)            !信号光的高斯强度分布 W02t6DW  
    loss_s := 0                            !信号光寄生损耗为0 a?NoNv)&  
    J^xIfV~ zt  
    R_oc := 0.70 {output coupler reflectivity (right side) }      !输出耦合反射率 Frd`u .I  
    8IQqDEY^  
    ; Function for defining themodel:   !定义模型函数,一定要有calc命令,否则函数只会被定义,但不会被执行 q Xj]O3 mm  
    calc 'X(Sn3  
      begin Tsc2;I  
        global allow all;                   !声明全局变量 Ae3=o8p  
        set_fiber(L_f, No_z_steps, '');        !光纤参数 DFvj  
        add_ring(r_co, N_Tm); L_(Y[!  
        def_ionsystem();              !光谱数据函数 $Ao iH{f  
        pump := addinputchannel(P_pump_in, l_p,'I_p', loss_p, dir_p);  !泵浦光信道 11Y4oS  
        signal_fw := addinputchannel(0, l_s, 'I_s',loss_s, forward);      !前向信号光信道 hha!uD~(  
        signal_bw := addinputchannel(0, l_s, 'I_s',loss_s, backward);    !后向信号光信道 37U$9]  
        set_R(signal_fw, 1, R_oc);                                 !设置反射率函数 TMAart; <  
        finish_fiber();                                   :3 p&h[M  
      end; U+7!Vpq  
    FrL ;1zt  
    ; Display someoutputs in the Output window (on the right side): !在Output aera区域显示输出 +h?Rb3=S  
    show "Outputpowers:"                                   !输出字符串Output powers: AY)R2> fW%  
    show"pump:     ", P_out(pump):d3:"W"  !输出字符串pump:和计算值(格式为3个有效数字,单位W) N&YQZ^o  
    show"signal:   ",P_out(signal_fw):d3:"W" !输出字符串signal:和计算值(格式为3个有效数字,单位W) dxk~  
    i^_?C5  
    dkI(&/  
    ; ------------- ^sb+|b  
    diagram 1:                   !输出图表1 -D^.I  
    UkzLUok]U  
    "Powers vs.Position"          !图表名称 _2p D  
    #Ab,h#f*7  
    x: 0, L_f                      !命令x: 定义x坐标范围 =+>^:3cCQ  
    "position infiber (m)", @x      !x轴标签;@x 指示这些字符串沿坐标轴放置 1_RN*M +#  
    y: 0, 15                      !命令y: 定义y坐标范围 XMi)PXs$  
    y2: 0, 100                    !命令y2: 定义第二个y坐标范围 yh{Wuz=T  
    frame          !frame改变坐标系的设置 5 cz6\A&  
    legpos 600, 500  !图行在图表窗口中的位置(相对于左上角而言) s{@3G8  
    hx             !平行于x方向网格 bG&vCH;}%  
    hy              !平行于y方向网格 T.B} k`$  
    M7. fz"M  
    f: P(pump, x),    !命令f: 定义函数图;P(pump, x)函数是计算x位置处的泵浦光功率 b+!I_g4P  
      color = red,  !图形颜色 LvbS")  
      width = 3,   !width线条宽度 6SVh6o@]  
      "pump"       !相应的文本字符串标签 oV,lEXz  
    f: P(signal_fw, x),  !P(signal_fw ,x) 函数是计算x位置处的前向信号光功率 68YJ@(iS  
      color = blue,     }&qr"z4  
      width = 3, D4;6}gRC  
      "fw signal" 5nh:S0M6V  
    f: P(signal_bw, x),   !P(signal_bw ,x) 函数是计算x位置处的后向信号光功率 OdHl)"#  
      color = blue, m[E#$JZtG  
      style = fdashed, eW.[M?,  
      width = 3, }v$T1Cw  
      "bw signal" r>.^4Z@  
    ]~J.YX9ST  
    f: 100 * n(x, 2),    !n(x ,2) 函数是计算x位置处激活粒子数在能级2上的占比 ^eHf'^Cvvu  
      yscale = 2,            !第二个y轴的缩放比例 ,W:Bh$%  
      color = magenta, }\wTV*n`X  
      width = 3, n1+,Pe*)  
      style = fdashed, jSMs<ox  
      "n2 (%, right scale)" 3E`poE  
    y jQpdO  
    f: 100 * n(x, 3),          !n(x ,3) 函数是计算x位置处激活粒子数在能级3上的占比 = }6l.9  
      yscale = 2, 81&5g'  
      color = red,  EWn\ ]f|  
      width = 3, m~U2 L  
      style = fdashed, XJ9l, :c,  
      "n3 (%, right scale)" [/Ya4=C@  
    w$)E#|i  
    GFmVR2z_+  
    ; ------------- 3GU JlFj  
    diagram 2:                    !输出图表2 xK;WJm"  
    L7 f'  
    "Variation ofthe Pump Power" 1h0cId8d  
    7\p<k/TS  
    x: 0, 10 @o6^"  
    "pump inputpower (W)", @x 7.DAwx.HYK  
    y: 0, 10 q)E J?-  
    y2: 0, 100 8jxs%N,aI  
    frame SYZS@o  
    hx c N^,-~U  
    hy hp6%zUR  
    legpos 150, 150 .f~x*@  
    ;.wWw" )  
    f: (set_P_in(pump, x);P_out(signal_fw)), !set_P_in(pump,x)改变泵浦信道功率;P_out(signal_fw)输出前向信号光 IsjD-t  
      step = 5, ~.'NG? %7P  
      color = blue, #w<:H1,4  
      width = 3, q9`!T4,  
      "signal output power (W, leftscale)",     !相应的文本字符串标签 \JX.)&> -  
      finish set_P_in(pump, P_pump_in) ob3Z I  
    kH10z~(e  
    f: (set_P_in(pump,x); 100 * n_av(2)),   !改变泵浦信号功率对能级2上激活粒子占比的影响 \%ZF<sV W  
      yscale = 2, 9azk(OL6  
      step = 5, SOPQg?'n=V  
      color = magenta, r\sQ8/  
      width = 3, Ikbz3]F^V  
      "population of level 2 (%, rightscale)", '5vgpmn  
      finish set_P_in(pump, P_pump_in) kb>/R/,9  
    DTw3$:  
    f: (set_P_in(pump,x); 100 * n_av(3)),   !改变泵浦信号功率对能级3上激活粒子占比的影响 Gj}P6V _  
      yscale = 2, L8zY?v(bG  
      step = 5, ]/y&5X  
      color = red, ;2=H7dq  
      width = 3, 1*!`G5c,}  
      "population of level 3 (%, rightscale)", Uhz<B #tj  
      finish set_P_in(pump, P_pump_in) E'J| p7  
    `1$7. ydQ  
    Wi?37EHr  
    ; ------------- BEv>?T 0  
    diagram 3:                         !输出图表3 'nFqq:2Xa  
    YLfZ;W|6u  
    "Variation ofthe Fiber Length" k^IC"p Uc  
    6k=ink-/  
    x: 0.1, 5 v!pT!(h4  
    "fiber length(m)", @x ~Z'3(n*9  
    y: 0, 10 PB :Lj  
    "opticalpowers (W)", @y ~X/1%  
    frame . d;XLS~  
    hx IaU  
    hy 7xOrG],E  
    Rj|8l K;,  
    f: (set_L(x);P_out(signal_fw)),     !改变光纤长度对信号光输出功率的影响 @ qfVt  
      step = 20,             yBPaGZ{f  
      color = blue, 45hjN6   
      width = 3, ~ZSP K;D[  
      "signal output" $Qv+*%c  
    9W{=6D86e  
    ;f: (set_L(x);P_out(pump)),                     !改变光纤长度对泵浦信号输出功率的影响 x"Hi!h)v  
       step = 20, color = red, width = 3,"residual pump" L.[ H   
    f@R j;R~Jp  
    ! set_L(L_f) {restore the original fiber length } I]]3=?Y  
    FX FTf2*T  
    J/j?;qx]j  
    ; ------------- T>&d/$;]  
    diagram 4:                                  !输出图表4 1r?<1vh:z  
    Fvy__ qcHi  
    "TransverseProfiles" D! 1oYr  
    O6^>L0'  
    I_max :=maxr(I(pump, -1, 0, 0), I(signal_fw, -1, 0, 0)) T #&9|  
    t&SC>8M<  
    x: 0, 1.4 * r_co /um X;7gh>Q'4  
    "radialposition (µm)", @x 1Z +3=$P  
    y: 0, 1.2 * I_max *cm^2 *N .f_s  
    "intensity (W/ cm&sup2;)", @y 8"4&IX  
    y2: 0, 1.3 * N_Tm n# %mL<  
    frame gsn3]^X  
    hx $/aZ/O)F  
    hy Y{TzN%|LV  
    X=Q)R1~6v  
    f: N_dop(1, x * um,0),      !掺杂浓度的径向分布 F#X&Tb{  
      yscale = 2, 4+od N.  
      color = gray, coHzbD~#H  
      width = 3, +s:!\(BM  
      maxconnect = 1, "r!O9X6  
      "N_dop (right scale)" ; /fZh:V2  
    ,%#FK|  
    f: I(pump, -1, x *um, 0) * cm^2,    !泵浦光沿光纤径向的强度分布 ] ZGP  
      color = red, 1nb]~{l  
      maxconnect = 1,           !限制图形区域高度,修正为100%的高度 7u!i)<pn  
      width = 3, 8F(lW)An  
      "pump" }&)X4=  
    bsw0+UY=9  
    f: I(signal_fw, -1,x * um, 0) * cm^2,  !信号光沿光纤径向的强度分布 0*MY4r|-  
      color = blue, Fu0 dYN  
      maxconnect = 1, 5Xu2MY=  
      width = 3, %vvA'WG  
      "signal" $DZ\61  
    \0iF <0oy  
    a$p?r3y  
    ; ------------- IWvLt  
    diagram 5:                                  !输出图表5 D9M<>Xz)  
    V,<3uQD9a  
    "TransitionCross-sections" lC Bb0k2  
    D.zEE-cGyb  
    I_max :=maxr(I(pump, -1, 0, 0), I(signal_fw, -1, 0, 0)) 5w)tsGX\  
    GndU}[0J  
    x: 1450, 2050 a9C8Q l  
    "wavelength(nm)", @x gwDQ@  
    y: 0, 0.6 Z,%^BAJ  
    "cross-sections(1e-24 m&sup2;)", @y D<5;4Mb  
    frame 5F{NPKa Q  
    hx I@a y&NNh  
    hy nqNL[w6{  
    j:# wt70  
    f: s12_Tm(x * nm) /1e-24,      !Tm3+吸收截面与波长的关系 +:3K?G -  
      color = red, o(GXv3L  
      width = 3, nFU'DZ  
      "absorption" JsohhkJNGi  
    f: s21_Tm(x * nm) /1e-24,  !Tm3+发射截面与波长的关系 3-z; pk  
      color = blue, {3F;:%$`c  
      width = 3, 3 f=_F  
      "emission" z^z_!@7v   
    w Q /IT}-  
     
    分享到
    离线lileisgsz
    发帖
    14
    光币
    69
    光券
    0
    只看该作者 1楼 发表于: 2021-09-28
    感谢,视频上有点看不清楚