RP 系列激光分析设计软件 | 示例案例：光纤中的非线性自聚焦

这里，我们研究光纤中非线性自聚焦的细节。首先，我们计算了由于非线性自聚焦的影响，大模面积光纤的基模如何收缩。

模式解算器实际上忽略了非线性效应。然而，只需几行脚本代码，我们就可以存储包括其非线性变化在内的折射率分布，然后重新计算光纤模式。重复这一过程，直到我们得到一个自洽的解:

dr := 0.05 um
defarray I[0, 200 um, dr]
n_f_nl(r) := n_f(r) + n2 * (if r <= r_max then I~[r])
  { nonlinear refractive index profile }
store_I(P) := 
  for r := 0 to 2 * r_co step dr do
    I[r] := P * I_lm(0, 1, lambda, r)
    { ignore index changes outside 2 * r_co, where the intensity is small }

CalcNonlinearMode(P) := 
  { Calculate the lowest-order mode with self-focusing for the power P. }
  begin
    var A, A_l;
    A := 0;
    repeat
      A_l := A;
      store_I(P);
      set_n_profile("n_f_nl", r_max);
      A := A_eff_lm(0, 1, lambda);
    until abs(A_l / A - 1) < 1e-6;
  end

考虑到光纤的非线性，可以对光束的传播进行数值模拟。为此，我们需要定义一个数值网格，并为光束传播设置各种其他输入:

x_max := 30 um { maximum x or y value }
N := 2^5 { number of grid points in x and y direction }
dx := 2 * x_max / N { transverse resolution }
z_max := 30 mm { fiber length }
dz := 100 um { longitudinal resolution }
N_z := z_max / dz { number of z steps }
N_s := 100 { number of sub-steps per dz step }

P_11 := 4 MW
A0%(x, y) := sqrt(P_11) * A_lm_xy(1, 1, lambda, x, y)  { initial field }

calc
  begin
    bp_set_grid(x_max, N, x_max, N, z_max, N_z, N_s);
    bp_define_channel(lambda);
    bp_set_n('n_f(sqrt(x^2 + y^2))'); { index profile }
    bp_set_loss('10e2 * ((x^2 + y^2) / (20 um)^2)^3');  { simulate loss for cladding modes }
    bp_set_n2('n2');
    bp_set_A0('A0%(x, y)'); { initial amplitude }
    bp_set_interpol(2); { quadratic interpolation }
  end